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PEEFAOE. 

In compiling this work I have freely drawn on all best 
known books of such authorities as Appell, Routh, Loney, 
Bouasse and many others. My object is to present the 
matter in an elementary form and to make it immediately 
intelligible to a reader possessing but a limited knowledge of 
mathematics. For this reason it was thought necessary to 
explain things, already almost evident; constantly to refer 
back to elementary books on mechanics, etc.; and to devote 
the whole first chapter to elements, which are no doubt 
already familiar to the average reader, but which he may 
find presented in a somewhat different form from that in 
which he has them fixed in his mind. 

Lagrange's method is like a slide rule: it has its limitations, 
yet, in many problems it enables us to write down the differ- 
ential equations of motion almost instantly. 

The wonderful beauty and power of this method will un- 
doubtedly appeal to the reader, engineer or student, and 
make him like the whole subject of dynamics, although his 
teachers may have completely failed even to interest him in it, 
as often is the case, beyond the painful necessity of memorizing 
a few distorted notions. 

However the primary object of the book is to be used in 
everyday practice; the writer, being but an average engineer, 
uses this method to great advantage in working out various 
problems of construction, etc. Why not suppose that others 
might likewise derive some benefit from this brief exposition 
of its principles? Those who want to know more are referred 
to Appell, MIcanique Rationnelle, Vol. II, and Routh, Dy- 
namics of Rigid Bodies, Vols. I and II. 

My thanks are due to Dr. Eric Doolittle, Director of the 

Flower Astronomical Observatory, for reading the MSS. and 

making many valuable suggestions. 

N. W. A. 

Philadelphia, 

December 5, 1916. 

111 

308816 



i 



Chapter I. 
Chapter II. 
Chapter III. 
Chapter IV. 
Chapter V. 



CONTENTS. 

page. 

Synopsis of principles of dynamics 1 

Lagrange's equations for a particle 86 

Lagrange's equations for a system 113 

Lagrange's equations for relative motion 144 

Small oscillations 160 



CHAPTER I. 

Brief Synopsis of Certain Principles of Dynamics. 

1. Constraints. By constraints special conditions are meant, 
limiting the motion of a particle in a certain manner, pre- 
scribed beforehand. For instance a particle may be free to 
move only along a certain curve (a small ring sliding on a 
curved wire, a car on the track, etc.); or, again, the particle 
may be compelled to remain, at all times, in contact with a 
certain surface (for instance if connected by means of a rod 
to a fixed point, about which it can, therefore, move on a 
sphere). Very often, the particle can move only on the 
exterior of a certain surface (imagine, for instance, a small 
particle sliding off a circular log, a well-known problem); 
or the distance between a particle and a certain fixed point 
may be prescribed to be equal to or less than a certain value 
(case of a stone on a string). All these are typical instances 
of constrained motion. 

Analytically, the constraints are specified by geometric 
equations. For instance, the surface on which the particle 
is comipelled to remain is usually given by some such equation 
as f{x, y, z) = 0; the curve, along which the particle can 
slide, would be given as the intersection of two surfaces such 
as fix, y, z) = 0, and F{x, y, z) = 0. 

It is of the utmost importance to note that the constraints 
may be either permanent or movable; that is, changing their 
position or even their shape. Consider, for instance, the 
motion of a particle constrained to move in a plane which 
itself is rotating, say, about a vertical axis with a certain 
angular velocity; or, the motion of a small particle of dust 
upon a soap bubble while it is being inflated. Conditions of 
this sort are characterized by the fact that the equations of 
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2 Lagrange's Equations. 

constraints contain the time; that is, the equation of a surface 
would then be f(x, y, z, t) = 0; and the constraining curve 
would be given by such equations as f{x, y, z,t) = 0; and 
F(x, y, 2, t) = 0; so that the derivative with respect to time 
would then not equal 0. If it is 0, this means that the con- 
straints are permanent, or independent of time. 

In the absence of constraining conditions the motion of a 
particle is termed free. 

2. Virtual work. The fundamental conception of virtual 
work and virtual velocity is known from elementary treatises 
(Bowser, Anal. Mech., p. 166). By virtual displacement 
we shall understand a very small displacement of a particle, 
conceived or imagined by us to take place in any direction 
whatsoever; it may or may not coincide with the displace- 
ment actually taking place under the action of the given forces 
and other conditions; the latter is called actital displacement. 
In case of constrained motion, certain displacements, called 
compatible or consistent with the constraints, can be con- 
ceived. For instance, in the case of a constraining curve the 
only compatible displacement would be either backward or 
forward, from some initial position, along the curve; in the 
case of a constraining surface, compatible displacements of a 
particle can be imagined to take place in a great variety of 
manners, but always subject to the initial condition, viz., 
adhesion to the surface. Other displacements cannot even 
be conceived without calling into play the idea of distorting 
the constraints; they are called inconsistent with the con- 
straints and will not here be considered; while under free or 
unconstrained motion the virtual displacements may be any. 

For the sake of clearness let us write down the few funda- 
mental principles and definitions established so far: (a) By 
virtual work of a force is meant the product of the virtual 
displacement of its point of application into the projection 
of the force upon the direction of the displacement, in other 
words the virtual work = P-dp' cos (P, Sp) ; (6) The virtual 
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work of a force for any displacement is equal to the sum of 
virtual works done by its components; in other words, 

R'Ss- cos (R, Ss) = SP-Sp-cos {P, dp); 

(c) For concurring forces the sum of virtual works done by 
the forces equals the virtual work done by their resultant; 
from this is derived the very important form in which virtual 
work is given in rectangular coordinates (x, y, z). Supposing 
that there is a force P referred to rectangular axes Xj y, z 
and that the projections of the force upon these axes are (say) 
Xo, Yq and Zq. Let the virtual displacement of the force 
be Sp, of which the projections upon the axes will be Sx, 
by, iz. Now in view of what has just been said, the virtual 
work of the force must equal the sum of virtual works of its 
components; that is, 

P-6p' cos (P, 8p) = Xobx + Yoby + ZqSz; 

we will represent this simply by bW\ (d) From (c) it also 
follows that when any number of concurring forces are in 
equilibrium the sum of their virtual works is = 0. {e) To 
the above the following principle should be added : In case of 
rotation the virtual work is the product of the moment of the 
force about the axis of rotation by the angular (virtual) 
displacement. 

All of this refers to free motion. So far as constrained 
motion is concerned the following remarks may be made: 
If the motion of a particle is constrained, this of course means 
that at any time the coordinates of the particle must satisfy 
the constraining equation 

f(x, y,z) = 0; or, f{x, y, z, t) = 0,. (1) 

otherwise the particle would not remain on the constraining 
curve or surface; and if a small virtual displacement, of which 
the projections upon the axes are bx, by, bz, be given to the 
particle, compatible with the constraints, then the new posi- 
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tion of the particle will still satisfy the equations (1). That 
is, if the original coordinates of a particle were x, y, z, satis- 
fying the constraining condition (1), then the new coordinates, 
X+ Sx, y + Sy, z + Sz, must also satisfy the condition (1) 
for the same instant t (that is t remaining constant). The 
f{x, y, z, t) must therefore be the same as f(x -{• Sx, y + by, 
z + Sz, t), and both of them must = 0; thus 

Six + bx,y+ dy, z + Sz, i) — f(x, y, z, t) = 0. (2) 

We can place the last equation f{x -{• Sx, y + Sy, z + Sz) 
in a somewhat more convenient form. The rules of ele- 
mentary calculus will enable us to develop this function into 
a Taylor's series, that is in terms of the virtual increments, 
dx, dy and dz. Confining ourselves to the first power of the 
small increments we have 

fix + 8x,y+ Sy, z + dz, t) = f{x, y, z, t) 

df df df 

+ ix''' + dy'y + dz''- (3) 

Comparing this with (2) we see that if the displacement is 
consistent with the constraints, then 

df dl df 

must necessarily be zero; this is often written thus 

8/ = 0. (4) 

In this deduction the time t was taken as constant since (4) 
must hold true for any moment, that is must be absolutely 
independent of the time variation; and of course the con- 
straints might have been independent of the time in the 
first place, that is instead of f{x, y, z, t) = 0, we might have 
had fix, y, z) = 0, in which case the independence of the 
time would have been implied from the first. This may 
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seem a very obvious remark, but it should be well mastered 
by the reader. 

Now actual displacement (not virtual) must take place con- 
sistently with changing constraints (that is, changing if they 
contain time). In general therefore the virtual displace- 
ments dx, dy, 5z, will not be the same as the actual displace- 
ments dx, dy, dz; they will be the same only if the constraints 
are independent of the time, at least explicitly so; that is, if 

df ^ dF ^ 

— =0* — = etc 

dt ' dt ' 

3. D'Alembert's principle. This perfectly general and very 
powerful principle is seldom satisfactorily explained; yet it is 
quite easy to grasp and to apply it in practice. As a simple 
illustration of it consider a material particle moving along a 
constraining curve (fig. 1) under the action of some external 




Fig. 1. 

force E; this is usually called an impressed or applied force. 
It is quite evident that only a certain part of the force, E, 
will really act upon the particle in a manner consistent with 
the constraints. Resolving the impressed force E into two 
components, C, tangent and R, normal to the curve, we can 
say that, owing to the reaction of the constraints, the force R 
is lost or wasted, so far as the motion is concerned; and only 
the remaining force C (which is called the effective, active or 
conserved force), is able to produce the motion of the particle. 
We can readily see, therefore, that the effect of the impressed 
force E and that of the effective force C is precisely the same; 
namely, they both are capable of producing the same accelera- 

2 



6 Lagrange's Equations. 

tion and this is equal to the conserved force C divided by the 
mass of the particle. The only difference is that the force E 
has a certain effect on the constraints in form of the force R, 
taken up by the reaction of the constraint, and does not in the 
least aflfect the motion; while the force C has no effect what- 
ever on the constraints. 

Let us now imagine that the particle is acted upon by the 
same force E as before, and that two equal and opposite 
forces (fig. 2), C and C, have been applied to it as well. These 




Fig. 2. 

two latter forces are in equilibrium, and therefore do not 
aflfect the motion in anv manner: under the action of the three 
forces the motion will be precisely the same as it would have 
been under the action of the force E alone. Since C and E 
are interchangeable so far as the motion is concerned, if C 
balances C\ it is quite evident that E will likewise balance 
that force C\ Now C is merely the reversed force C and 
equals •— C (it is often called the force of inertia or the kinetic 
reaction). This is stated in D'Alembert's principle: In view 
of the constraints, the impressed forces are in equilibrium with 
the reversed effective forces. 

In extending this principle to a system (instead of a mere 
particle) the same reasoning can be applied to each separate 
particle of the system (or of a rigid body); and what we 
meant by constraints and represented as a curve to which 
the motion of the particle is confined, need not necessarily 
be a material path or track; it may also be construed to mean 
the action of one or several surrounding particles, or, for 
instance, the action of the rest of the system upon the specific 
particle under consideration. But the action of the sur- 
rounding particles upon the given particle is necessarily 
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equal to the action exerted by the latter upon the former; 
so that in a system all internal constraints, molecular actions, 
etc., automatically balance each other. The only forces which 
are unbalanced, and therefore produce the motion, are the 
impressed forces, or rather their effective components. In 
order to balance these we mentally apply the forces of inertia^ 
that is, forces equal but acting opposite to the effective forces; 
this will produce equilibrium and thus reduce the problem 
from one of dynamics to one of statics. In other words, 
given all the impressed forces, we have but to calculate what 
forces, consistent with the system, would balance them; and 
then, reversing the latter, we will have the forces actually 
producing the motion. That is, given E we calculate C to 
balance the former; then the reversed C gives us the force 
that actually produces the motion and to find which is our 
ultimate object. Hence another way of expressing D'Alem- 
berths principle: If, at any time, the moving system be stopped 
and all impressed forces, as well as all forces of inertia^ applied 
thereto, the system will remain at rest. 

Expressing this principle analytically, we can obtain the 
necessary equations of motion — our immediate aim. This is 
done in the following manner. If the impressed (or applied) 
forces are to be in equilibrium with forces of inertia (or re- 
versed effective forces, that is — mj for each particle, j being 
the acceleration of motion), we can simply express this by 
stating that the virtual work done, by such forces combined 
is for any displacement consistent with constraints. The 
expression of virtual work in rectangular coordinates, as we 
have seen under Virtual Work, is Xoix + Yoiy + Zoiz; but 
in our case any force Xq, for each particle, will consist of the 
projection of the corresponding impressed force, that is X, plus 
the projection on the same axis of the corresponding force of in- 
ertia, that is reversed effective force or mass times acceleration 
upon that axis; in other words 



Xq — X "T" 
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8 Lagrange's Equations. 

and of course the same applies to y and z. We must write = 0, 
not the work of any one, but of all the particles of which the 
system is made up; this is indicated by the sign S, so that the 
final form of D'Alembert's principle, expressing analytically 
that there is at all times equilibrium between the impressed 
forces and the forces of inertia, is 

(1) 

+ (z-mg)fe]=0. 

This is the fundamental equation of dynamics. It can be 
applied to any system of particles or bodies; to free as well as 
to constrained motion (if free, 5a:, by, iz are wholly arbitrary; 
if constrained, they must satisfy the equations of constraints). 
This equation was given by Lagrange. The differential 
equations of motion are directly derived from it. Let us 
suppose, for instance, that the system consists of n particles 
and that there are k constraining conditions 



• / 



fl{Xy y, Zy XXy J^l, %Xy ' '''''{)'= 0] 

f2{x, y, z,xu yu zi, • • • = 0; 

etc. 






connecting the particles in a certain manner. These equa- 
tions of constraints will, in general, contain 3w coordinates, 
a:, y, z, xi, yu ^u • • • of the particles, as well as the time, t. 
The virtual displacements 8x, Sy, dz of the fundamental equa- 
tion, in order to be consistent with the constraints, must 
satisfy such equations as Sf = 0; in other words 5a:, dy, 8z 
must satisfy the following k equations 

_5^ + _8y + _8, + _8a:,+ ...=0, 
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(compare with (4) under Virtual Work). It will be observed 
that t has been assumed constant, since our discussion applies 
to any instant we wish to choose. From these equations (3) 
we can find any k such variables as 5x, 5y, 5z, • • •, etc., in 
terms of the other Sn — k; introducing them into the funda- 
mental equation we shall have only 3n — k variations to deal 
with, but they are now perfectly arbitrary since k charac- 
teristics due to constraints have been eliminated and there- 
fore the constraints no longer enter into play. Now the 
remaining 3n — k variations being absolutely arbitrary, 
subject to no condition whatever, we can put each of their 
coeflScients equal to 0, which is really the only way to satisfy 
the fundamental equation (1) for any arbitrary value of these 
variations. Therefore, equating to 0, all of the coefficients of 
the remaining Sn — k variations will give 3n — k differential 
equations of the second order, which, added to the k original 
equations of constraints will enable us to find the desired 3n 
coordinates as functions of time /, and this is precisely our 
object. In integrating the 3n — A: differential equations of 
second order we shall introduce 6n — 2k constants which will 
be determined from the initial conditions of each problem. 
In eliminating the k variations, 5a:, by, dz, etc., from (1) by 
means of (3) it is much more convenient, instead of direct 
elimination, to use the following artificial method, which was 
also given by Lagrange: Each of equations (3) is multiplied 
throughout by a certain indeterminate factor such as Xi, 
X2, •••, etc.; then they are all added to (1); factoring the 
result of such addition by dx, Sy, dz, etc., we have 

+ (2-«f + X,f + X,|+...)5.+ ...]-0, 
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10 Lagrange's Equations. 

which equation involve 3n such variations as dx, 5y, 5z, 
Sxi, Syi, etc. Here the coefficients of the first k variations 
Sx, by, Sz, Sxi, iyiy etc., can be made equal to because the 
k arbitrary factors Xi, X2, X3, may be chosen so as to satisfy 
this condition; and the remaining Sti — k variations being 
perfectly arbitrary, their coefficients must be all equal to 0, 
or else there could not be equilibrium, as given by the equation 
(1) under all conditions. In other words, all such expressions 
as 

y-m^ + Xi^ + X2^+ •.., 



etc. 



(of which we have 3w), must equal 0; these 3n equations to- 
gether with the k equations of constraints will enable us to 
determine 3w coordinates and k multipliers Xi, X2, • • • , etc. 
Taking one of such equations as for instance, 

^x dfi a/2 

X-m^+Xi^ + X2^+ ..- = 0, 

we can put it in the following form 

of which the extreme importance should not be overlooked. 
Remembering that force equals mass times acceleratiouy we can 
easily see that the right-hand side represents the sum of the 
forces (projected upon the corresponding axes, x, in this case) ; 
one of these forces is X, that is the projection of the given or 
impressed force; other terms represent forces due to the pro- 
jected reaction of the constraints; so that this method is not 
only a powerful means for determining the motion proper, 
but gives the lost forces (or the corresponding reactions) as 
well. 
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We shall illustrate the application of both the fundamental 
equation and of the method of Lagrange's multipliers in the 
following easy example: 

A particle of unit mass is constrained to move in a vertical 
plane which itself is made to revolve about a vertical axis with a 
uniform angular velocity co. Find the motion. 

Let the axis of rotation (fig. 3) be the axis of Z, directed 
downward, and the other axes— as shown. In the beginning 




-^3r 



z 
Fig. 3. 

of motion, when the time / is 0, let the plane coincide with 
the plane x, z; then the plane will always be given by its 
angle cat from its initial position, and its equation in rect- 
angular coordinates will he y = x tan cat. The equation of 
the constraint in our case will hence be 

/(•P* y>^>t) = y — X tan cat = 0. 

(This equation does not contain z, because the plane is parallel 
to the Z axis, but contains t, as a movable constraint should.) 
The only impressed or applied force is that of gravity, which 
for unit mass is = ^r; therefore we have in the fundamental 
equation (1) X = 0; 7 = 0; Z = ^r. 

In order to write down the equations similar to (4), into 
which the fundamental equation (1) resolves itself, we have 
to assume a certain Lagrange's multiplier X (only one, since 
there is only one equation of constraint), so that the pro- 
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jected reactions of the constraints will be 

.c^ .dl df 

^dx' ^dy' ^dz' 

Differentiating f{x, y, z, t), that is y — x tan cot, with respect 
to the variables x, y, z, we have 

x/=-Xtanco<; x/=X; X^=0; 
ox ay oz 

so that the equations (4) will become (the mass being = 1) 

cPx ^r , ^ ^f 

Integrating the last equation twice, we have the velocity 
and the height traveled in the time t: 

where Ci and C2 are arbitrary constants. The first two 
equations can now be freed from X, which can be done by 
means of the equation of constraints, as was suggested in 
the theory. Differentiating the first equation twice and sub- 
stituting into the result the values {d^xjdf) and {(Py/df) 
from (a) we have after easy reduction 

dx 
X = 2co ^ + 2co^a; tan cot 

which, in the first equation (a) gives 

(Px dx 

-7^ + 2co tan cot- -ij+ 2co^ tan^ cot-x = 0. (b) 
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It is of advantage to introduce here a new variable x — r cos (d, 
whence 

dx U/T 

~n^ 'jj cos cot •— rco sin cot; 

and 

d^x d^r dr 

^ = ^ cos cot — 2-77 CO sin cot — ror cos cot. 

Substituting the values of dx/dt and d^x/dt^ into (6) we have 
after reductions 

^^ = ^^ > 
of which differential equation the usual solution is 



but 



X 

r = 



so that 



cos cot 



X = cos w^CCs^"* + C^e-""^). 



Having thus found both z and x in terms of the time we have 
completely determined the motion. 

Another easy example illustrating D'Alembert's principle is 
as follows: Two particles, P and P' {of mass m and m') are 
fastened to a light rod and made to swing like a pendulum. 
Find the motion. 

Let (fig. 4) AP = a and AP' = a'. According to D'Alem- 
bert's principle we must imagine the system stopped; then 
there will be equilibrium under the action of the applied forces 
and the forces of inertia. In our case the apphed forces are 
the weights mg and m^g of the particles; their moments about 
A are — mag sin d and — m^a^g sin d (minv^, because they 
decrease the angle; moments, instead of forces because one 
point of the system is fixed, which is always a hint). The 
forces of inertia are m{d^s/dt^) and m\d^s'/df); and their 



14 



Lagrange's Equations. 



% 



moments are ma{(Ps/dP) and m^a^{cPs^/df) {s and s^ being the 
arcs described). Therefore 

cPs d?s' 

^^ 77/2 "I" ^'^' j^ + (^^ + 'ni'a')g sin S = 0. 

But cfe = ads and ds' = a'dOy so that 

(ma^ + mV^) jj!^-}- (jfia + m'aO^f-sin S = 0. 

Comparing this with the problem of a compound pendulum 
(Bowser, Anal. Mech., p. 458) we see at once that the system 



^qA 




Fig. 4. 



will oscillate precisely as a simple pendulum of the length 



/«/ 



ma + m'a 

D'Alembert's principle and Lagrange's fundamental equa- 
tion (1) are perfectly general and can be applied to any 
problem involving either equilibrium or motion. Several 
useful principles can he derived from these, which themselves 
help us to solve a great many problems, without applying 
the general equations. We shall mention only two such 
principles; that of kinetic energy, and of areas. 
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4. Integral of kinetic energy. The principle of vis viva or 
kinetic energy is akeady familiar from the elementary course 
(Bowser, Anal. Mech., p. 489), but it may easily be deduced 
from the fundamental equation of dynamics as follows: In 
estimating the value of the virtual work of a system we always 
referred it to some particular moment of time t, that is, we 
imagined the constraints momentarily fixed for that instant. 
The displacements in that case were denoted by 5x, 5y, 5z, 
etc. But, if the constraints contain the time, then the real, 
actual, displacements dx, dy, dz, cannot, generally, be equal 
to the virtual displacements, Sx, dy, 5z, which were subject 
only to one condition, to be consistent with the constraints, 
temporarily fixed for that moment. For instance, if the 
constraints are represented by a surface A (fig. 5), we can 




Fig. 5. 

imagine a number of displacements a, b, c, d, e, • • • , on that 
surface for a given instant, during which the surface is sta- 
tionary; while in reality, owing to the motion of the surface 
itself, the actual displacement may be such as m/, that is 
entirely different from any virtual displacement we imagined 
for the case of a fixed surface. Omitting the case of movable 
constraints, however, that is supposing that the equations of 
constraints do not contain the time t, we can always select 
the virtual variations 5a:, 5y, 5z, etc., so that they will be 
equal to the actual displacements, dx, dy, dz, • ". In this 



16 Lagrange's Equations. 

case (and in this case only) we have the fundamental equa- 
tion of dynamics in the following form 

2[(jr-«^)&+(y-».g)* 

which can be rearranged thus 

i:m(^^dx + ^,dy + ^dz^ = li{Xdx+ Ydy+Zdz). (1) 

The second member of this equation represents the work, 
not virtual but actual, done by the system (compare with 
the expression of virtual work given above). It is often possi- 
ble to calculate a special function TJ of the coordinates (that 
is only of Xy y, z, and not of the time or of time derivatives, 
etc.), of which the partial derivatives, with respect to axes, 
will give forces acting along these axes. Such a function is 
called the potential function and exists in a great many cases. 
If such a function exists, the right-hand side of (1) can be 
further simplified. 

Indeed, from the very definition of the potential function 
we have 

dJL-x- ^-y. ^JL-z. 

dx~ ' dy~ ' dz~ ' 
which in the right-hand side of (1) gives 



(ll'^ + ^'^2/ + ^^) 



and this a total diflFerential, dU, of the potential function U 
(Bowser, CalciJus, p. 122). This will reduce the equation 
(l)to 



^ fd'x^ , d^y 
^"^[l^^+df 



2 



dy + -jT^dz j = dU. 



General Principles op Dynamics. 17 

In order to further simplify the first member we will use the 
well-known expression of velocity (Bowser, Anal. Mech., 
p. 244), 

diflFerentiating this with respect to time t we have 

so that finally we have \i^tw? = dJJ\ integrating this we 
have 

Sm^- Smy= f7- V^ (2) 

where the initial velocity i?o and the initial value J7o are given 
by the initial conditions of the problem. This result is known 
as integral of kinetic energy; it shows that the kinetic energy 
acquired (or lost) is equal to the diflFerence of the initial and 
the final values of the potential (or force) function, no matter 
what path was followed in changing the position of the system. 
In other words, if the potential function exists, the increase 
(positive or negative) of the kinetic energy, between any two 
positions of the system, will be independent of the manner in 
which the change was performed. 
Substituting 

in equation (1) we have the expression 



f=X<" 



Sm ~ Sm ^ = I {Xdx + Ydy + Zdz), (3) 



that is, the increase of kinetic energy equals the work done 



/ 
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by all the forces acting upon the system in the interval of 
time between to and /; this applies also to the more general 
case, when there is no potential function. 

By all the forces we mean all external, internal, applied or 
reaction forces; this latter expression applies whether the 
constraints are moving or not, to any system consisting of 
solid or liquid matter; it should be committed to memory by 
the reader merely as the kinetic energy equals the work done. 

In our future work we shall often be called upon to calcu-. 
late the kinetic energy of motion which is given, not in x, y, z, 
but in some other, less familiar, system of coordinates. In 
order to obtain this, all we have to do is to find x, y, z in terms 
of the new coordinates, then differentiate such expressions 
with regard to t and finally to substitute the results in 



v" 



=( 



dt ) ^\dt 



)H^)'- 



For instance, let the motion be given as taking place on a 




Fig. 6. 



x 



sphere of variable radius, r, and characterized, say, by lati- 
tude, d, and longitude, <p. We can easily see that (fig. 6) 

X = r cos d cos <f>; y =" r cos 6 sin <p; z = r sin 6, 
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and substituting the time-derivatives from these formulae 
into 1? we have, after easy reduction, 

7? ^m dr^ -\- rMd^ + r^ cos^ dd(p^ 
^2 "" 2 de • 

In the case of plane polar coordinates, r, ^, this reduces 
{e being 0) to 

m dr^ + r^d(p^ 

2 df ' 

or, in case of ordinary rotation (about an axis, when r = const.) 
it further reduces to 



2^ \dt J ~ ^ 2' 



m 
2 



where I is the moment of inertia about the axis. This is a 
well-known expression of kinetic energy for rotation about a 
fixed axis. 

To return to the integral of kinetic energy: it is called 
integral because it is a so-called first integral of motion; if put 
in the form 

Xm^= U+h 

(where Zm{vo^/2) and Uq are absorbed in one constant, h), we 
see at once that the second derivatives, dJ^xjdt^, etc., which 
appear in the fundamental equation, have vanished, and in 
the integral of kinetic energy we have only the variables 
^> y> z> " '> (dx/dt), (dy/dt), (dz/dt), • • •, and one constant, h; 
such an expression itself is a partial solution of (1). 
It is customary to denote kinetic energy by T, so that 

dT = i:{Xdx + Ydy + Zdz) 
and, if there is a potential function, then 
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This form of the integral of kinetic energy will be constantly 
employed. It has no meaning when the constraints are vari- 
able, because then there is no potential function U; but it 
exists in a great number of cases, — in •all sorts of molecular, 
attractive (gravitational) forces, etc. 

Let us take a few es^amples: For a particle under the action 
of gravity, the axis z being directed downward, the potential 
function will be Z7 = mgz, because 

dx dy ' 
as they should; while 

dU 

-dz = ^^' 

that is the weight of the particle is the only moving force Z. 

In all problems involving suspended or pivoting systems z 
is generally the vertical distance of the center of gravity 
from the fixed point (or axis). 

If the axis z is directed upward, U would be = — mgz. 

For a particle which is repelled inversely as the cube of 
the distance we can at once write down 

because the derivative of this expression with respect to the 
line of action, x, will give mk/x^. 

For a particle moving in a vertical plane, under the action 
of gravity, along a constraining curve, given in polar coordi- 
nates, r, d, we have U = mgz = mgr sin S (if ff is reckoned 
from the horizontal axis x). 

In problems involving turning moments (man turning a 
crank with a uniform moment M) the potential function will be 

U= ilf tan-i- 

x 

(if y is the axis of rotation). 
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Indeed we have the following partial derivatives 

dU^ Mz 

dz ^ «^ + J^ Z 
and 

dU Mx 



. ,. I 



dx «^ + J^ * ir- 

Since a^ + ^ ^ f^ (the radius of the crank) we see that 
these partial derivatives are equivalent to — (M/r) sin <p and 
(M/r) cos <p, where <p is the angle of the radius with the axis x; 
in other words they represent forces along the axes z and x. 
This remark will be of practical interest in engineering work. 

For a particle subject to the action of a force radiating from 
a center and proportional to the distance, we have 

u=- — 

2 ' 
where a is a constant; because then 

dU 

as required. 

Example. The motion of a system is given by the following 
equation 

-Tg = (cos X — a) sin x, where a > 0. 

f^ind the potential function U. 
Multiplying by 2{dx/dt) we have 

^dx cPx d / dx\^ ^. ■ . dx 

Now, since the motion is evidently rectilinear, we have 

__ mv^ __ m / dx\^ 

^~X""2v5ry' 
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or, from the above, 

T ^ f (cos X ^ a) sin xdx + c = J sin^ x + a cos x + c. 

But, since T = U + h, we have the required solution 

U = i sin^ X + a cos ar + 6, 

where 6 is a constant, absorbing c and A. 

We shall say a few more words regarding the force func- 
tion U and the practical methods of finding it. 

Since by the very definition of the force function 

dU dU dU 

therefore, in the general expression of elementary work we 
have 

dW= Xdx+ Ydy + Zdz = ^dx+ Y-dy+^dz^dU; 

that is elementary work = complete differential of potential 
function. 

From this we should not conclude, however, that any 
function, of which the partial derivatives with respect to 
axes give projected forces on these axes, is the potential 
function, for this is not always the case; the following is the 
criterion: we have 

^=x- — =y- — =z. 

dx ' dy * dz ' 

Differentiating the first equation as to y and the second as 
to X we have • 

d^U ^ dX 

dxdy dy * 
also 

d'^U dY 

dxdy dx ' 
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so that 



Similarly 



also 



dX 


dY 


dy 


dx' 


dX 


dZ 


dz ~ 


dx ' 


dY 


dZ 


dz 


dy ' 



If these equations are identically satisfied, then (and then 
only) the potential function U exists and can be found by 
integrating the expression Xdx + Ydy + Zdz, giving, as we 
have already seen, 

"2 2"= ^"^«' 

which of course is nothing more than 

-$-""f=jiXdx+Ydy+Zdz). 

Example. If, for instance, the acting force is merely gravity, 
the axis z being directed downward, then X = 0\ 7=0; 
Z = mg. Here our criterion is identically satisfied and 
dU = mgdz, hence U = mgz + c; the constant, c, as a rule 
is taken = 0; although later on it will be shown that at times 
it may be convenient to assign it a certain value, in accordance 
with the problem (see small oscillations). 

Example, A particle, x, y, z, is acted upon by a force, 
emanating from a center, the intensity being a function /(r) 
of the distance from the center. Find the potential function. 

Let us take the center as origin of coordinates. Then, the 
distance being = r, the cosines of r with the axes will be 
— xfr, — y/r, — z/r; and if the force is = /(r), then its pro- 
jections on the axes will be of course 

X=-/(r)p F=-/(r)f; Z=-/(r)^; 



24 Lagrange's Equations. 

and therefore 

Xdx+ Ydy+ Zdz^ -■^(xdx + ydy+zdz). 

T 

However, a^ + y^ + z^ = r^, so that xdx + ydy + zdz = rdr) 
and 

Xdx+ Ydy + Zdz ^ ^^-^ .rdr-- - f{r)dr = dU; 

so that 

D^= - f f{r)dr 
Jo 

is the potential function, if only 

dX_dY 

dy'^ dx' ^^^'' 

which we can easily prove is the case. Indeed, the condition 

dX_dY 

dy dz 
means 

or, since this is partial differentiation, where x is considered 
constant when deriving with respect to y; also y is considered 
constant when deriving as to x, we have simply 

dy r ^ dx r ' 

or, applying the rule of differentiating functions of functions 
we have 

dr\ r J dy ^ dr\ r Jdx' 
But from r^ = x^ + y^ + z^ we have 

dr y dr X 

■r- = — and ^~ = - > 
dy r ox r 



s 
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hence our conditions become 



r dr\ r J r dr\ r J ' 



an identity. The same would apply to the other two condi- 
tions. Hence U is the potential function sought. 

Example. The force is such that its components upon the 
axes, X, y, z, are functions of the coordinates: X = fi{x); 
Y = f2{y); Z = /sCz). Find the potential function. Here 

Xdx+ Ydy +Zdz=d[ ffi{x)dx + ff2{y)dy + fMz)dz] = dU; 

Hence U = [ jfi(x)dx+ ff2{y)dy + J Sz{z)dz]y provided that 

the criterion is satisfied. 

We should not feel in any way compelled to use rectangular 
coordinates in compiling dW or dU. Thus, instead of 

Xdx + Ydy + Zdz = dU 

we might have Pdp + Qdq + Bdr = dU, if the position of 
the particle is defined by some special coordinates p, q, r. 
For instance in polar coordinates, r, 0, we have 

Mr + Qrde = dU. 

In semi-polar or cylindrical coordinates, r, 9, 2, we have 
Bdr + QrdB + Zdz = dU. It is easy to see why in both 
cases we have Qrdd and not simply Qdd; because the latter 
(force times angular displacement) is not work. Now 6r is a 
moment and moment times angular displacement is work. 
So that the method of finding the potential function is always 
the same: find expression of elementary work and reduce it, 
if possible, to a complete differential of some function. If 
the criterion is satisfied, then the indefinite integral of such a 
function is J7. 

One more example will illustrate the mechanism of such 
reduction; it forms part of a problem given in Chapter V 
(see fig. 52). 



\ 
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An inspection of the working sketch of this system supplies 
tiie following expression for elementary work, corresponding 
to elementary displacement dtp 

dW= - \AMd(AM) - fiBMd(BM) + Md{c sin ip)] 

now this does not look like a complete differential of any 
function XJ\ but, arranging it in the form 

dW = Mc cos (pd<p - idiUAMy + ix{BMf), 

also taking into account the general relations 

{AMY = a^ + (^ — 2ac cos (p, 

(BMy = a^ + (^ + 2ac cos <p, 
we have 

dW = c[a(jjL — X) sin ^ + jJf cos <p]d(p = dU; 
hence U is the indefinite integral of dU 

U = c[M sin <p — a(ji — X) cos (p] 

(the constant does not matter for the present; later we shall 
be very specific in selecting it). It will be noted that in this 
case our only coordinate is angular, cp, and that the derivative 
of f7 as to ^ gives us the force corresponding to that coordinate, 
that is a certain moment, not an ordinary force; because a 
Tooment multiplied by an angular displacement gives work. 
A force does not. The force understood in this broader sense 
is called generalized force, and in our future work we shall 
have a great deal to do with such forces. 
Remembering that 

it will be seen that the integral of kinetic energy affords a 
very powerful principle for solution of problems of dynamics. 
Example. A small material particle of mass m, slides, 
without friction, off a cylindrical log; to determine the point 
at which it leaves the log. 
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Instead of applying the fundamental equation, we can 
simply use the integral of kinetic energy. We know that the 
pressure exerted by a particle on a frictionless curve is equal 
to the centrifugal force, m • {^Ir) minics the radial component 
of the applied force (Bowser, Anal. Mech., p. 348). It is 
clear that the particle will leave the log at a point where 




Fig. 7. 

this pressure will vanish. The equation of the pressure is 
P = miv^lr) — mg cos <p (fig. 7). Now, when this is = 0, we 
have 

t)2 ^ y 

r " ^ r 

or D^ = gy. Since the acting force is that of gravity, we have 
the potential function f7 = — 'mgy, so that 



v" 



^o' 



^2"^y== 



— mgy + mgyQ) 



substituting ir^ = gy we have 



y = 



2gyo + ^0^ 

3(7 



If, at the beginning of motion, the velocity was 0, we would 
have had y = fyo, that is the particle would have left the log 
after having descended one third its initial elevation above the 
center of the log. 



^ 



1 
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It will thus be seen how a judicious introduction of the 
integral of kinetic energy simplifies the work; it is itself a 
partial solution of the problem. Another example will be 
given in the next article. 

It is important to remember that the integral of kinetic 
energy was evolved from the fundamental equation, that is 
from the most general case that can be imagined, with the 
only limitation, that the constraints were independent of the 
time. Therefore the integral of kinetic energy can be applied 
to constrained as well as to free motion. This is evident 
from the fact that when we consider the work of all forces 
to be equal to the increase of kinetic energy, we imply the 
idea of the reactions being normal to the constraints and 
therefore producing or absorbing no work (we do not consider 
the effect of friction). So that in the case of constrained 
motion we also have 

dS-„- = S(Xcfo; + Ydy + Zdz), 
if there is no potential function and 

if there is such a function, in which latter case we might 
also say that 



mi? 



= U+h; 



which is what we represented by T. 

Synopsis of notations relative to work and energy. Virtual 
work SW ^ force times displacement, Fds, that is force times 
any displacement projected upon the direction of the force, 
= Pdp cos (P, dp). It also equals the sum of virtual works 
done by the components of the force upon their respective 
axes; that is, = Xdx + Ydy + Zdz. In polar coordinates. 
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virtual work — moment X angular displacement^ Md<p. The 
difference between 8 and d is the arbitrary character of the 
former and the actual nature of the latter. D'Alembert's 
principle and the corresponding fundamental equation apply 
to any motion under any circumstances whatsoever, but can 
be simpMed for special cases; thus for permanent constraints 
(in both shape and position), the formula (1) can be split 
into (2), of which the right side is elementary work, and the 
left side can be reduced to the form dSm(v^l2); this was 
briefly expressed thus: kinetic energy = work done, by which 
we mean that the increase of kinetic energy can only take 
place owing to external forces, and equals to the work done 
by them. We then agreed to denote the kinetic energy by T 
so that dT = ^{Xdx + Ydy + Zdz) is the same as (2) only 
in differential form. A still further reduction was made for 
the special case when l^{Xdx + Ydy + Zdz) could be shown 
to be an exact differential of some potential function U (of the 
coordinates only), that is, a function whose partial derivatives 
with regard to any coordinate equaled the forces corresponding 
to that axis. For this special case we had dT = dU or 

idi:mv^=dU; or T =^ U + h; or S-|-=f7+A; 

each of which expressions is called the integral of kinetic 
energy. Such expressions exist only when the constraints are 
permanent and when, in addition to this, there is a potential 
function U; if there is no such function, there is no integral 
of kinetic energy but the general expression 

dT = S(Zda: + Ydy + Zdz) 

still holds true. In our future work we shall use various 
expressions of kinetic energy or of its integral in almost every 
problem; they should be well understood and memorized. 

5. Integral of areas. This is another "short cut" principle, 
supplying a convenient method for beginning the investigation 
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of a problem; we shall give it in a form, slightly different from 
that adopted in Bowser's Anal. Mech., p. 487. 

Let us write down the fundamental equation of dynamics 
for the case when a system can receive a virtual displacement 
only about the axis z; this will be an angular displacement, 
say da, and therefore it will be well to adopt polar coordinates, 
r, a, where r will be kept constant for each particle of the sys- 
tem (that is, no motion, except rotation about z). The usual 
formulae of substitution will be a: = r cos a; y = r sin a; 
z = 0, we therefore have 

dx = ^ r sin ada and dy = r cos a8a, 

or else, Sx = ^ y8a; by = xha, and these values in the funda- 
mental equation will give 

which is usually written in the form. 



dt 



2-(4-^f)=^(^^-^^)- (i> 



The right-hand member is well known from statics; it repre- 
sents the sum of the moments of all forces about the axis z; the 
first member involves quantities [m{dy/dt)]x and [m{dx/df)]y; 
the expressions such as m(dy/dt), etc., are known as the 
momentum of a moving particle projected upon the axis y; 
they are known from elementary courses (Bowser, Anal. 
Mech., p. 7) simply as mv, or mass times velocity. Each of 
these expressions, when multiplied say by x, is known as 
moment of momentum, or angular momentum, about the axis z 
(the axis which it does not contain) ; the same notations apply 
to the expression [m(dx/dt)]y; so that (1) can be expressed 
as follows: The time-derivative of the sum of the moments 
of momenta of all particles of a system about any axis, is 
equal to the sum of the moments of external forces about 
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that axis. The choice of the axes of coordinates is purely 
arbitrary and therefore we can select any line whatever as 
the axis about which the estimate of the moments of external 
forces is made; then the sum of angular momenta of the 
moving system about that axis, if differentiated with respect 
to time, will give the exact value of these external moments. 
If, however, there are no external moments, the time-deriva- 
tive is and S is constant. Thus, 

(dy dx\ 

Substituting polar coordinates x = r cos cp and y = r sin <pj 
we shall have 

^mr^ -37 = a 
at 

or, integrating, 

S I mr^dip = a(f, — <o)j 

whereby the sum of areas described by the projected radii- 
vectors of each particle, each multiplied by the corresponding 
mass, is proportional to the time. This is the so-called 
integral of areas and is often used in either of the forms 

^ ( dy dx\ ^ c,d<p .^. 

Sm(^x^-y-^j = a; or Smr^^=a. (2) 

Similar expressions obtain for the other axes. This integral 
can be used when there are no external forces acting on a 
body; or when the resultant of all forces passes through or is 
parallel to the axis about which it is desired to apply the 
integral of areas; in other words when there is no moment 
about that axis. 

The following interesting problem will illustrate the appli- 
cation of this principle. 

Example. A particle of unit-mass is constrained to move on 
the surface of a paraboloid of revolution oi? + y^ = 2az, whose 
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axis z is vertical. The initial velocity of the particle is horl- 
zontaL It is required to calculate the reaction of the surface in 
terms of z. Also to discuss the case when the initial velocity is 
= y^2gzo, where Zq is the initial value of the coordinate of the 
particle, and finally to show that the particle will move in a 
circle in a horizontal plane and that the period of revolution will 
he the same for all such parallel circles. 

The impressed forces reduce themselves merely to the 
weight of the particle: Z=0; y=0;Z= — g. The poten- 
tial function is f/ = — gz. The integral of kinetic energy 
(m being = 1) is 

2 " y = f^ " f^o = - 5^(2 - zo); 

but from the general expression of velocity in rectangular 
coordinates 



therefore 



"-(ly+dr+d)'^ 



(iy+(i)"+(i) 



(^\ = 



2g{zo— z) +j)Q^. 



The surface being one of revolution about z, it is evident 
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that whatever the reaction is it must be normal to the surface, 
that is that it must pass through the axis of z, and therefore 
the reaction will give no moment about z. Neither will the 
impressed force (gravity), for it is parallel to the axis z (fig. 8). 
This is, therefore, a typical case where the integral of areas 
should be applied, viz., 

dy dx 

where a is to be found from initial conditions. Let, for in- 
stance, the particle start from a position in the a; — z plane, 
so that yo = 0; we also have, the initial velocity being hori- 
zontal, (dyo/dt) = vq; therefore the integral of areas for these 
initial conditions reduces to aro^o = o = i^o^o, where r© will 
represent the initial distance of the particle from the z axis. 
Having thus found the constant a, we can write the general 
expression of the integral of areas as 

dy dx 

Substituting polar coordinates (x ^ r cos <p and y = r sin ^), 
we can transform the equation of the surface into r^ = 2az; 
the integral of kinetic energy becomes 

and the integral of areas will be 

2d(P 

'-di = ''''' 

Differentiating the equation of the surface, r^ = 2az, we 
have 

dr __ dz 

''di''''dt' 

substituting the values of {d(p/dt) and (dr/dt), found from the 



i 
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last two equations, into the integral of kinetic energy we have 

7 \dt) + ""r^+ U/ ^ ^^^'" " '^ + '"'• 

In order to eliminate r^ and r^, we can substitute their value 
from the equation of the surface, that is r^ = 2azo, and 
r^ = 2az, after which, solving for dz/dt, we have 



\dt) " 



2z[2g{zo - z) + Vq^] - 2i?o^2o 
a+2z 

4/7 



a+2z 
This is satisfied by putting 



(z-Zo)(|-z). 



fo^ 



indeed, in that case, 






2 = 77-= const. and tt = 0. 
25f dt 

Now z = const, means that the motion will be in a plane 
parallel to x — y. Substituting z just found in the equation 
of the surface we have 

whence the angular velocity 

^0 [g 

r ^ a 
and the period 

m 27r ^ fa 

CO yg 

in other words the period is independent of z. In this typical 
example it was possible to use the integral of kinetic energy 
as well as that of areas. In the example given under D^Alem- 
berfs principle (a particle in a plane rotating about a vertical 
axis) no such application could have been made; the moving 
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constraints (containing t) made it impossible to use the 
integral of kinetic energy (see how this was derived for the 
case of constraints involving no time) ; likewise, although the 
impressed force mg does not give any moment about z, the 
reaction of the plane does (being perpendicular to the plane), 
so that the integral of areas cannot be used either. 

In problems of that sort it is necessary to apply the general 
methods, the fundamental equation, or at least D'Alembert's 
principle in its original form, as given above. 

6. Coordinates of a rigid body. When we speak of a system, 
we mean two or more material particles, in some manner 
connected together. This connection will prevent the abso- 
lute freedom of the motion of each particle and will therefore 
constitute what we have characterized as constraints. 

We shall principally have to do with internal, material, 
kinematical constraints, such as for instance a string by means 
of which the distance between two stones cannot exceed a 
certain value; or a car track; or, say, a rigid body, all the 
particles of which are subject to the condition that the original 
form of the body is permanent; etc. But some dynamical 
constraints can be imagined, acting in precisely the same 
manner as would kinematical constraints, but involving no 
material means of any kind; the solar system is an example 
of this sort, where the planets under the action of the so-called 
forces (the true nature of which we do not understand), 
describe paths just as definite as if they were running on 
tracks. 

At any rate the idea of constraints necessarily implies our 
knowing all, or at least something of the relative position of the 
particles constituting the system. If we have n particles, of 
which the coordinates are, say, x, y, z, Xi, yi, zi, • • • , etc., 
then any geometric relation of such form as 

f(x, y, z, xi, 2/1, zi, ••0 = 0, 
will mean some sort of a condition, in view of which the 
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particles are not altogether free to move; it is called^ as we 
have seen before, an equation of constraint. We may have 
several such constraining equations /i = 0; /2 = 0, • • •, etc., 
except that their number cannot be equal to 3n, the total 
number of coordinates. If such were the case, we could 
definitely establish, from these constraining equations, all 
the 3n coordinates of all the particles; in other words, the 
system would then be permanently fixed, and there could be no 
question of either equilibrium or motion. If the number of 
constraining equations is 3n — 1, that is one less than the 
number of coordinates, then each particle can move only 
along one definite path, and moreover, such motion will cause 
perfectly definite displacements of all other particles of the 
system. 

This would be the so-called system with one degree of free- 
dom; its position can be fully determined by some one charac- 
teristic, for instance by one given point on any of the paths 
of the particles. Such characteristics, which are not neces- 
sarily coordinates in the usual sense, but may be distances 
from fixed points or angles (as will be shown later) are called 
parameters. If the constraining equations are 3n — 2, that is 
two less than the number of coordinates, then each particle, 
instead of moving along a definite curve, will be able to move 
upon a certain surface, and two characteristics will be required 
(not necessarily coordinates in the old-fashioned sense, but 
distances or angles from certain fixed positions) to fully locate 
one of the particles on its surface; and that will fully determine 
the position of the system. In other words there will be 
two degrees of freedom, and two parameters necessary to char- 
acterize the system. In general, if the number of particles 
is n, and the number of limiting conditions (constraints) is k, 
then the number of parameters necessary to fix the position 
of the system will he Sn — k; they may be either coordinates 
or other typical characteristics in the sense indicated above. 
They are to be determined from the equations of motion; so 
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that, in general, the problem will involve 6n — k variables 
to be determined: 3n coordinates of the n particles, and 3n — k 
parameters. 

In order to do this in a definite way we must have the same 
number of equations; now we have available 3n — k equations 
of motion, also the like number of equations connecting the 
parameters with the coordinates of the particles, and the k 
equations of constraints, or 6n — k equations, in view of 
which the problem will be definite. By way of illustration, 
let us take a crude example of two particles given by their 
coordinates Xi, y\, Zi, and 0:2, ^2, Z2; the fact that they are two 
in number does not constitute a system. If we impose some 
condition, for instance, fixing the distance between the par- 
ticles, which is expressed very easily by rules of elementary 
space geometry, then this connection will change two separate 
particles into a system of two particles; here we have 6 coordi- 
nates and one equation of constraint 

(xi — X2y + (2/1 — 2/2)^ + (zi — 22)^ — P 

= fi{xi, 2/1, zi, X2, 2/2, Z2) = (1) 

and the system will thus have five degrees of freedom. That 
is to say one of the particles can have any motion whatsoever 
(three coordinates) and the other can be anywhere on a 
sphere, radius /, described about the first (two more coordi- 
nates, angular or any other). Taking the other extreme 
case suppose now that we have five constraining conditions, 
such as for instance Xi = a; yi = b; Zi = c; 2/2 = d and 
finally /i = 0, as (1). The first three equations permanently 
fix the first particle; the fourth equation, 2/2 = d, means that 
the second particle can move only in a plane parallel to the 
plane x — z; and /i = 0, maintains the distance between the 
two particles equal to /. Therefore our system will be capable 
of moving in a manner shown in fig. 9; the apex of the cone 
will be the point Xi, 2/1, Zi; the length of its generator will be I 
and the end of the latter, that is the movable particle X2, y^y z^, 
4 
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will be capable of assuming any position on the circle in the 
plane B, parallel to the plane x ~ z. Since we have five im- 
posed conditions and six coordinates of the particles, the 



Fig. 9. 

system will have one degree of freedom, that is to say, only 
one parameter will be required to definitely locate the whole 
system. This parameter may be selected in a great variety 
of ways; for instance we might take the distance MO; or the 
angle $ of the line MO with the y axis; or the coordinate MN; 
or PN, or PM, etc. Our problem will then be to express 
with the aid of the equations of motion, the parameter as a 
function of time, since this will determine the position of the 
system for any given time. Of course great care will be 
exercised in selecting the parameter, so that the labor involved 
in solving the equations of motion will be a minimum. 

We shall now take another supposition, namely, that there 
are only four imposed conditions, that is that the system has 
only two degrees of freedom. Let, for instance, xi = a; 
pi = b; ifi = d and /i = as before. The difference between 
this and the preceding case is that the apex of the cone itself 
will here he indeterminate, except that it will move on a line, 
parallel to the z axis; therefore it will be necessary to have two 
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separate parameters to characterize such a system. One of 
them may for instance be the same as before, say MP, and 
the other simply the missing coordinate z; or again MP 
together with PN can be chosen; or two angles (p and d, etc. 
A rigid body may be said to consist of a very great number 
of separate particles, and we might think therefore that 
three times as many coordinates might have to be determined 
in our problem in order to characterize the whole body in 
its motion. Yet, in view of the constraints (mutual action, 
maintaining the original shape of the body) only six param- 
eters are necessary. What they are, is implied by the six 
equations of elementary statics: in order to have equilibrium 
(that is the number of parameters equal to the number of 
coordinates) we must have six conditions, three of which 
refer to the center of gravity of the body and the other three, 
to the angular equilibrium about the axes x, y and z. Any 
coordinate of any point of the body can be readily expressed 
through these six parameters. By way of reviewing the 
above it is useful to mention the following: 

1. A free particle has three degrees of freedom, and three 
parameters (coordinates) are necessary to identify it. 

2. A particle constrained to remain upon a given surface 
has two degrees of freedom. 

3. A particle moving upon a curve has one degree of freedom. 

4. A rigid body, if free, has six degrees of freedom. 

5. A rigid body having one point fixed has three degrees of 
freedom (rotations about the three axes). 

6. A rigid body rotating about a fixed axis has one degree 
of freedom, and only one (angular) parameter is required to 
completely determine the position of the body. 

We shall return to this important subject presently and 
show how Lagrange has based his wonderful method upon 
this fruitful conception: the possibility to identify a system 
by means of some characteristic parameters, the number of 
which is less than the number of separate coordinates of the 
particles. 
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7. Motion of rigid body. Many extended treatises have 
been written upon this subject; in the following few pages we 
can only briefly mention a few main principles involved in 
the study of the kinetics of rigid bodies. The reader should 
carefully review his stock of knowledge already acquired 
from elementary treatises (for instance Bowser, Anal. Mech., 
Chapters VI, VII and VIII) in addition to which the following 
remarks can be made. 

(a) Moments of inektia. We know that there exists, 
about every point of a rigid body, a momental ellipsoid, which 
is the locus of all points, such that for each point the expression 
1/r^ (the reciprocal of the square of the corresponding radius 
vector) is equal to the moment of inertia about the axis coin- 
ciding with that radius vector. In every such ellipsoid there 
are three rectangular axes, called principal axes of inertia, 
about which the products of inertia vanish, so that the equa- 
tion of the ellipsoid contains only the principal moments of 
inertia A, B and C 

Axi^ + Byi^^- Czi'= 1. 

Knowing A, B and C we can readily express the moments 
of inertia of the body about any line passing through the center 
of the ellipsoid (whether this is the center of gravity of the 
body or not) and given by its angles a, /3, y, with the principal 
axes. Such a moment of inertia will be 

H = A cos^ a + B cos^ /3 + C cos^ 7. 

In this case the corresponding radius of gyration can be 
shown to be 

P = a^ cos^ a + 6^ cos^ /3 + c^ cos^ 7, 

where a, h and c are the radii of gyration corresponding to 
A, B and (7. 

From this equation the following deductions can be made: 
(1) The axis of the greatest moment of inertia coincides with 
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the smallest axis of the ellipsoid and inversely. (2) If two 
of the principal moments of inertia are equal, say A = B, 
the ellipsoid is one of revolution about the third axis C; in 
this case all axes forming the same angle with C will have the 
same values of the moments of inertia; and all axes perpendicu- 
lar to C will be principal axes; with reference to any of these 
the moment of inertia will = A = B. (3) In order that a 
line may be a principal axis of inertia relatively to any of its 
points, it is only necessary to have the products of inertia 
vanish, that is ^myz = and ^mxz = (if the proposed 
line be taken as the z axis). In this connection it might be 
well to caution the reader against thinking that a principal 
axis must necessarily pass through the center of gravity of 
the body; principal axes may have as an origin any point 
within the body; however a principal axis passing through 
the center of gravity possesses a curious property, viz.: it 
will remain the principal axis of any set of axes built on any 
of its points. 

A second important theorem, is the well-known one that 
the moment of inertia of a body about any axis is equal to 
the moment of inertia about any axis through the center of 
gravity, parallel to the given axis, plus the mass of the body 
multiplied by the square of the distance between the two 
axes. Hence the fact that the moments of inertia about the 
axes passing through the center of gravity are the least of all 
parallel axes; also, that the moments of inertia about all parallel 
axes, equidistant from the center of gravity, are equal. 

Combining the above we may easily find the moment of 
inertia about any axis whatever, no matter where located and 
how inclined, if the principal moments of inertia about the 
center of gravity are known. 

(6) Body rotating about a fixed axis. The moment of 
inertia can now be introduced into the formula already known as 
the integral of areas 
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Here the first part (as has already been mentioned) is known 
as the sum of the angular momenta (or moments of momen- 
tum); the sum 'SmT^{d<p/dt) consists of two factors, of which 
the first 'Emr^ is the moment of inertia (about the axis z), and 
the second {d<p[dt) is the rate of change of the angular dis- 
placement, known as angular velocity. For a rigid body 
instead of the sum of angular momenta we might simply use 
the expression angular momentum; so that finally the angular 
momentiun, H, of a body is equal to the product, Zco, of the 
moment of inertia into the angular velocity (all taken about 
the same axis). This theorem is already known (Bowser, 
Anal. Mech., p. 454) and will be constantly used in our future 
work. 

Another interesting theorem gives the kinetic energy of a 
rotating body as half the product of its moment of inertia 
into the square of its angular velocity: the velocity of a rotat- 
ing particle is equal to the radius multiplied by its angular 
velocity v = rco; therefore 

^ mtF ^ mcoV co^ ^ „ ' co^ 

which also may be written = Mlc^{oi^l2) (where k is the 
radius of gyration). 

Still another theorem of equally great importance follows 
immediately from the equation (see under Integral of areas) 
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We have just seen that l!,m[x{dyldt) — y(dxfdt)] is the angular 
momentum, H, of the body and = Zco; therefore its time- 
derivative dH/dt will be 



in other words 
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so that the moment of forces about the axis of rotation equals 
the moment of inertia times the angular acceleration. All 
these theorems have been given in the elementary course, 
but it is of great advantage to see how easily and quickly 
they can be derived either from the fundamental equation of 
dynamics or at least from one of the great principles, following 
therefrom. (Remark. — ^The last three theorems have been 
derived only for the case of a body moving about a fixed axis. 
That of a body moving about a fixed point will be considered 
presently.) As an example let us take a ring (of a rectangular 
cross-section) whose outer and inner radii are a and b; its 
moment of inertia about its central axis will be immediately 
found by an easy integration (Bowser, Anal. Mech., p. 448) 
and its moment of inertia about any diameter of the ring 
will be one half of that (Ibid. Remark on Polar moment of 
inertia, p. 436). It is often customary to denote, in a body of 
revolution, the equal moments of inertia by A and B, while 
the moment of inertia about the axis of revolution is generally 
referred to as C. Therefore for the ring in question the 
three moments of inertia A = B, and C are all known: 

A = B = im{a^ + V), and C = ^^(a^ + V). 

If we now spin the ring about its central axis, we will have 
something like a gyroscope. Let the angular velocity of the 
spin be co; then, according to the formula if = 7co (in our 
case = (7co), we can calculate the angular momentum H of the 
rotating body; also its kinetic energy 7(co^/2) (in our case 
= C(co^/2)); also the necessary spinning moment, which will 
bring the system from rest to the angular velocity (of spin) 
= 0), in, say, t seconds, the acceleration being constant and 
equal, say, a. The following easy reasoning will be applied : 
since the acceleration is constant, multiplying it, that is 
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by dt and integrating we have 

where A; is a constant of integration, to be derived from initial 
conditions; in the beginning of motion, that is when ^ = 0, 
there is no rotation, that is 



hence k = 0, and finally 



dt "' 



d<p 



but d(pldt is our (given) angular velocity co; therefore co =at, 
whence a = co// is the constant acceleration d^ipldf. Ac- 
cording to this the spinning moment which can, with this 
constant acceleration, bring the ring from rest to the angular 
velocity co in given time, must be = I{d^(pfdf) or, in our 
notation, = C{o)ft). 

It is of advantage to represent the angular momentum about 
its axis of rotation as a vector, directed along this axis, and 
in sui^ch a manner that the right-hand rotation of the body 
will advance the end of the vector like a cork screw. With 
this understanding we shall now explain the meaning of the 
theorem: the time-derivative of the angular momentum about 
any axis is equal to the sum of external moments about that 
axis (see (1) under Integral of Areas) that is 

about any axis whatever. Let us consider the same ring of 
the moment of inertia, C, the angular velocity (of spin) 
being = co, as before. The angular momentum H about the 
axis of rotation will be represented (fig. 10) by the vector ab; 
now, if the ring, together with its axis, be displaced into a 
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new position ac in some plane E-F, we can easily see that the 
projection of the vector on the axis F will be decreased by a 
small amount m; while the projection of the vector on the 
axis E, which was originally 0, will now be n; therefore n/dt 




Fig. 10. 

and m/dt will be the time-rates of change of the angular 
momentum-vector on these axes, that is, the time-derivatives 
of the angular momentum about these axes (taken at random). 

Hence an elementary explanation of the action of a gyroscope 
(or top): a mere translation does not cause any particular 
effect because this does not alter the projection of the angular 
momentum upon any axis; but a rotation about an axis A, 
for instance, will immediately be felt by the observer: even a 
small rotation about A will change the projections of the 
angular momentum upon any axes such as E or F, and, as a 
consequence, external moments will arise, acting about these 
axes; the end A will have a tendency to rise, the end B to go 
down, precisely as if there were another moment applied about 
an axis P-Q, parallel to E. 

(c) Body with one point fixed. (Consult Bowser, Anal. 
Mech., pp. 493-505.) In order to understand what follows 
the reader should thoroughly master the following remarks: 
(1) A rigid body having one point fixed can be displaced into 
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any desired position by means of a rotation about some axis, 
passing through the fixed point (Euler's theorem); the proof 
is left to the reader (Hint: Describe a sphere about the fixed 
point; mark two points corresponding to the initial position; 
also two points corresponding to the final position of the body, 
upon the sphere; try to find a general method by which the 
first two points can be displaced into their new position). 
In view of this theorem any motion whatever of a body 
about a fixed point may be regarded as a rotation about some 
continuously changing axis (we assume it changing, since if 
it were stationary, this case would degenerate into the one 
just considered under 6), and with a definite (for that moment) 
angular velocity co. 

This axis is called the axis of instantaneous rotation, corre- 
sponding to the time t. It is well to note that the rotation, oj, 
itself is finite, although considered only during a very short 
period of time, after which it may change in both axis and 
magnitude. Now the rotation co, being a quantity given in 
both magnitude and direction, is a vector, that is a directed 
quantity, and as such it may be represented, to a certain 
scale, by a portion of a straight line upon its axis; also it may 
be resolved into components, coi, C02, C03, upon any three fixed 
axes, such as x, y, z. These component rotations simply 
mean that three axes, fixed in a body, coinciding at the time t 
with three axes fixed in space, would get away from the latter, 
diu*ing the next element of time dt, by the corresponding 
amounts due to coidt, c»)2dt, co^dt (this is a so-called rotation of 
axes about themselves). 

Suppose now that we have a point M (fig. 11), belonging to 
the body, and given by its coordinates x, y, z. It will tend to 
describe a circle about the instantaneous axis co, although in 
reality it may describe only a very small arc thereof, owing 
to the motion of the axis itself. Nevertheless we can deter- 
mine the linear components (upon x, y and z) of the velocity 
of the point M, during this very brief interval of time di. 
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having"at]our disposal only these data: the coordinates x, y, z 
of the point and the corresponding components, wi, C02, C03, 




Fig. 11. 
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These are given by Euler's 
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— 0)lZ, 



— co2ir. 



The angles establishing the direction of the instantaneous 
axis by, or in terms of, its components coi, C02, C03 are found in 
elementary works on this subject; they are as follows: calling 
<^f jSj 7> the angles of the instantaneous ^xis with the fixed 
axes X, y, z, we have 

COi CO2 CO3 

cos a = — : cos p = — ; cos y = — : 

CO CO CO 

CO is evidently the resultant of its own components 



CO = l^COi^ + C02^ + COs^. 



This establishes the position of the axis with respect to any 
axes for which the components are given as coi, C02, C03). The 
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reader will note that the Euler's formulae can be solved for x, 
y and z, should this be required, hut these solutions would involve 
time-derivatives (velocities). This is very important. 

Having thus found the component velocities of any point of 
the body, we might derive a formula of the kinetic energy 
2m(i;^/2) of the body for the same instant. We shall do it in 
a much shorter manner, however, remembering that the radius 
of gyration (see under Motion of rigid body) is given by 

k^ = €? cos^ a + 6^ cos^ jS + c^ cos^ 7, 

where a, 6, c are the radii of gyration corresponding to princi- 
pal axes of inertia -4, B and C\ so that 

ME} = ^ cos2 a + J5 cos2 jS + (7 cos^ 7. 

On the other hand we know that the kinetic energy of a 
rotating body is given by I(co^/2), and furthermore that the 
cosines of a, jS, 7, are 

cos a = — : cos p = — ; cos 7 = — ; 

CO CO CO 



so that, finally, 

Miej- 



= r = i(^coi2 + Bco22 + Ccoa^) 



which is generally written 

2r= ^coi2 + Scoa^ + Ccos^. 

The reader will understand that heretofore the axes ar, y, z 
were taken at random; but in this last calculation of T they 
were selected along the principal axes of the body, in order 
to do away with the products of inertia, which would appear 
for any other choice of axes. 

Another very important vector will now be mentioned, 
which will probably be rather unfamiliar to the reader. 



Genekal Principles of Dynamics. 49 

Under the Integral of areas we had the expression 

d ^ f dy dx\ ^, ^r x^ 

which was interpreted as follows: The time-derivative of the 
sum of the angular momenta is equal to the sum of moments 
of the external forces; by the sum of angular momenta (or, 
which is the same thing, moments of momenta) we mean the 
expressions similar to 



^ ( dv dx\ 



Imagine, now, that we are considering a body moving in 
any manner about a fixed point; let us draw a set of principal 
axes of inertia through the fixed point, about which the 
moments of inertia will be A, B and C. The body, in general, 
will be rotating about some instantaneous axis and the com- 
ponents of rotation upon these axes will be say coi, C02, C03; 
by means of these it is easy to calculate the linear velocities 
dxidt, dyjdt, dz/dt, according to Euler's formulae given above. 
Substituting these values of dx/dt, dyjdt, dz/dt, into 

^ ( dy dx\ 
we have 

^^ (2/37~^j7) = (az^m{x^ + y^) — coiSm.TZ — o)2^rnyz; 

but Sm(a:^ + y^) is the moment of inertia, C, about the axis z, 
while Sm.Tz and Xmyz are products of inertia; they vanish 
because our axes are principal axes; we have, therefore, 



^ { dy dx\ „ 
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and two similar expressions 

^ ( dz dy\ . 

^ f dx dz\ _ 

The vector, whose components upon the principal axes are 
the corresponding sums of angular momenta, that is = Ccos, 
etc., is called the impulse-axis; in our future work we shall 
use this vector rather freely; it is ihfi resultant or total angular 
momentum. In other words we shall use the expression 
impulse axis precisely as an equivalent to total angular momen- 
tum; it will be denoted by P. To interpret the meaning of P 
graphically, let the point M (fig. 12), under conditions similar 



s 




Fig. 12. 

to those shown in fig. 11, be one of the points of the system, 
in temporary rotation, during a very short period, about the 
instantaneous axis co. The momentum (1) due to this in- 
stantaneous rotation will of course be = mv, or, what is the 
same thing, mcor. If at the fixed point, 0, we add two equal 
and opposite momenta, (2) and (3), we can readily see that 
the effect of the momentum (3) will vanish owing to the reac- 
tion of the fixed point, while the other two momenta, (1) and 
(2), will form a couple, of which the moment will be the 
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moment of momentum of the particle M in regard to the 
point 0; this moment of momentum will be numerically 
= mvX OM (momentum times the arm) and can be repre- 
sented (as would an ordinary couple) by a vector of a certain 
length, laid off upon a certain axis ^; this axis must be: (1) 
perpendicular to the arm OM, and (2) must be in the plane 
MOo), Since r = OM sin a, it will also be equal to OM cos 
(w, ^), where by w, ^ we mean the angle between these two 
axes. So that if the angular momenta are taken, about the 
point 0, of all particles of the body, we can, adding such 
vectors as ^, obtain their resultant which will be the total 
angular momentum about the fixed point 0. This is precisely 
what we have just referred to as impulse axis and have denoted 
by P. Now if 

^ = mvOM = mcar X OM = mwOM^ cos (co^), 

we can also say that 

SJ (which is = P) = S(mcoOi/2 cos (cof); 

the projections of this resultant moment of momenta upon 
the principal axes are, as we have said before, Ao)i, Bw2 
and Ccoz. 

It is of interest to note that the kinetic energy, which we 
know to be ^{Awi^ + Bco2^ + Cws^) is also 

' ^ (tp- '^mr^fJ' CO ^ „ 
= J 2" = -^— = 2 Smr^« 

CO CO 

= 2 l^mcjOM^ cos^ (co ^) = ^P cos (co^). 

In other words the kinetic energy (body with fixed point 
only) is equal to the half -product of the total angular momen- 
tum (impulse axis) by the instantaneous angular velocity and 
by the cosine of the angle between their axes. 

Taking again the general expression of kinetic energy 
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2T = A(a-? '\- B(a^ -\- Co3^y let us differentiate it partially, 
with respect to each of the variables coi, C02, C03, so that 

i^T f^T f^T 

2 J- = 2^coi; 2 T— = 2BC02; 2 J- = 2Cco3; 

OCOi 00)2, 00)3 

multiplying the first of these expressions by coi, the second 
by C02 and the third by C03, and adding the results we have 

coi|^+ co2|^+ co3^ = Au)i^ + Bw2^ + Cwi" = 2r, 
acoi 00)2 ocoz 

which is another expression for the same kinetic energy. Com- 
paring the coeflScients of the like varables we have 

dT ^ dT ^ dT 

Awi = ;r— ; Bo)2 = t— ; Co)z = t— ; 
ocoi 00)2 O0)z 

that is, the projections of the total angular momentum upon 
any principal axis is equal to the partial derivative of the 
kinetic energy taken with respect to the velocity component 
upon that axis. 

Remembering the very important theorem (see under 
Integral of areas) that the time rate of change of the total 
angular momentum about any axis is equal to the sum of 
external moments about that axis, we shall consider the 
causes to which such change of angular momentum may be 
due: there can be but two causes: (a) The vector itself, being 
not constant, but variable (with the time), may change by a 
small amount, so that its projection, say Ao)i at the instant t 
may vary by the amount d(Acai) in the next instant, so that 
the time rate of change will be (d/dt){Ao)i); and (6) there is 
also rotation about the instantaneous axis, in which the whole 
body is taking part, including its principal axes as well as all 
vectors associated with them. We may express the time rate 
of change of the coordinates of any point x, y, z, upon any 
axes (otherwise called simply velocities) ; by Euler's formulae, 

dx 
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etc. ; so that the projection of the end of the vector P (which 
is the total angular momentum), will experience the effect of 
this double change in the very brief element of time, dt 

^ {Ao)i) + C02Z — CO32/. 

This, then, is the total time-derivative of the total angular 
momentum, and as such, is equal to the corresponding moment 
of the external forces; in other words, calling the external 
momenta about the axes x, y, z, simply L, M, N (to simplify 
writing) we have 

^^^coi + W2Z — CO32/ = L; 

and similar expressions for other axes, y and z. Now A is 
constant, so that the differentiation extends only to coi; and 
then, instead of any arbitrary point we will consider the end 
of the vector P, because it is its motion upon the axes that 
interests us, so that, substituting the coordinates of the end 
of P (that is, the projections of P itself upon the axes z and y 
which are: Bco2 and Ccoz); we finally obtain 

^^'+(C-JB)co2co3 = L; 

the two other similar expressions for the other two axes are 
derived in the same manner: 

B ^^ +{A- C)<^sc^i = M, 

C^+{B- ^)C01C02 = N. 

These are the famous Euler's equations of the motion of a 

rigid body with one point fixed. They give, for every instant, 

the time variation of the angular velocity components coi, 

C02, C03, about the principal axes, in terms of the components 
5 
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of external moments on these axes. In other words, if the 
external moments can be found about the principal axes, 
passing through the fixed point, we can find the rotation com- 
ponents coi, W2, W3 on these axes (and, therefore, the instan- 
taneous axis); or, again, if the latter components are given, 
we can find the moments which are producing them. 

Remembering that the angular momentum is the product of 
the moment of inertia by the angular velocity (see Bowser, 
Anal. Mech., p. 454) and that the components of the angular 
momentum upon the three principal axes are Acai, JBco2, Ccoz, 
we have the following pretty explanation of the meaning of 
Euler's equations (Bouasse, Mecanique, p. 563) : 




a ^6 ag,< tfi 



Fig. 12a. 

The velocity of the tip of the vector /x, representing the angular 
momentum, is equal to the vector representing at that moment the 
applied couple (Fig. 12a). 

Taking for example the first Euler's equation, let us find the 
velocity of /x as projected on the axis x. Such a velocity would 
consist of three parts: 

1. A(dwi/dt); this would be the result of the change of the 
vector Acoi with the time, and would be represented, say, 
by ab. 

2. Owing to the rotation of the vector Ccoz under the action 
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of the component W2 there would be a second component 
CCO2CO3, represented, say, by cd, 

3. Under the action of the component wa the vector JBco2 
would likewise contribute a component JBCO2CO3, shown by fc 
and directed backward (since all rotations are clockwise, 
looking toward the center). Hence Euler's first equation, 

.4^+(C-JB)co2C03 = L. 

The other equations can be explained in the same manner. 

The integration of these equations is rather diflScult and we 
shall limit ourselves to the case when there are no external 
moments, that is when L = M = N = 0. This can happen 
only in the following cases: (a) When the particles of the 
system are under action of no forces; (6) when the particles 
of the system are subject to some attractive or repelling 
forces, emanating from the fixed point, 0; (c) when the parti- 
cles of a system are subject to mutual attraction according to 
any law, or (d) when the system is under the action of gravity 
and its center of gravity is at the fixed point, 0. 

Let us consider the last assumption, that is the case of 
gravity being the only external force and the center of gravity 
being at the center 0* There will then be no external moments ; 
that is L = M = N = 0; so that Euler's equations will be 

^ ^ + (C - 5)C02C03 = 
B~^+(^-C)c03C0l=0, 

C^+ (JB - ^)coico2 = 0. 

One of the first integrals is easily found by multiplying 
these equations respectively by coi, C02, C03, and adding 

Aci)idcjoi + Bo)2do)2 + Cwzdo)z = 
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or, integrating, 

A(ji^ + -Baj2^ + C(jiz = const. = A. 

But this = 2r, or the kinetic energy which is, therefore, 
constant as it should be since the outside forces perform no 
work. This is the integral cf kinetic energy. 

Again, multiplying the same equations by Acoi, JBco2, Ccos 
and adding we have 

or, integrating, 

AW + BW + C^cos' = const. = P. 

Remembering that Ao)iy JBco2, Co)3, are the projections of the 
resultant moment of momentum (or impulse axis) upon the 
principal axes, we see that P is the square of the impulse axis 
itself, and that in this case it is constant, as it should be, 
since the outside moments L= M = N = 0; and the resultant 
angular momentum remains unchanged. The above expres- 
sion is the integral of areas for our case; we shall return to 
Euler's equations when studying the motion of a rigid body. 
8. Relative motion. It is useful to think of a coordinate 
system, in general, as of a mere temporary scaffolding, by 
reference to which motion can be identified by the observer 
or investigator; the axes possess no character of permanence 
and can be drawn in any convenient manner; moreover, in 
the most general case, the system, together with the axes to 
which it is referred, can have a bodily motion, translatory or 
rotational or both, in relation to some other set of axes* 
Therefore the observer, having determined the motion in 
relation to the original (so-called "reference'') axes, but not 
having taken into consideration the motion of the reference 
axes themselves, would be surprised to see that the results 
of his calculations might be altogether different from what was 
actually observed. The earth, for instance, possesses a 



General Principles of Dynamics. 



57 



number of very intricate motions in space, whose effect is 
hardly felt in our practical life; and yet a gyroscope or a 
pendulum can easily be made, which will actually reveal at 
least one of such motions, namely the rotation proper of the 
earth. In general, therefore, we shall consider two sets of 
axes (fig. 13); in the first place the reference axes x, y, z, to 




Fig. 13. 

which the motion of a particle (or system) can be referred; 
and another set of immovable axes, X, Y, Z. The motion of 
the given system in relation (fig. 14) to the reference axes, 
X, y, z, is called relative motion; the motion of the reference 
axes themselves, in relation to the fixed, immovable set. 




Fig. 14. 
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Xy Y, Zy is called "mouvement d' entrainement/* for which 
there is no English equivalent; we shall refer to such motion 
as MAR (motion of axes of reference) ; and, finally, the result 
of these two motions would be the true motion of the system 
(or particle) in relation to the fixed axes, which is termed 
absolute motion. So that we have to deal with paths, veloci- 
ties and accelerations, each of which may be (1) true, or 
absolute; (2) MAR, or of reference axes; and (3), relative, 
that is true only in relation to the relative axes, x, y, z. In 
this connection it is well to remember that, in general, relative 
equilibrimn will mean absolute motion; absolute equilibrium 
will mean relative motion. 

Very often the motion is such that the origins of both rela- 
tive and absolute system of axes coincide; in this case we have 
merely rotation of one of the sets, x, y, z, in relation to the 
other X, Y, Z, In a more general case, however, the relative 
(or reference) system has both translatory and, at the same 
time, rotational motion in regard to the absolute system, 




Fig. 15. 



X, Y, Z. If, for instance (fig. 15) Jf is a particle, either taken 
by itself or belonging to some system, its relative path (that 



% 
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is relative to the reference axes, x, y, z) may be such as Jkf-5; 
at the same time the MAR motion of the axes proper would 
bring this same point into D; while the resulting, absolute, 
motion (that is referred to the fixed axes X, F, Z, will result 
from these two motions and wijl be M-C. On the other hand 
another viewpoint can be taken: a particle M has absolute 
motion M-C; how will this motion appear to an observer, 
connected with some reference system x, y, z, which itself 
possesses a certain, known, motion in regard to X, Y, Z1 
The latter question is considered the real problem of relative 
motion, although both viewpoints are closely connected. 

In view of what has just been said, the main problem of 
relative motion can be framed in the form of the following 
question: Assuming a certain system, referred to some axes, 
X, y, z, to be under the action of certain forces and also assum- 
ing that the motion of such a system can be determined in 
relation to these axes, according to the rules of dynamics, 
already established; will there be any change in that motion, 
if the axes themselves, together with the system, are moving 
in relation to another set of fixed axes, X, Y, Z1 The answer 
is: the motion will be entirely different; it will be as if 
two new forces had been added to the already applied forces, 
and the motion then determined according to principles estab- 
lished above. The observer belonging to the original (refer- 
ence) system, x, y, z, will record precisely this, corrected, 
motion; what these two additional forces are will appear 
presently. 

We shall first establish the simple relations, connecting 
together the absolute and relative velocities, and also the 
absolute and relative accelerations. If, during a very short 
interval of time, A^, a particle, in its relative, apparent, 
motion, moved from M to B, and, if, during the same period, 
the whole relative path moved, in its MAR motion, to a new 
position, 2, so that the particle M would be displaced from 
M to D, then the resulting position of the particle at the end 
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of the period A^ will be C, found from a parallelogram, built 
on these two displacements, relative, M-B, and MAR dis- 
placement, M-D; so that the distance actually passed by M 
will be M-C. In other words, speaking geometrically, 

(MC) = {MB) + (MD), 

where the brackets indicate that the addition will be not 
algebraic, but meant in the same sense as when adding forces 
according to the law of triangle of forces. Dividing by A^ and 
passing to the limit we have 



-(f)='-(^) 



+-(^) 



but these limits all represent the time-rates of certain dis- 
placements, that is velocities of certain motions: 

lim ( -Tj ] = Va = absolute velocity; 

lim ( -^ 1 = Vmar = vclocity of the reference axes; 



and 



lim • ( -7-7- 1 = Vr, is the relative velocity; 



so that finally 

(Va) = (Vr) + (V)mar, 

that is the absolute velocity is the sum of the relative and 
(mar) velocities. From this, of course, the relative velocity 
will be (Vr) = (Va) — {v)mar all of which can well be repre- 
sented by the triangle of velocities (fig. 16). Having this 




Fig. 16. 
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construction in mind, we can easily derive the analytical 
values of velocities by projecting them upon the desired direc- 
tions (compare Bowser, Anal. Mech., p. 23). 

The problem of finding the absolute acceleration in relative 
motion is much more intricate. It will be remembered that 
the usual definition of acceleration (fig. 17) is as follows: 




Fig. 17. 

Let, at the time t, the velocity of a particle M be represented, 
in magnitude and in direction, by Vi; after a short interval of 
time. At, the new position of the particle will be M\ and the 
new velocity will be, say, V2; if we draw M-b equal and parallel 
to V2, then a-b will represent the variation of velocity during 
the interval A^; the limit of the ratio ab/At, of the variation 
of velocity to the corresponding time is what we call accelera- 
tion of the particle corresponding to the time t. (It is really 
the velocity of velocity, in other words the first time-derivative 
of velocity; or, the second time-derivative of space, dv/dt or 
d^s/dt^.) Its direction is along the acting force and its magni- 
tude is F/M, the acting force divided by the mass of the 
particle. 

To return to our problem of relative acceleration. Let 
(fig. 18) the relative path be (1) at the time t; let the bodily 
motion, MAR be of a translatory nature, so that after a 
short interval of time, A^, the new position of the relative 
path will be (2) parallel to the former position (1); if at the 
time t the relative velocity is p and the (mar) velocity is r, 
then, after the interval A^ the new relative velocity will be, 
say, q, and the new (mar) velocity will be 5; so that the corre- 
sponding accelerations will be, as explained above, ar and amav 
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Constructing an acceleration triangle, in precisely the 
same manner as we would a force triangle, we can say that 
under these conditions the resulting (absolute) acceleration 




Fig. 18. 

is the geometric sum of the two accelerations, of relative 
motion and of (mar) motion of the relative axes 

(tta) = (ttr) + {Omar) 

(where the brackets are meant to show that the addition is 
geometric, not algebraic) ; this, indeed, is a simple rule, but it 
applies only to problems in which the (mar) motion is trans- 
latory; such problems are not general and not interesting; 
in general the (mar) motion is both translatory and rotational 
and for this reason a further investigation must be made. 

Let the (mar) motion consist of a certain translation and 
of a certain rotation about some instantaneous axis (fig. 19) ; 



> 




Fig. 19. 



in other words let the relative path MMq be moved bodily 
into PMx, and then receive a rotation in the direction w about 
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the axis P-0, so that the final position of the path will be 
JPM2. A displacement of this double nature is the most 
general that can be imagined. (Owing to the very short 
period of time during which the displacement takes place, 
there can arise no question as to the order in which these 
displacements are to follow; the result will be the same what- 
ever the order.) Knowing how to express the space traveled 
by means of acceleration and time (Bowser, Anal. Mech., 
p. 235: 8 = ^ff) we can write M1M2 = iaddty (where a*, is 
some sort of acceleration, producing a given displacement 
M1M2 in given time dt). On the other hand the radius of 
rotation can be expressed from the triangle POM\] in other 
words MiO = Vrdt sin a, and the small arc M1M2 will then 
be M1M2 = v^t sin a-wt; comparing with M1M29 just found 
above, we have Oc = 2vr'0)*sina. This is an entirely new 
acceleration (called Coriolis^s acceleration) which gives us the 
desired correction. We can express it in the following easy 
manner: in order to find Coriolis's acceleration, project the 
relative velocity upon any plane perpendicular to the instan- 
taneous axis of rotation and multiply by 2w; finally turn 
through 90° in that plane, consistently with the direction of 
mstantaneous rotation. The result will be the Coriolis's 
acceleration, which is also known as compound centrifugal 
acceleration. Indeed by projecting the relative velocity on 
any plane perpendicular to the rotation we have Vr sin a; 
by multiplying by 2w we have the numerical value of ac; 
and, since such acceleration is acting from Mi toward Jkf 2, we 
see that it must be perpendicular to both the axis of instan- 
taneous rotation and to the relative velocity; hence — ^rotation 
through 90° in the plane perpendicular to P-0. This is the 
final solution of our problem: the acceleration of absolute 
motion is made up of accelerations of relative and (mar) 
motions and of the Coriolis's acceleration 

(tta) = (Or) + {dmar) + iflc) 
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(where the brackets indicate the geometric nature of this 
addition). If the value of the relative acceleration is desired, 
we have, from the same equation 

or, since the absolute acceleration is equal to the resultant of 
forces actually applied, R, divided by the mass of the particle, 
we have 

m 

Substituting in the above and multiplying by m, 

m(ar) = (R) — miomar) — m{ac); 

But the last two terms also represent certain forces (each 
being a mass X acceleration) ; let us denote, therefore, the 
force — m(amar) by Fmar and the force — m{ac) by Fc; when 
we will have 

miar) = {R) + (Fmar) + (Fc) = , Say, (Ri). 
Hence 



(«r) = 



m 



SO that the relative acceleration is directed along the resulting 
force and is equal to the latter divided by the mass. In other 
words the relative acceleration can be found precisely as if 
there were no (mar) motion of relative axes; except that the 
force is not merely the applied force, R, but consists of the 
latter plus two other forces, one of which equals minus mass 
times acceleration, due to motion of relative axes (— m{amar)) 
and the other of which equals minus mass times Coriolis's 
acceleration (— m{ac)). These are the two additional forces 
of which mention has been made in the beginning of the 
problem; introducing them we can altogether neglect the 
relative nature of the motion and apply our ordinary principles. 
The motion thus found will be precisely such as may be ob- 



General PMNaPLEs op Dynamics. 65 

served by an investigator, connected with the moving axes 
of reference. 

We have indicated by brackets the geometric nature of the 
additions but it is equally easy to put the matter in an analytic 
form. Let Xo, Yo, Zo, be the components of R, the resultant 
of the forces actually applied, upon the reference axes, x, y, z; 
also let Xi, Fi, Zi, be the projections of the additional force 
(Fmar) and X2, Y2, Zii those of the other additional force {F^ 
upon the same reference axes. Then we have merely three 
relations, one for each of the axes x, y and z, each expressing 
that mass times acceleration equals the sum of all forces upon 
that axis. That is, instead of 



we have 



m{ar) = {R) + {Fmar) + (^c) 



d^x 
m^ = XQ-\r Xi+ X2; 



and two other equations for the other axes. Integrating them 
we have six constants, which are supplied by the initial data 
of the problem. 

These two additional forces (Fmar) and (Fc) will vanish if 
the (mar) motion of the axes of reference is uniform, recti- 
linear and translatory, as in this case both accelerations (amar) 
and (ttc) are 0; the force (Fc) will vanish if there is no relative 
motion (vr = 0) ; if the instantaneous axis is parallel to the 
relative velocity, (a = 0); or if co = 0; that is if the (mar) 
motion of the axes x, y, z, is translatory (as we have already 
seen). These deductions can readily be made from the 
formula ae = 2vrCo sin a. 

The following typical example will show the application 
of the above: A material particle of mass m (fig. 20) is moving 
inside of a tube inclined at an angle a to a vertical axis, F, and 
in uniform rotation, co, about the latter. It is required to deter^ 
mine the motion. 
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Here the relative motion is that of the particle inside the 
tube, as viewed relatively to the tube; the (mar) motion 
reduces itself merely to the rotation of the tube (which is the 
relative path) about the vertical axis. Let s be the distance 
of the particle, at the time t, from the intersection with the 




>m(i)V 



Fig. 20. 

axis. The acting force is the weight mg; one of the two addi- 
tional accelerations to be found is minus the (mar) acceleration; 
since this (mar) motion in our case means simply rotation, 
we have that its acceleration is w^r (centripetal force) and this 
is directed toward the axis of rotation. The first additional 
force will thus be mwV and directed away from the center; 
the second (Coriolis's) additional acceleration acts in the 
following manner: it is located in the plane perpendicular 
to the axis; is tangent to the circle described by the radius r, 
and is directed forward, in the sense of rotation; therefore the 
force (Fc), which is minus mass times (cc), is also tangent to 
the circle but directed backward. They are both perpendi- 
cular to the relative path (tube) and therefore do not aflfect 
the motion (no friction being considered) in any manner; 
the latter force, however, (Fc), represents the action of the 
particle upon the tube. In order to determine the motion 
along the tube the most natural way would be to project all 
forces upon its direction; having taken care of the additional 
forces, we can consider the relative motion (inside the tube) 
as absolute and apply our ordinary methods. The accelera- 
tion along the direction of the tube will be cPs/dP, and the 
projected forces will be mg cos a and moci^r sin a or moj^s sin a 
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(since r = « sin a) ; therefore 



\ 



^* I 2 • 2 

m T^ = mg cos a + moi^a sm^ a, 
whence 

^* . 2 • 

-p = gf cos Qj + 0)^5 sm a. 

This is the final equation; in order to integrate it we first put 
it in the form ^ 

d? / <y cosa \ 2 • 2 Z' I flL^osa \ 

dp \ OCT sin^ a / \ w^ sin^ a / 



whence 



s + 4^?^ = ^e"' •'*"*+ 5e-"'«*°*, 



0)^ sin^ a 



which gives the distance in terms of time; and, diflferentiating, 

ds 

^ = w sin aC^e"'^*'*'* - Be'-^''''"'), 

which gives the relative velocity. These equations contain 
two constants of integration, A and B, which will be deter- 
mined from initial conditions. If, for instance, at the be- 
ginning of motion (when t = 0) the particle is at V and has 
no initial velocity, we have 

ds 
.0 = 0; ^ = 

whence 

g cos a 



A = B = 



2(j? sin^ a ' 



substituting into the expressions for a and isldi just found, 
we have all that is necessary to specify the motion. Should 
the projection be required, on the horizontal plane, of the 
path of the particle in its absolute motion, we have but to 
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note that the polar coordinates of the particle will be 

r = s sin a and = (at, 

which in the final equation for s will give a spiral. 

A most important expression will now be established for the 
absolute velocity in relative motion. We have seen that this 
is made up, geometrically, of the relative velocity and of the 
(mar) velocity, of the motion of the reference axes themselves. 
An analytical expression will be derived as follows: Since any 
motion whatever can be decomposed into translation plus 
rotation about a certain instantaneous axis, let us take for 
origin (fig. 21) any point in the moving body and let us 




Fig. 21. 

imagine the actual, absolute, motion resolved into two mo- 
tions: the (mar) motion of the axes x, y, z (together with the 
body), which is to be purely translatory; and the relative 
motion of the body about some instantaneous axis w. The 
latter motion will be the same as of a body with a fixed point 
(see under Body with one point fixed) , since we have fixed the 
coordinates, x, y, z, in the point of the body, so that there 
can be no other relative motion but rotation about some 
axis (w) through this point. Since the (mar) motion is trans- 
latory, the axes x, y, z, being at all times parallel to the fixed 
axes X, Y, Z, it is evident that the (mar) velocity of all points 
of the body will be the same, and will be equal to that, say. 
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of the origin, 0. Analytically therefore the projections of 
the absolute velocity on fixed axes, X, 7, Z, will be made up, 
(1) of the corresponding projections of (mar) motion, and 
these which are the same as projected velocities of any point 
will be designated by Vmarx, Vmary, Vmart] and, (2) of the corre- 
sponding projections (on the axes x, y, z, which are parallel 
to the fixed axes) of the velocities of the particle we are con- 
sidering, due to the rotation of the system about the fixed 
point 0. From the theory of rotation of the body whose one 
point is fixed the values of these velocities are (Euler's formu- 
lae) 

dx 

^ = 0)23 — cozy; etc., 

so that finally the absolute velocities on fixed axes, X, Y, Z, 
are 

Vax = "Omarx + ^2^ — 0)3!/, 

Vay = Vmary + (^^X — WiZ, (1) 

VaM = Vmart + COiJ/ — (a2X. 

These equations give the absolute velocities of a point, located 
by its coordinates, within the moving (reference) system a:, y, z, 
of which the motion is the same as that of a body with one 
point fixed. 

But we can imagine a still more general case, where the 
point itself will have a motion of its own within the moving 
system; let the velocity of this point on the moving axes be 
dx/dt, dy/dt, dz/dt Since the velocities of this same point, 
stationary within the moving system, and moving with the 
latter, are given by (1), we shall now have 

dx 

Vax = j7 + f^marx + (^2^ — 0)32/, 



dt 

dy 
dt 



Vay =17+ 'f^mary + CO^X — 0)12, (2) 



Vat "= ^ + '^mars + ^i^ — 0.^2^, 
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as projections of the absolute velocity of a particle, moving 
within a moving system of reference, a very broad expression, 
which we shall need in om* future work. 

To return to the equations (1), which are of exceedingly 
great importance, we can write them in the condensed form 

Let us suppose that instead of the arbitrary point 0, we have 
chosen, as the origin of moving axes, the center of gravity, G. 
Let us, under this supposition, write down three equations 
similar to (3) for each point of the system, square them, 
multiply each by the corresponding mass and add the results 
together; this will give 

Now the first term ll>m{va:? + Va^ + 'Ca^F) is evidently = Hiiwo^y 
«a being the absolute velocity of the center of gravity; likewise 

the term ^miVmarx + rimary + 'TlmaTz) = ^mVmarg^y Or, simply, 

Mvmarg^ whcrc Vmarg is thc total {inar) velocity of the motion 
of the center of gravity in its translatory movement; also 

-[(t;+(i)+(i)i--^' 

because dxjdt, dy/dt, dz/dt are merely projections of the velocity 
on the reference axes x, y, z (in our case projections of the 
relative velocity Vr). The last term can be rearranged as 
follows: Omitting the brackets we have three terms 

dx dy dz 

2Zi7)lVmarx "TT , ^^'t^'^mary TT } 2^inVmarz "77 ; 
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but, since the motion is translatory, Vmarx is the same for all 
particles, including the center of gravity, G (and the same 
applies to Vmary and Vmarz)'y SO that these velocities can be 
placed before the S sign: 

dz dy dz 

But, for the center of gravity, 'Emx = 0; Smy = 0; Smz = 0; 
(Bowser, Anal. Mech., p. 109) and therefore also 

dx dxi dz 

2"^^=0' 2m^=0; Sm5^=0. 

The whole last term of (4) is therefore 0, and we have 

or, dividing by 2, 

^ 2 " 2 "*"^ 2 • 

This theorem (called Koenig's theorem) is of greatest im- 
portance in practice and will be constantly used in our future 
work; it tells us that the total kinetic energy of any system is 
equal to the kinetic energy of its center of gravity (with the 
whole mass concentrated in it) plus the kinetic energy of the 
whole system, in its relative motion about its center of gravity 
(considered as fixed). Caution, — ^In calculating the kinetic 
energy of the center of gravity, mv^/2, absolute, not relative 
velocity should be used; it is geometrically compounded of 
relative and (mar) velocities. 

Let us take a few examples. 

1. A thin horizontal rod, of the mass m, can slide in a holder y 
which can turn about a vertical axis. The holder has negligible 
mass and there is no friction. Find the general expression of 
kinetic energy in term^ of the angular velocity of rotation (fig. 22). 

Let C be the middle of the rod (the center of gravity). 
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According to Koenig's theorem we have to calculate, (1) the 
kinetic energy of the middle point, C (the whole mass being 




Fig. 22. 

assumed concentrated in it), in its rotation about the vertical. 
This of course 

■"T"/ "2^"2" 'di^ 

(compare Bowser, Anal. Mech., p. 243, also, Bowser, Calculus, 
p. 178). To this we must add (2) the kinetic energy due to 
the motion of the rod about its center of gravity, in other 
words, to its rotation, of which the instantaneous value will 
be the same as d<p/dt, the rotation about the vertical support. 
(This is often misunderstood: a well-known theorem of kine- 
matics tells us that the axis of rotation can be transported 
(parallel to itself) to any point, provided that a certain trans- 
latory motion is added; but the angles of rotation will be 
precisely the same for all such parallel axes; hence, instead of 
rotation about V we may simply consider that about C. 
This remark will apply to almost every problem of this sort; 
we shall transfer centers of rotation freely, into any desired 
position.) 

Now, the kinetic energy of a rotating rod = | moment of in- 
ertia X the square of angular vel ocity, that is (ml24:)P{d<pldty 
(Bowser, Anal. Mech., p. 431); so that the total kinetic energy 
will be 

We note that this result is in terms of variables dr/dt, r and 
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d<p/dt, of which we know nothing, so far; it will be shown in 
the next chapter how, from expressions of T, found in the 
general form, similar to the expression just calculated, the 
equations of motion can be easily derived by Lagrange's 
beautiful method. 

2. Find the kinetic energy of a ladder (mass m) the top of 
which is resting on a smooth lamp-post, while the bottom is 
slipping, in some manner, on the icy sidewalk (no friction). 
(By saying that the ladder is slipping " in some way '' we 
mean to say that the motion of the bottom of the ladder will 
not be necessarily radial, away from the post, but may have a 
rotation, at the same time, about the post, as shown by the 
angle (fig. 23).) 

Here according to Koenig's theorem we have: (1) the center 
C of the ladder is moving on a sphere, described from with 




Fig. 23. 



the radius = 1/2 (check this statement by elementary methods 
of Anal. Geometry) ; therefore the velocity of C can be made 
up of very small displacement- velocities: along the meridian, 
Z/2 • d<pldt and along the arc of a parallel circle, (1/2) sin (p 
(dd/dt) ; hence the expression 



m 



im^'Mm- 



Now the relative motion of the ladder, that is to say, about 
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its center as fixed point, may be about a horizontal axis, if 
sliding outwardly, and, also, about a vertical axis, in its rota- 
tion about the post; so that we have to consider both rota- 
tions. Taking the ladder merely as a rod of the same mass 
and length, we have its moment of inertia = mP/12, there- 
fore this part of kinetic energy will be 

(Note that (dd/dt) sin (p indicates that the rotation about the 
post, dO/dt, has been resolved, along, and at right angles to 
the ladder; and only the latter component taken into con- 
sideration; draw a sketch, directing dd/dt vertically along the 
lamp-post and illustrate this remark.) 
Finally we have 



r-mh 



sm' 



<m 






(Remark, — Instead ot T = S(miJ^/2) it is customary to write 
2T = ^mv^, and we shall make frequent use of this notation; 
likewise, in dealing with rods, as was the case in the last two 
examples, the length is generally denoted 2a or 21 to simplify 
writing.) 

3. A rod (mass m, length 21) can slide, without friction, on a 
horizontal plane; its end A (fig. 24) cannot leave the axis x, 





while its end B is acted upon by a repelling force F, perpendicular 
to 0-x and proportional to BD, the ordinate of B, Find 2T 
in terms of (p. 
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This is a somewhat more intricate case. Let a and b be 
the coordinates of the center of gravity of the rod. The only 
forces applied to the rod consist of F and of corresponding 
reaction of the axis x; they are both perpendicular to the 
axis X, hence, with regard to this axis, the equation of motion 
can be written 

(Pa 
m iTo" = 0, or a = kt -\- n; 



dt 

where h and n are the constants of integration. This means 
that the motion of the center of gravity, parallel to the axis a;, 
will be uniform. The general expression of the velocity of C 
will be 

'-=(5)'+(fr=''+'-'^(f)'. 

(since 6 = Z sin ^), 

which gives the velocity of the instantaneous bodily motion 
of the rod. We now must calculate the kinetic energy due to 
the rotation of the rod about its center with the angular 
velocity d^/d^, which will he^{mP/3){d(p/dty and the total 
kinetic energy will be 



2T = mV^ + 



f(^)'=»[^+''«+-'^>('^)']- 

9. Euler's angles. We shall terminate this chapter by 
deriving Euler's angles, which will be often used in our work. 
We have seen (under Body with one point fixed) that the motion 
of a body about a fixed point resolves itself simply into a 
rotation, with a certain angular velocity, w, about some 
instantaneous axis, whose projections upon any rectangular 
axes, 0)1, 0)2, wa, are" connected together by the relation 

0)^ = 0)1^ + 0)2 + W3^. 

These projected rotations give us an idea as to the rotary dis- 
placement, which is going to take place, during the next 



76 



Lagrange's Equations. 



instant, about each axis; but they do not help us directly to 
characterize or identify the position of the body at a given 
time. Indeed, in order to do this we would have (1) to imagine 
a system of axes, x, y, z, say, the principal axes, fixed in the 
body, that is to say, moving with it; and (2) a system of axes 
Xy Y, Z, fixed in space. (3) We would then have to find the 
cosines of each axis x, y, z, with each of the fixed axes, X, Y, Z 
(by rules of analytic geometry), that is nine cosines for the 
whole set of movable axes x, y, z, and, finally, (4) we would 
derive expressions of time-derivatives of all these cosines 
which would enable us to watch the motion of the axes x, y, z, 
that is of the body, in relation to a fixed set of axes, X, Y, Z. 
All this can be done, in a much simpler way, by Euler's angles 
and the formulae that can be derived therefrom. 

According to Euler's method, in order to change the posi- 
tion of axes from any initial position, X, Y, Z, to any current 
position, X, y, z (fig. 25) only three rotations are necessary: 



>^i-. 




Fig. 25. 

(1) about Z, through an angle ^; (2) about the new position 
Xi-0, through an angle 0; and (3) about the new axis, z-0, 
through an angle (p. This completely identifies the new posi- 
tion of the axes, and therefore, the new position of the system; 
therefore, if we could watch the motion of the system by means 



I 
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of the components wi, W2, C03, which, through the change of 
the nine direction-cosines could tell us where the axes would 
be during the next instant, we certainly could do so, with 
much greater ease, by following the changes of Euler's angles 
of which we have but three. Indeed, the resultant of the 
three small variations of Euler's angles will give the body the 
same instantaneous twist as would result from the three com- 
ponent rotations wi, W2, W3. Our problem is then to replace 
wi, 0)2 and W3 by expressions involving only the angles xp, 
and <p, as well as their derivatives (angular velocities), dxp/dt, 
dOjdt and d<p/dt; in other words to follow the motion of the 
(moving) principal axes, along which it is customary to 
direct the moving axes x, y, z, entirely through the angles rp, 
and <p, reckoned from some fixed system of axes X, Y, Z, 

Now, since the instantaneous rotation w can be resolved 
into any desired components, let us resolve it into three com- 
ponents, along the axes 0-Z, 0-Xi and 0-z (fig. 26). This 




Fig. 26. 



means that, instead of old components of instantaneous rota- 
tion, wi, 0)2, 0)3, we shall have the following new components 
(which, combined together, will give, of course, the same 
effect, 0)) : About 0-Z • • • d\l//dt, about 0-Xi • • • dd/dt, about 
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0-z • • • d<p/dt. All we have now to do is to project dyl/jdt^ 
dO/dt and d<p/dt upon the moving axes x, y, z (remembering 
that the projection of w is equal to the sum of projections of 
its components). This gives 

# . . dd 

wi = -^ sin ^ sm ^ + ^ cos <p, 

dyp . dd . 

0)2 = -^ sm cos ^ — ^ sm <p, 

C03 = ^COS9+^. 

It is very advisable for the reader to derive these formulae 
several times; the mechanism of projecting the components is 
clearly indicated on the drawing; in practice, as has been said 
before, for axes z, y, z, we generally take the principal axes of 
the body, which is thereby fully identified in reference to some 
fixed axes, such as X, Y, Z, It is well to remember that 
Euler's angles, ^, and (p are positive in the direction indi- 
cated by the arrows; rotations in an opposite direction would 
be assumed negative. The reader will probably see at once 
why the elementary rotations drp/dt, dO/dt and d<p/dt have 
been projected upon the moving axes (or, rather, upon direc- 
tions, coinciding at that instant, with the moving axes, 
^f Vf z)f while it would be just as easy to project them upon 
the fixed axes, X, Y, Z; the real advantage is this: by pro- 
jecting the elementary rotations upon the moving axes we 
have means of locating, for any instant, the position of the 
instantaneous axis in relation to these moving axes; while the 
moving axes themselves are located at any instant, in respect 
to the fixed axes, X, Y, Z, by Euler's angles ^, and (p. In 
other words, as soon as Euler's angles are given (in terms of 
the time), the reader will picture to himself two things: 
(1) by means of these angles, the position, for any time, t, of 
the moving axes (and therefore that of the body itself), can 
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be instantly established; and (2) within these moving axes, 
the instantaneous axis of rotation, can be established by its 
three projections wi, C02 and C03 on these same axes, which, 
furthermore are given in terms of these same Euler's angles. 
The latter are, therefore, suflScient to completely identify the 
position of the body at any given moment of time. The 
reader will no doubt understand that in problems involving 
spinning bodies (tops, gyroscopes, etc.) the axis of spin and 
the instantaneous axis co are two entirely different things; in 
using Euler's angles for such problems it is customary to 
direct the Euler's axis z along the axis of spin, in which case 
the angular velocity d<p/dt is the constant velocity of spin; 
while the instantaneous axis w may be entirely outside of the 
rotating body; these two axes can coincide only in a very few 
problems, involving, for instance, spinning of a top about its 
steady vertical position. 

By way of illustration we shall take the following very 
simple exercise, although it is of extremely great importance: 
Supposing that a point is given in reference to a new system, 
^f Vf 2, by its coordinates, x, y and z; and that this new system is 
characterized by Euler's angles, ^, and <p in relation to the old 
{or fixed) system, X, Y, Z; how can the position of the point be 
found in direct reference to the old coordinate system? In other 
words, how to find X, Y, Z, in terms of x, y, z, and of the 
angles, rp, and ^? In order to find these we have only to 
find expressions connecting together the two sets of coordinates 
and then to solve them for X, Y, Z. The reader is advised 
to draw two sets of axes to a large scale and then to proceed 
as follows: the new system was obtained from the old one by 
three successive rotations rp, 6 and <p; now, after the first 
rotation, rp, when the axis X became Xi, any point, Xi, yi, 
in the system thus obtained, would be expressed as follows 
in the X, Y, Z system: 

X = Xi cos ^ — 2/1 sin ^; Y = Xi sin ^ + 2/1 cos yp; 
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where Xi and Yi are the temporary coordinates of the point 
. on the axes Xi and Yi (this is evident from plain geometry). 
After the second rotation (through 0, about 0-Xi), the axis Z 
took its final place z and the axis Yi moved up to, say, Y2; 
under these conditions any point Y2, 2, Xi will be expressed 
as follows in the system Fi, Z, Xi: 

Yi = Y2 cos — z sin 9; z = Y2sm6 + z cos 6, 

where Xi remains the same as before, but Y2 and z are the new 
coordinates in the system just obtained. Another, final, turn 
will now be made about z, through the angle <p, which finally 
gives any such point as x, y, z, in the system Xi, Y2i z as follows, 

X\ = X cos <p — y sin <p\ Y2 = x sin <p -{- y cos <p 

{z being the same as before). Eliminating Xi, Fi, and Y2 
from these three pairs of equations we have 

X = a:(cos ^ cos ^ — sin ^ sin (p cos 0) 

— 2/(cos ^ sin ^ + sin ^ cos <p cos 6) + z sin 6 sin xp; 

Y = a:(cos ^ sin ^ + sin <p cos xp cos 0) 

+ yisin ^ sin ^ — cos ^ cos <p cos 0) — z cos ^ sin d; 
Z = a; sin ^ sin + 1/ cos ^ sin + z cos d. 

These equations in themselves are not interesting, but the 
following deduction should be well retained by the reader for 
his future work: it is possible to find the coordinates of any 
point for the old system, solely through the Euler's angles and 
the new coordinates; but, in finite form, that is involving no 
derivatives (compare this with the result of our attempt to 
express a point in terms of the rotations, wi, 0)2 and C03 (see 
under Euler's formulae); it was shown then that x, y, z can 
also be given, in terms of wi, W2 and W3; only such expressions 
would necessarily involve velocities, time-derivatives). 

Euler's angles are of great importance, therefore another 
example will be given: A solid body, with one point fixed, is in 
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motion which can be described as follows; the instantaneous axis 
is rfiomng in space, and the projection of the constant angular 
velocity, 12, about it {as laid off on that axis) upon a certain 
straight line, D, belonging to the body (and passing through the 
fixed point), is constant and = co; furthermore, the projection 
of the angular velocity, Q, upon a plane perpendicular to D, is 
turning, relatively to the moving axes, x, y, z, with a constant 
angular velocity = — co. Required to determine, in terms of 
time, t, Eukr's angles, with which the system can be located in 
relation to some fixed axes. 
Let the line D (fig. 27) be the axis z of the moving system, 




X, y, z (solid with the body) ; also let 12i be the projection of Q 
upon the plane perpendicular to D; let the moving axis x coin- 
cide with the projection of Q at the beginning of the motion 
(that is when i^ = 0) ; finally let the fixed axes represent the 
position of moving axes at the beginning of motion, ^ = 0. 
From Euler's formulae just found we have 

# . . do 

coi = 12i cos o)t = -jT sm d sm ^ + "iT cos <p. 



dil^ . 



de . 



C02 = — Qi sin w^ = IT sin 6 cos <p — -jrsin <pj 



(1) 
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We have thus resolved the instantaneous velocity into its 
components along the moving axes, x, y, z; as bl matter of 
fact the instantaneous velocity was given to us not by itself 
but by its projections; nevertheless, the formulae (1) will 
enable us to find d\l//dt (from the first two, multiplying by 
sin (p and by cos (p respectively, and adding) 

dxf/ __ Qi sin ((p — a)t) 

df " ^bl ^^^ 

and also to find dd/dt (multiplying the first by cos <p, the second 
by sin (p, and subtracting) 

de 

-jT = Qi cos ((p — o)t). (3) 

Also from the last equation (1), 

d<p ' dxl/ - ,. ,_.. _ . . - cosff ... 

CO — ^ = ^ cos * = ('ro°^ (2)) = Qi sm {<p — o)t) -^-g (4) 



whence (multiplying (3) by (4)) 



de 

T. COS 



COS (<p — cot) f dip \ \_^^ _ A 

sin {ip — o)t) \dt J sin 6 ' 

which immediately integrates into a logarithmic expression, 
from which 

sin {(p — cot) sin 6 = const. = (5) 

(the constant being = 0, because of our choice of coordinates, 
in view of which, when t = 0, 6 = do = 0). The equation (5) 
can be satisfied only, when (1) 6 = which means merely 
rotation about a fixed axis (z the same as Z); or, (2) when 
<p = cot; this means d(p/dt = co and (see (4)) d\l//dt = 0; 
whence rl/ = const. = (in view of our conditions). Finally 
from (3) we have dO/dt = 12i and d = ilit it appears, there- 
fore, that the line D remains in the plane X-Z (since ^ = 0), 
while revolving about Y with constant angular velocity Qi 
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(since d = Qit); and, meanwhile, the body is revolving about 
that axis with a constant velocity w (since (p = cot). 

We give one more example of the application of Euler's 
angles: Consider the behavior of a gyroscope on the surface of 
the earth (fig. 28). 

Let Z be the vertical axis and X, Y, Z, the fixed axes; is 
one of Euler's angles, z — the axis of spin, and x and y are 




reference axes moving together with the body; let co^i, C0e„ co^g 
be the components along the same axes, x, y, z, of the rotation 
of the X, Y, Z system proper (owing to the rotation of the 
earth); the rotation of the gyroscope proper is cos; the rotation 
o)e of the earth gives an additional component co^g = — coe cos a. 
All acting forces (gravity and centrifugal) have resultants 
passing through the fixed point, so that their moments are 0, 
and we can apply Euler's third equation 

C^+ (5-^)coia)2 = 

C is here the moment of inertia of the gyroscope about its 
axis of spin; A = B are two other moments of inertia, equal, 
arid at right angles to the axis of spin; so that 



C-T.ii^z — coe cos a) = 
or 0)3 — coe cos a = const. If 0)3. and ao are the initial values 
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of 0)3 and a, we have 

W3 = cosq + a>«(cos a — cos ao) 

in the absence of external work the kinetic energy (remem- 
bering that 2T = Acai^ + Bo)2^ + Cw3^ and that ^ = J5) is 
constant: 

A{0>i' + 0)2^) + C0>z' = ^(0)1,2 + C02,2) + Ca)3,2 

substituting wi and a>2 in terms of d and ^ (see traru^ormation 
formulae), as well as C03 just found, we have, neglecting ta^ 
which is small, 

4(i)*+(fy-'] , 

= 2Cco3oCo«(cos ao — cos a) + A(ju)\^ + 0^20^). 

Now, in order to utilize the fact that there are no external 
moments about the du-ection N, fixed m space, we can state 
that the sum of the projections of the impulse axis on this 
direction is constant; these projections (instead of Acaiy JB«2, 
C(t>z)y are -4(coi + (^c^ cos (iV, z), etc., so that 

A{(jii + a>e,) cos (iV, x) + 4(w2 + coe,) cos (iV, y) 

+ C'(a>3 + cocg) cos (iV, z) = const. 

But we can consider, at some given instant, the axis x as coin- 
ciding with the projection of N upon the equatorial plane of 
the gyroscope; in that case 

cos (iV, x) = sin a; cos (iV, y) = 0; cos (iV, z) = cos a 

also 

coei = "~ «e cos {Nx) = — CO siu a\ coe, = 0; 

Weg = — ■ C0« COS a 

so that, finally 
A{(ai — o)e sin a) sin a + c(w3 — coc cos a) cos a = const.; 

or, substituting C03 from above 

A((ji)i — We sin a) sin a + 0(0)3^ — coe cos ao) cos a = const. 
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In order to find this constant let us re-write this equation for 
t = 

A(o)i sin a — coi^, sin ao) + -4coe(sin^ ao — sin ao) 

+ C(o)Zq — coe cos ao)(cos a — cos ao) = 0; 
whence 

^(coi sin a — wi^ sin ao) = (cos ao -- cos a)[^coc(cos ao 

+ cos a) + c(co3o "" ^« cos ao)]. 

This equation together with the integral of kinetic energy, 
found above, are the solutions of the problem. This is a 
rather abstract problem introduced here solely for the purpose 
of illustrating the application of Euler's angles. In what 
follows several easy problems will be given on gyroscopic 
motion. 

{Caution, — ^The reader will understand that conceptions 
like force, energy ot work are really never relative; for instance, 
the kinetic energy can be calculated merely as Sm(»o^/2) 
where Va is the absolute velocity; although Koenig's theorem, 
which gives such an easy working rule, might be misinter- 
preted to mean that there is such a thing as relative kinetic 
energy in itself; we shall often refer to kinetic energy calculated 
relatively to a certain point or axis, as if they were fixed, but 
this should not mislead the reader; while of course relative 
velocity or relative path are perfectly real conceptions.) 



CHAPTER II. 

« 

Lagrange's Equations for a Particle. 

The position of a free material particle in space has so far 
been characterized by three coordinates, x, y, z, referred to 
three rectangular axes. When asked to determine the motion 
of a particle, we had to form the equations of motion and to 
derive from them the expressions of coordinates in terms of 
the time. Yet, mention has been made (see under Coordinates 
of rigid body) of a manner in which the number of coordinates 
can be reduced, owing to constraints; each constraining equa- 
tion reducing the number of coordinates by unity; thus for 
example to determine the position of a particle in space three 
coordinates are needed, but, to determine the location of a 
point on the surface of the earth we need only two charac- 
teristics, longitude and latitude, the third characteristic being 
implied in the condition that the point is on the sphere, given 
by the equation of the latter. To make one step further, a 
point on the equator can be characterized by only one co- 
ordinate, such as an angle from some fixed point on the 
equator. What became of the other two characteristics? 
They have not been lost sight of or dropped altogether; they 
are implied in the doubly constraining condition; that the 
point is on the surface, and that it is in the equatorial plane, 
intersecting that surface. Hence it is easy to see that every 
additional constraining condition reduces the number of 
necessary characteristics or coordinates by one. 

The characteristics which specify the position of a particle 
upon the given curve or surface are called generalized or inde-^ 
pendent coordinates. They are independent in the sense that 
in forming equations of motion involving these coordinates 
we can totally disregard the equations of constraints, unlike 

86 
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in methods dealing with ordinary position coordinates. 
These independent coordinates can be of quite different 
nature from the usual position coordinates: for instance^ 
in the example of the particle on the surface of a sphere, we 
have two independent coordinates, longitude and latitude, 
both being angular; in the example of a particle on the equator 
of a sphere we have only one such independent coordinate, 
the azimuthal angle from a given meridian; but we could 
likewise take as independent coordinate, say, a linear distance 
along the equator, from a given point. So that we can take 
as independent coordinates any characteristics, angular or 
linear, provided that they fully describe the position of the 
particle on its constraints. It is important to realize, how- 
ever, that these independent or generalized coordinates, al- 
though less in number than ordinary position coordinates (one 
or two instead of three) lead to no ambiguity whatever in 
characterizing the position of a particle. In fact any one of 
the ordinary position coordinates, x, y, z, can at all times be 
expressed in terms of the new, independent, coordinates; and 
such expressions, if substituted in the equations of constraints, 
will render them identical. Take, for instance, the equation 
of the sphere a^ + y^ + z^ = r^; it is easy to see that the 
coordinates of any point, x, y, z, can be given in terms of the 
generalized coordinates, longitude and latitude, by the equa- 
tions (fig. 29) 

X = r cos p cos q, 

y = r cos p sin q, 
z = r sin p, 

where q is the longitude and p the latitude; substituting these 
equations into that of the sphere we have r^ = r^, an identity. 
To recapitulate; generalized coordinates imply independ- 
ence of the constraints; they are less in number than the posi- 
tion coordinates, x, y, z, yet the latter can always be expressed 
through them, and such expressions must necessarily satisfy 
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the constraints. . The number of these generalized coordinates 
is called the degree of freedom, as we have already seen. In the 
case of a constraining curve, which is given as the intersection 
of two surfaces, we have only one degree of freedom; the 




Fig. 29. 



particle can move only along the curve and we need only one 
characteristic, that is, only one generalized coordinate, to 
completely locate its position on the curve; for instance in 
the case of a particle moving on the meridian of a sphere we 
had two equations Q!^-{-y^-{-:s?--r^=0 (the sphere) ; in 
general, say, fi(x, y, z) = 0; and z = (the equatorial plane) ; 
in general, say, /2(ir, ?/> 2) = so that only one generalized 
coordinate was necessary; we saw that it can be either an 
angle, reckoned from a given meridian, or a distance along 
the equator, from a certain initial point. And if we knew 
how to express this only variable, the generalized coordinate, 
which we shall call q, in terms of time, by some such equation 
as g = f(t), we would know all about the motion of the particle, 
since the latter is fully located by this one variable, q. In 
other words, while in three-coordinate system we would have 
to know three such expressions, connecting coordinates, 
Xf y> 2, with time, t, x = Fi(t); y = F2{t); z = Fz{t), with our 
new method we need only one equation q = f(t) ; afterward. 
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if desired, we can express the old position coordinates in terms 
of q, therefore obtaining the same equations, x = Fi(t)y etc. 
Limiting ourselves, for the present, to the case of a particle 
moving along a curve given by two equations, fi(x, y, z) = 
and f2(x, y, z) = (which indicate permanent constraints, 
containing no time, t), let us see how the introduction of q, 
our new independent or generalized coordinate, will aflPect the 
shape of the fundamental equation of dynamics (from (1), 
under D'Alembert's principle), which, written for one particle 
of mass m, is 



+ (z-mg)6. = 



or 



Xdx+Ydy+Zdz=m(j^dx + ^6y + ^,dzy (1) 

the first part of which represents virtual work (see under 
Virtiml work), and was denoted by &W. Knowing that 
Xy y and z can ultimately be expressed in terms of our new 
variable q (see example above), by some such equations as 

x=^F,{q)', y=F2(q); z=F^(q); (2) 

let us differentiate these equations with respect to time t; 
this will give velocities dx/dt, dyjdt, dz/dt in terms of q and 
of derivative dq/dt; in order to simplify writing we shall adopt 
the following notation; the time-derivatives will be indi- 
cated by (0 for each variable; so that dx/dt, dyjdt, dz/dt, dq/dt 
will be denoted simply as x', y\ z' , g'; and the values of these 
time derivatives will be (from (2)) 

, d-Fi dFi dFz 
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or simply 

/ 5^ / , dy , , dz , 

*=a^^' V=rq^' '=d-q^- 

(The reader is advised to review his stock of knowledge of 
partial differentiation; see Bowser, Calculus, pp. 120, 122, 125.) 

These values of x\ y\ 2', if substituted in the general ex- 
pression of kinetic energy 



r-^ or 



will give 

2T = mix'^ + y'^ + z'^) (4) 

which otherwise may be written thus (in view of (3)) 

-="'[(f)'-(f)'-(fl]'^. 

which is a quadratic function in q\ whose coefficient, [ ], 

contains only q or its functions and constants; in other words 

2T= *(g) X g'' (40 

which means that T is a function of q of constants and of the 
square of g' but contains no first power of q\ 

Remembering that q itself is a function (as yet unknown) 
of the time, as is its derivative g', we can differentiate (4) with 
respect to q' (it will be partial differentiation, during which 
all variables other than g', or its functions, remain constant) 

— ^m\^x-,+ y—,+ z^); (5) 

and also, with respect to q (partially), 

dT ( ,dx' ,dy' , ,dz'\ ,.^ 
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To return to our fundamental equation (1): remembering 
how the equation (3) under Virtual work was derived, we can, 
by a similar application of Taylor's theorem, obtain the 
following (from (2)): 

dx dy dz 

in view of which 

XSX+ YSy+ZSz=[xf^+ r|+z|)sg; 
or, putting 

we have simply 

SW = XSx + YSy + ZSz = Qdq. 

Here SW is virtual work, while Sq is the generalized displace- 
ment, and perfectly arbitrary, of our new independent coordi- 
nate; therefore Q cannot be other than some sort of a force, 
corresponding to such a displacement, since the product of 
the two represents work. Let us call Q generalized force. Its 
expression can be transformed as follows: In 

let us substitute, instead of the forces X, Y, Z, their equiva- 
lents, m{d?xldP), m(d?y/d^), m{(Pz/df), and extend our simpli- 
fied notation, so that 

this will transform (7) into 



m 



('"!-:+ ^'l+^'l)-=«- <« 



> 
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But this can be simplified still further; in the first place 

which can be immediately verified by diflPerentiation of 
x''(dxldq) with respect to time [the same remark applies to 
y''(dy/dq) and z'^dzjdq)]] on the other hand, diflPerentiating 
(3), partially, with respect to q' we have 

dx' dx ,^^^ 

with two similar equations in y and z. Also, diflPerentiating 
(3) with respect to q we have 



dx^ _ d^x 
dq dq 



= —2Q' (11) 



(because, since q = /(O, q' is a function of t only), and two 
similar equations in y and z. Also, evidently 



|(|)=ii9'.thatis = i^ (12) 



d^x Sx^ 

dt\dqj-d^^^''^^^^'''=d^ 



(from (11)). Therefore, substituting (10) and (12) into (9) 
and, finally into (8) we have 



m 



V dq^y dq^^ dq) 



d ( ,dx' , ,dy' , ,dz'\ 

( ,dx' , ,dy' , ,dz'\ ^ ,,„, 
-m{^^'-+y'^ + ,'-j=Q, (13) 



which in view of (5) and (6) means 



d fdT\ dT 
Equations of this type are called Lagrange's equations of 
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motion; it will presently be shown that for the motion on a 
surface, where two independent coordinates are necessary to 
specify the position of a particle, there will be two such equa- 
tions, one for each independent coordinate; also in the next 
chapter it will be explained how this same equation can be 
extended to a system. 

The deduction of these equations is not diflScult and not 
unnatural: all that we have done was to substitute new nota- 
tions in the equation expressing that kinetic energy equals 
work done; this, in the first place, gave us a new expression 
of kinetic energy; on the other hand the " work done " has 
the same value, no matter what the choice of coordinates; it 
is always the product of force by displacement, and if, for 
some reason, we have chosen, instead of, say, 5a;, a peculiar 
displacement dq, we must expect that the corresponding 
(generalized) force will also be of a special nature, subject, 
however, to the condition that its product Q8q is equal to the 
kinetic energy. 

Now, we can imagine several kinds of displacements: 

1. Linear displacement, the product of which by the corre- 

sponding force gives work. 

2. Angular displacement, the product of which by the 

corresponding moment gives work. 

3. Increase of volume, the product of which by the corre- 

sponding pressure gives work. 
Etc. 

So that what we called generalized force may be not only 
force, but moment, pressure, etc., according to the choice of 
the independent coordinate; it will be easily found from the 
conditions of the problem. If there are no external forces Q 
will be = 0; if there are external forces, the force Q must be 
found such, that it will correspond to the virtual displacement 
of our new (independent) coordinate (force for linear, moment 
for angular displacement, etc.). As a final result of these 
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transformations (and this was our primary object) we have 
only one equation of motion (14), instead of three; it can be 
written down as soon as the expression of kinetic energy T 
has been found; and has nothing to do with the constraints. 
Before illustrating the application of this equation by an 
example we can observe the following: in order to simplify 
the derivation of these equations we have made the assumption 
that the constraints fi(x, y, z) = and f2{x, y, z) = were 
permanent, that is that they contained no time, i. In general, 
however, the constraints may be variable in either shape or 
position, and therefore, as a general rule, their equations will 
be, say, fi{x, y, z,t) = and fzix, y, z,'t) = 0, and conse- 
quently the equations (2) will also contain the time, 

X = Fi(q, t); 2/ = F2{q, t); z = Fz{q, t). 

By way of illustration imagine a particle on a sphere of varying 
diameter — a soap bubble. Here it is not possible to charac- 
terize the position of a particle by longitude and latitude 
alone, but time t must also be specified. We can easily prove, 
however, that moving constraints do not alter the form of 
Lagrange's equations; the kinetic energy will no longer be a 
function of g'^ alone, but will have another term depending 
upon g'; but, although the equations (3) will have an addi- 
tional term: 

, dx dx , 

the equations (10) will still hold true, and in deriving the 
equations (12) we shall simply have 

dx' _ d^x d^ , 
dq""didq^d^^' 

but also {d/dt){dx/dq) will now be (d^x/dtdq) + (d^x/dq^)q' that 
is = dx'/dq, from what has just been obtained; so that the 
equation (14) will remain unchanged. 
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Having thus extended the method to the case of variable 
constraints we shall make another remark regarding the 
problem involving two parameters. The proof just given con- 
templates a particle with only one degree of freedom, that is 
to say, the motion is limited to a curve, given by two equa- 
tions. Suppose now we have only one equation of constraints, 
J{Xy y, z) = 0; this will represent a surface, so that two param- 
eters or independent coordinates will be necessary to locate 
the position of the particle on that surface. In this case we 
can apply the same reasoning, except that two generalized 
coordinates will mean two generalized forces, and then there 
will be two Lagrange's equations, one for each independent 
coordinate. This can be applied to any number of coordi- 
nates, as we shall see in deriving Lagrange's equations for a 
system; but a particle can have only one, two or three degrees 
of freedom (in the latter case it is absolutely free of con- 
straints). We shall rapidly review the derivation of equations 
similar to (14) for two independent (generalized) coordinates 
(such as latitude and longitude, for instance), which we shall 
call qi and 92* 

Let the constraining equation be f{x, y, z) = 0, and assume 
that the coordinates «, y, z, can be expressed in terms of the 
two independent coordinates thus: 

X = Fi(qu q2); y = -^2(^1, ^2); z = ^^3(^1, ^2). (2') 

Differentiating and substituting in the equation of kinetic 
energy we have 

that is a quadratic function, homogeneous in the generalized 
velocities {qi and g2')> where E, F, G, contain only generalized 
coordinates or constants. 
Differentiating 

2T= m(a;'' + / + 2'') 
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with respect to gi and qi; also to ^2 and ^2', we have two 
sets of equations like (5) and (6). But in this case the virtual 
displacement is 

_ dx ^ , dx ^ 

(by Taylor's theorem) with similar expressions for dy and dz. 
Finally we have 



+(^:i+>'i+^af>'- 



Here again we have generalized displacements dqi and 6^2; 
so that {8W being work) the coeflScients of dqi and 6^2 must be 
generalized forces. Let 

dqi dqi dqi 

so that bW = Qidqi + ^28^2; repeating the method used above 
we obtain two equations like (8), one for qi and one for ^2- 

Simplifying by the introduction of equations exactly similar 
to (9), (10) and (12) we finally have 

dt\dqi') dqi'^'' 



d /ar 



These are Lagrange's equations of motion for a particle with two 
degrees of freedom. For an entirely free particle we would 
have three such equations with three generalized forces, 
corresponding to the three generalized displacements; that is 
one equation for each independent coordinate. It will be 
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seen that they would be precisely alike in form and all derived 
from the same T; but the partial differentiation is performed 
with respect to the three different q's or g"s and the corre- 
sponding generalized forces, Q, will be different for each of 
these equations (15). 

The case of varying constraints, for the reason given above, 
does not alter these equations; in that case instead of (2') we 
shall have x = /(gi, 72, t), etc. 

Example, Form equations of motion of a free particle in a 
horizontal plane, under action of any force, P, 

Here the only constraining equation is that of the plane, 
2=0; the particle evidently has two degrees of freedom, and 
we shall choose as independent coordinates what otherwise 
are known as polar coordinates, r and <p; so that gi = r, and 
q2 = (p- The kinetic energy in polar coordinates (see under 
Integral of kinetic energy) will be 

dr^ + r'^dip^ 



2T = m 



de 



or, in our simplified notation, 2T = m{r^^ + ^^' ). 

Differentiating 2T with respect to r and r'; also to <p and ^', 
we have 

dT 2 dT dT dT 

Substituting in Lagrange's equations 

d {dT \ dT d (BT \ dT 

dt\d^' ) ~ d^^- ^' """^ ^\d^')~d^,-^' 

we have 

mr" — mr<p^^ = Qi, and m j. {r^(p') = Q2. 

These are the general equations of motion, except that Qi 
and Q2 (generalized forces, corresponding to our generalized 
displacements Sr and 5^) must be found from the original 
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conditions and substituted in these equations. Resolving the 
given force, P, into its (generalized) components, one, R, 
directed along the radius vector, the other, X, at right angles 
to the radius, we can observe the following: The force R is 
evidently = Qi, because it is directed along the radial dis- 
placement Sr, multipUed by which it gives Qiir, that is work; 
the force X itself, if multiplied by the corresponding displace- 
ment, gives XS<p, which is not work, therefore A" is not Qt; 
we can see at once that Qt will not be a force but a moment, 
since the corresponding displacement is angular, and it is 
necessary to have a moment, in order to have tcork, when 
multiplied by an angular displacement; so that Qt is not N, 
but Xr. 

The final equations of motion are 



Example. Form general equations of motion of a free particle 
in space. 

We take for independent coordinates r, and <p (see fig. 6 
and corresponding remark on kinetic energy in polar coordi- 
nates). Here we have 

and there will be three Lagrange's equations for the three 
independent coordinates, r, and ip, 

mV - r(e'2 4. cos2 ^.^/2)] ^ q^^ 
m I j^ (r^eO + r2 sin e cos ff-i^'M = Q2, 

m T. {r'^ cos^ 0'(P^) = Qz- 

It remains only to substitute the corresponding values of the 
generalized forces Qu Q2 and Qz, in accordance with the 
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external forces as given in the problem. Suppose we have a 
force F which can be resolved into the following components: 

R along r, 

M perpendicular to r in the meridian plane^ 

P perpendicular to R and Jf , 

then Qi will be the force R itself, because the displacement ir 
along it will be linear; but Q2 will be Mr (because of the angular 
nature of the displacement q^y or bS) ; also Qz will not be simply 
P but Pr cos By because the angular displacement i(p will 
result in a certain work iW done by a force P acting on a 
lever r cos 6. 
Finally 



m 



"'[i''f+'"»°'"»»«(f)"]-^'' <"' 

m -t; r^ cos^ B -jr = Pr cos d. 
at at 

With a little practice the reader will be able to write down 
equations like (16) or (17) at once. Now the integration of 
equations thus obtained, or, in fact, of all but a very few 
simpler equations of dynamics, is very difficult, and often 
altogether impossible; but this difficulty is of a mathematical 
nature, while the formation of the equations is the real problem 
of dynamics. And in this connection Lagrange's equations 
are often of the greatest help. 

An interesting case presents itself when the acting forces are 
the derivatives of some common force-function XJ. This 
symbol may be used in a broader sense than before, and may 
contain not only the coordinates but the time, t, explicitly 
(which the potential function we had before, could not), so 
that 

-^=g^ ^=f^ -=v^ <«) 



^ 
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If there is such a function, we can introduce (18) into (7), 
and obtain 

^ dx dq dy dq dz dq dq 

so that instead of Q we will have, in equations (14), simply the 
partial derivatives, such as dUjdq, with respect to the 
corresponding independent coordinate. The same remark 
will apply to a greater number of parameters, so that, for 
instance, the equations (15) will then be re-written thus 



dt\dqi') dqi~ 



dq 



dt \ dq2 ) dq2 dq^ 



and so on. 



The general rule of procedure when a force function can be 
found is as follows: first, we select the independent coordinates 
in such a manner that the x, y, z, coordinates can be expressed 
by the former in connection with the equations of constraints, 
which may also contain time; this will give equations (2) 
or (2'). The next step is to form T by the proper substitution 
of x\ y', z' y derived from (2) or (2') ; the function T will thus 
contain time, ^, as well as the independent coordinates and 
their time-derivatives; after this it remains to form the ex- 
pression of generalized forces Q, or to substitute proper values 
of (2) or (2') into the force function t7, which will then con- 
tain only independent coordinates and, possibly, i. Thus 
we have everything required for writing down Lagrange's 
equations, one for each independent coordinate. 

A few specially selected examples will now be given. 

Example 1. Find equations of motion of a particle, of mass m, 
• moving without friction on a straight line, inclined at an angle a 
to a vertical axis (fig. 30) and rotating about the latter ivith the 
angular velocity co. 



• • 
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Comparing this with the equation of pendulum (Bowser, 
Anal. Mech., p. 458) we conclude that the motion of the par- 
ticle will be pendular about ^ as a lowest point, and the time 
of an oscillation will be = tt/co. 

Example 3. A particle starts with an initial velocity ^o (fig. 
32) along a circumference whose radius increases with the time; 
find the motion. 

This is a typical case of varying constraints (where x = /(g, t), 
etc.). If the initial radius is a, the constraining equation will 




)/. • 



Fig. 32. 

be a^ + y^ = a\l + kty^ where A: is a constant coefficient. 
Of course z = 0; also, introducing polar coordinates, r and 
<p, X =^ a(l + kt) cos <p and y = a(l + Jet) sin <p; differ- 
entiating with respect to both the independent coordinate 
and the time, t, we have 

a;' = — a(l + kt)<p^ sin <p -{- ak cos <p, 
y' = a(l + kt)ip' cos ip-\- ak sin ^, 

which in 2 r will give 

2r = ma\{l + A:OV'' + 1^]=^ mv^ 

(incidentally it follows that a^[(l + ki^ip^'^ + A:^] = i)^ at any 
time, or at the beginning of motion, when ^ = 0, a^^' -\x\k^= V(? 
whence aV' = «o* — lA^). From 27 we have 



dT 



d<p 



7 = ma%l + ktf<p', 
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so that 

-T. [ma^{l + kty<p'] = 0, or ma^{l + k{fip' = const. = c. 

Let us determine the constant c for the beginning of motion, 
that is when ^ = 0; we have 7ad?<p' = c; but a^^' = ^o — <^, 

so that a^p' = V^i^o^ --«^P; whence t^miV' — mal/i?o^ — Jfc^ = c. 
Finally our equation becomes 

d(p 



ail + kifip' = Vn^ -o%2 or ^= ^ 



dt {1 + kty ' 

where 



-4 = 

a 

Separating the variables we have, 

dip = A 



(1 + ktf ' 
and finally 

But for < = 0, ^ = 0, and Ci = (-4 /A:), so that 



^■"i V 1 + A:^/ A 1 



A:^ ^^ 



+ A:^ l + Ar 

Substituting this into x and y, formed in the beginning of the 
problem, we have 

x= a{l + kt) cos YJ^t ' 
y = a(l + kt) sin j-:^-^^ , 

these are the final equations of motion (and completely express 
the coordinates in terms of the time). 
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Remark 1. Due mention has been made of the ease with 
which generalized forces can be found if there is a force func- 
tion U; the latter may be a potential function, that is a func- 
tion of the coordinates only; or, it may contain the time; but, 
in any case, its partial derivatives, with respect to any co- 
ordinate, must give corresponding forces, so that 

^— =Z. etc.; 

instead of generalized forces Q, etc., we then simply have 
dU/dq, etc. 

But if there is no force function, the generalized forces Q 
can be found from general considerations as follows: We 
have seen (equation (7)) that 

so that it is only necessary to calculate the partial derivatives 
dx/dq, etc. from the equations (2) or (2'), express x, y, z through 
the independent coordinates, and substitute in (7); the result 
will be Q for that specific displacement. Partial differentia- 
tion means assuming other generalized displacements, all = 0, 
and the virtual work, Q5q, deriving only from the displacement 
we have under present consideration, say 8q; similarly for Q2 we 
shall find dx/dq2y etc.; and so on; X, F, and Z, being of course 
the projections of external forces upon the axes, x, y, z. 

Remark 2. Should it be desired to calculate the reactions 
of the constraints, this can be easily done by referring to equa- 
tion (4) under D'Alembert's principle, where Xi, X2, etc., are 
Lagrange's multipliers, or as they are often called, constraint- 
coefficients. The working of this method will be illustrated 
presently. 

Remark 3. We have seen, under Integral of kinetic energy, 
that, for constrained motion, if the constraints are permanent 
or independent of the time, the integral of kinetic energy 
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r = f/ + A holds as true as for non-constrained motion; this 
integral can also be shown to be an immediate consequence, of 
and therefore an equivalent to, of any of Lagrange's equations. 
Considering the simplest case of motion, with one degree of 
freedom, we have 

dT 
dq 

multiplying by q' we have 

,d {dT\ ,dT ^, 

But q\d!dt){dTldq') is the same as (d/dl)(,q'{dT/dq')) - 
q"{dTldq'), as can immediately be verified by differentiating 
q'(dT/dq') with respect to the time; therefore 

^ , d / ,dT\ „dT ,dT 

Remembering the remark made under (4') we can dijflferentiate 
2T = ^(q)q'^, with respect to g', and the result (which holds 
true for permanent constraints only) will be 



whence 



2^=20((7).g', 



(z'f^=0(^)-g'' = 2r. (21) 



From the same expression, for 2T, not containing time ex- 
plicitly, we have 

dT_dT ,,dT „ 

dt~dq^ '^dq'^ ' 

which substituted in (20) and then in (19) will give 

dt^'^^' dt dt '^'^ ^dt 



I 
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whence dT = Qdq, which is the same as (2) under hdegral of 
kinetic energy. 

In the case of several degrees of freedom, we have, in the 
same manner, 

which likewise leads to the integral of kinetic energy. If 
there is a force function U of qu q^, qz, of which Qiqi + $2^2 
+ Qiqs is an exact differential, dU, then we simply have 
dT = dU, so that T = U -{- h; this is the case when there is a 
potential function, U{x, y, z), not containing the time. Hence * 
in problems involving the potential function it is perfectly 
legitimate to replace any Lagrange's equation by its equiva- 
lent, the integral of kinetic energy. 

The most diflScult equations can thus be done away with 
and the remaining set, including the integral of kinetic energy, 
will then be the final equations of motion. 

Examvle. By way of illustration of these remarks let us 
take the case of a particle moving, under gravity, on a sphere 
of constant radius r. 

We had a similar problem (see equations (17)) for a more 
general case of motion; as before, let us take as independent 
coordinates the latitude d and the longitude <p. Through 
these, our original coordinates can be expressed as follows: 






X = r cos d cos <p; y = r cos d simp; z = r sin 6; 

arid, r being constant, the kinetic energy will be 

2T= mir^e'^+r'cos^e.ip'^). 

Writing Lagrange's equations, one for d and one for <p, we 
have 



mr{d'^ + sin 6 cos fl^'^) = 0, and m t. (cos^ B^p') = $; 



dt 
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where and $ are generalized forces which we shall now 
proceed to determine (see remark 1). The axis z being 
directed upward, the potential function will be 

J7 = — mgz = — mgr sin 0; 

forming dU/dd we have 

-^ = 6 = — mgr cos 0; ' 

also 

dU ^ ^ 

so that the final equations will be 

-p-+ sin S cos d(p'^ = — ^f cos 6, and t; cos^ ^ ^ ~ ^' (^2) 
From the latter we have 

which is the integral of areas for S = (that is for the plane 
x-y). 

If instead of a potential function U we had a force whose 
projections upon the axes were X, Y, Z, we would find and $ 
as follows: In view of the remark 1 we would have 

but 

X = r cos 6 cos (p; y = r cos 6 sin (p; z = r sin 6; 

so that 

■— = — r sm cos (p; T7: = — r sm fl sm <p: ^ = r cos 8. 
ad 06 du 

Substituting in the expression of we would find it in terms 
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of the (known) projections X, Y, Z, and of the parameters 
and <p; $ would be found in the like manner. (Instead of 
this we might at once have written down and $ in terms of 
the components of F, tangent to the sphere in the meridian 
plane, M; and along the tangent to the parallel circle, P (see 
equations (17)). But our object was to indicate a general 
method for obtaining the Q's.) 

According to Remark 3, we can replace the first equation 
(22) by the integral of kinetic energy, which is T = U + h. 
In our problem J7 = — mgr sin and 

T = Yir^0'^+r'cos^0'ip'^) 

so that the final equation T — U = h will be 

f^0'^ + ^ cos2 0'ip''^-\- 2mgr sin = 2h 

(where A is a constant); this equation, together with the 
second equation (22), will determine the motion; it is only 
necessary to substitute ^' from the latter into the former and 
to integrate the result. Such integration is not easy, although 
it can actually be performed. 

The reaction of the constraints (in our case, of the spherical 
surface on which the particle is constrained to move, or else 
of the rod (length = r), fixed in the center of the sphere, which 
can perform precisely the same function as would the sphere), 
according to Remark 2 can be found as follows: The result of 
integrating the equations of motion (22) will give us and <p 
in terms of time, t; as soon as we have and <p, we can substi- 
tute them in the expressions of x, y and z, which will give 
X, y and z in terms of t] then the constraints' coefficient X can 
be calculated from the general equation 

d?x ^ , , df 
(with similar equations for y and z), where X is the projection 
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of the acting force and / is the constraining equation (in our 
case the equation of the sphere, a^ + y^ + s?-— r^= 0). 
Since gravity is the only acting force, we have X = 0; F = 0, 
and the only projection is Z = — my ; so that X can be derived 
from any of these three equations 

(1/X u/V U/Z 

mj^ + \x=0; m^ + \y = 0; m-j^ — g + yiZ^O. 

If the motion were on a curve instead of a surface, we would 
have two constraining surfaces; and two coeflScients Xi and X2 
could be derived from any two equations of motion; but in 
our specific problem a much easier method can be used. 
Mention has already been made (see the problem of a particle 
sliding ojBf a cylinder, under Integral of kinetic energy) of the 
fact that the pressure exerted by a particle describing a curve 
is equal to the centrifugal force minus the radial component 
due to the applied force; in our problem the centrifugal 
force is mif^/r, and t^ can be derived from the equation of 
energy, T = U + h, or 7f^ = 2{h '— gz); on the other hand it 
is easy to see that the component of the weight, mg, upon 
the radius will be = mgz/r so that the reaction 

P = - (2/i - Zgz), 

h being a constant. 

Remark on relative motion. We need not dwell on the 
subject of relative motion in studying the application of 
Lagrange's equations to a single particle, for the reason that 
here the kinetic energy, which is the first item to be calcu- 
lated, can easily be found. Indeed, we have seen that Cori- 
olis's acceleration, which is the most characteristic attribute 
of relative motion, does not perform any work, being at right 
angles to the relative path; hence it does not appear in the 
expression of kinetic energy, which is based upon absolute 
velocity, or, according to Koenig's theorem, can be derived 



[ 
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as a sum of the kinetic energy that would be due to (mar) 
motion alone, and of the kinetic energy that would be due to 
relative motion alone. In studying the motion of a system 
(two or more particles, rigid body, etc.) which will be done in 
the next chapter, we shall have several new principles to estab- 
lish in connection with relative motion; but the word system 
implies such things as principal axes, moments of inertia, etc., 
all of which will complicate our problem quite considerably 
and, incidentally, make it most fascinating. But in dealing 
with a particle we have no such complications; for instance, 
the problem of a particle, moving freely in a circular tube, 
revolving in a horizontal plane, which is really a typical 
problem on relative motion, was solved by us without any 
reference to relative motion: having selected cp as our inde- 
pendent variable, we wrote down the expressions of absolute 
coordinates, x and y, in terms of (p. After this, an easy 
dijflferentiation and substitution in the expression of 27 gave 
us the desired expression of kinetic energy, and the solution, 
which showed that the motion of the particle would be pendu- 
lar, gave us merely the dependence of <p upon the time, which 
made it purely relative, that is such as it would appear to an 



'V^i 



^"^.-^ M 




Fig. 33. 



observer, connected to (and traveling with) ' the moving 
system, or ring. But of course it is just as easy to treat the 
problem from the standpoint of relative motion: we can find 
the absolute velocity Va as the geometric sum of two velocities. 
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tmar and Tri after which Va can be substituted in mCo^/2; thus, 
IV is evidently the velocity along the circumference, = r X ang. 
vel. = r<p'; while the Tmar is the velocity of the point M (see 
fig. 33) rotating about with the angular velocity a> of the 
reference system (of the circumference itself). Adding these 
we obtain v,? and therefore the kinetic energy. We shall ask 
the reader to check this in accordance with the drawing and 
to show that the kinetic energy so obtained is precisely the 
same as before. 



CHAPTER III. 
Lagrange's Equations for a System. 

By a system we mean two or more particles, connected in 
some manner; or, more particularly, a rigid body, a great 
multitude of particles subject to the one essential condition, 
permanence of form; or, a combination of such bodies. In 
deriving Lagrange's equations for a system we shall be gov- 
erned by the great principle of virtual work, in view of which, 
if we know, for any virtual displacement, the total amount 
of work done by all forces, we have all that is necessary for 
finding the motion. Indeed, the conception of work enables 
us to fonn the expression of kinetic energy of the system, and 
from the latter, as we have seen before, and shall again re- 
establish presently, in a different manner, the equations of 
motion can be derived. Of course it is almost evident that 
what has been said regarding a particle can be extended to a 
system, but we shall derive Lagrange's equations without 
reference to what has been said in the last chapter. However, 
in following the mechanism of derivation of these equations, 
it is well to constantly bear in mind Lagrange's fundamental 
idea: the general equation of dynamics can be split into two 
parts, one of which, representing virtual work, consists of a 
coefficient (say Q) multiplied by the virtual displacement. 
Of course for the same virtual displacement, dq, the coefficient 
Q must be the same, regardless of what particular system of 
coordinates x, y, z, we have chosen. Therefore, the other 
side of the equation will likewise be independent of the co- 
ordinates in the same sense. But, kinetic energy can be 
expressed in terms of g, • • • , etc., q\ • • • , etc. (and of t, some- 
times) and is independent of the choice of coordinates x, y, z. 
Hence, the natural conclusion that the first member of the 
equation just referred to must be made derivable from T. 

t 113 
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Lagrange's beautiful transformation shows how this can be 
done. 

Consider a system consisting of n separate particles whose 
coordinates are given in rectangular coordinates by xi, yi, zi, 
X2, y2, Z2, ' ", ocn, Vn, ^\ each particle is actuated by forces, of 
which let the components upon the axes be Xu Yu Z^ Xz, 

Y29 Z29 • • • , Xn, Yn, Zn* 

Then the fundamental equation of djTiamics, for a system 
comprising all such particles, will be 

+ (z-mf)«.] = 0; 
which can be re-arranged thus: 

S(Z5a: + Y5y+ Z6z) = Sm ( Jfix + ^Sy + ^5z^ . (1) 

Adding to this equation the equations of constraints we have 
all that is necessary for finding the motion of the system. 
Let us select k independent variables gi, q2, •••, Qk (this 
number k will be equal to 3w — m, where n is the number of 
particles, and m of constraining equations, as we have seen 
before) by which any coordinates xi, yi, zi, X2, ^2, Z2, • • • , etc.> 
can be expressed in some such manner as 

xi = /i(gi, ^2, •••, qk, t), 



yi = Mqi, q2, 


"f qk, t). 


2i = Mqi, q2. 


'", qk, t), 


^2 = Mqi, q2. 


• • ', qk, t), 


2/2 = fbiqi, 92, 


'", qk, t), 


Z2 = Mqi, q2, 


•••5 qk, t), 







(2) 



Xn = Szn^2{qi, q2, " ', qk, t), 
y« = /3n-l(?l, ?2, •••,qk.t), 

2n = faniqi, Qi, '••, qk, t). 
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In other words every coordinate (three for each particle) can 
be expressed, in terms of the independent variables, and this 
may and, generally does, also involve time, t (that is, including 
the most general case, that of movable constraints). If the 
independent coordinates have been chosen in accordance with 
the equations of constraints, the latter will be satisfied by a 
substitution of (2) therein. Our object is to transform (1) 
to new variables q^ q2, • • • , etc. The first part of (1) can be 
easily transformed thus 

2(Z5a: + YSy + Zdz) = Qidqi + Q^hq^ + • • • Q*5g*. (3) 

Indeed, from (2) we have (by developing a: + 5a: by Taylor's 
theorem in the same manner as has been done before), 

dzi dxi dxi 

dj/i ^ , dj/i ^ , , dj/i . ^^ 

. dz\ - , dz\ - , , dz\ . 

_ dxi . , dxi _ , , dxt . 

5a:2 = ^^5gx + ^5g2+---+^^5gfc, 

^ ^2/2 . , ^2/2 . , , 3^2 . 

hzi = ^ 5gi + HT 5^2 + • • • + J7 Sgjfc, (4) 



dq\ ^ dq2 ^ dq 



k 



dXn^ , dXn^ , , dXn 

oqi ^ dq2 aqk 



r 
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^ 



substituting these in (3) and factoring by 8qi, Sqi, • • • , dqt we 
have 

^(Xbx + YSy + Zdz) 

/ xr 5ari ayi dzi , ydx2 dy^ 

^ V-'d^. + ^^d^+ ^^d^i+ ^*d7x+ ^^d^x 



+ 



■+^-^+''.^+^.l>. 



+(^'l'+'''l'+^. 



3zi 

3^2 



+ 

or, denoting the coefficients of hq\, hq^, • • • , S*?*, by Qi, Qj, 
• • • > Qt. we have simply 

SCJTSx + Yhy + Z5z) = QiScji + q,ihqi + • • • + Qt^?*, 

that is (3), where 

dxi ai|^ dzi dx2 .7^. 

^' = '^'di,+ ^'d^,+ ^'6^2+ '■' + ^'d^.' ^^^ 

We know that (3) represents virtual work and we had a 
special notation for it, 5W. (It is well to observe, however, 
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that, in general, there is no such function as W, the exact 
variation of which is (3), so that dW is merely a symbol, 
introduced in order to simplify writing.) 

From (3) it follows that all terms like Q • dq are elements of 
work; and, whatever the nature of q, if dq is a generalized 
displacement, its coefficient, Q, must be a corresponding forces 
which we shall agree to call generalized force. (As we have 
seen this may be not only force but moment, pressure, etc., 
according to the nature of the independent variable q itself.) 

It only remains to transform the second member of (1), so 
as to make it involve only the kinetic energy T, In the first 
place, using the notation cPx/df = a:", etc., and considering 
(4) we can rewrite the second member thus (where the sum- 
mation includes all g's, from qi to qt inclusive) 

Sm(."|+y"|+."|)5,= 5PF. (6) 

On the other hand the expression of kinetic energy is 

2T= 2m(a;'' + / + z'') (7) 

but, from (2) we have (remembering that the variables are 
<lu ?2> • • • > Qk, and also the time 

dxi dxi dxi dxi 



from which it also follows (by partial differentiation as to 
Q'l^ 92, • • • , qk) that 

dxi dxi dyi dyi dzzn dzzn 



dqi dq^ ' dqi dqi ' dg*' dqk 

9 



. (9) 
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Diflferentiating (7) with respect to any q we have 
dT ' / ,dx' , ,dy' , ,dz'\ 

while the derivative of (7) with respect to any q' will be (in 
view of (9)) 

dT ^ / ,dx , ,dy , ,dz\ 
This can be diflferentiated with respect to t [remembering that 



dt\dq J dq 



for any x; which can be easily seen by diflferentiating say (8) 
with respect to any q 

dxi _ d^xi d^j, 
'dq^didq'^'d^^^ '"• 

Also diflferentiating dxi/dq (from (2)) with respect to t: 



d dxi _ d^Xi d^Xi , 
dt dq dtdq dq 



+ -^q' + 



and comparing both results]; so that (for any q) 



d dT ^ / ,dx' , ,dy' , ,dz' 



„dx „dy „dz\ 



dq 

or, subtracting (10) and in view of (6) and (3) 

Assuming that all qi, g2, \", are independent, the only 
manner in which (12) can be satisfied is by equating the 
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coefficients of each 5q. Whence 

d^ dT dT 
dt dqi dqi 



T4 ^7~/ — IT = Qb 



dt dq^' dq2 ~ ^'' (13) 



dt dqk{ dqk " ^** 



There will be k such equations^ one for each independent 
coordinate. The generalized forces Qu Q2J " ', etc., are given 
by the k equations (5). Any desired generalized force, say Qm 
can be calculated from the general expression (see (5)) 

oqm oqm oqm 

by assuming that all other displacements, 6q, corresponding 
to other forces Q are temporarily made equal to 0, so that 
all the work done by the forces Xi, Yi, Zi, • • •, Xn, Yn, Zn, 
is represented exclusively by QmSqm- 

All we have to do is to find all the necessary derivatives 
(dxi/dqm), • • •, (dzn/dqm), from (2) and to substitute them in 
(14). We shall return to this in solving examples. However, 
it often happens that the given forces derive from some force 
function U; this may be the potential function as we had 
before, or, again, it may contain time (which the potential 
function does not). At any rate such a function possesses 
the property of having its partial derivatives, with respect to 
any coordinate, equal to the corresponding forces; so that for 
instance 

Y dU ^ dU „ dU 
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in that case instead of (5) or, say, (14) we have 

^" dxidqm'^dyi dqm'^ ""^dZn dqj 

which is simply = dUjdqm] so that, instead of the equation 
(13) we shall have, in this special case, 

d_ dT dT dU 
dt dqi dqi dqi ' 
(15) 



d dT dT ^dU 

dt dqk dqk "" dqk ' 

Both sets (14) and (15) will be simultaneous diflferential 
equations of the second order, and their integration will 
result in furnishing k independent coordinates or parameters 
in terms of time t There will be 2k arbitrary constants. 
The general order of procedure will be precisely the same as 
before: (1) Selecting k independent coordinates, such, that 
the coordinates x, y, z, etc., can be expressed in terms of them, 
in connection with the constraints. (2) Formation of the 
expression of kinetic energy and replacing in it the old vari- 
ables by their new equivalents qi,."', etc., qi, •••, etc. 
(3) Formation of U, if any; if not, computing the generalized 
forces Q; (4) formation of the equations (14) or (15) proper; 
and (5) integration of the latter. For the reason outlined 
above the integral of kinetic energy, if it exists, can replace 
any (generally the most complicated) Lagrange's equation. 

A few examples will now be given; the reader is advised 
not to skip any of them as " uninteresting " or " useless." 
They have been selected with great care and are being pre- 
sented in a certain order, which is designed to make matters 
absolutely clear. 

Example 1. A bug, of mass m, can crawl with a given 
velocity upon a wire (fig. 34) of the length 2a and mass m, the 
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ends of the wire sliding upon a circumference (in a horizontal 
plane) of radium R. Form the equations of motion and calculate 
the angle through which the wire will have moved when the bug 
reaches the end B, 

B 




Fig. 34. 

The notations will be as follows: the angle COx will be <p, 
and the distance mC = r. The bug starts from C and its 
velocity, v = r/t, is constant; also let Om = p, and mOx = a 
(polar coordinates of m). Here we have only one parameter 
or independent coordinate, the angle <p, to be found in terms 
of t, because the position of the bug upon the wire is always 
known for any time, t, through its velocity, so that r = vt. 
The kinetic energy of the system will be the sum of kinetic 
energy of the bug (mtr^/2), moving forward, and of the wire, 
rotating backward, about its instantaneous center 0[= /(w^/2), 
where / is the moment of inertia of the wire about and is 
equal to mP, k being the proper radius of gyration; of course 
CO will be simply = d<p/dt], so that the kinetic energy of the wire 
will be = (mA:V2)^'^• that of the bug will be (m/2)(p'^ 
+ p^a'^) (from the ordinary expression of velocity in polar 
coordinates, p and a). 

Hence 2T = ml{?(p'^ + m(p^ ^+ P^o^'^)- W e shall now elimi- 
nate the angle a; since p = ViZ^ — ^^ + '^^, it will readily be 
seen that p' = vHjp] also 

, , v'<B?- a" 

(from a = ^ + tan~^ ^[m~^ — 2 5 ^^^ triangle m OC I 
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Finally 



,2 



2r= mJfc^^'^H 2"+^ 



mi^ . / . . rV/P-o*V 






where p is a function of time only. There being no external 
forces, we shall have as our only Lagrange's equation, for the 
parameter (p, 



that is 






^ [mF^' + m(pV + ^ VF^^2)] = 0, 



whence iV + pV' + ^ ViP — o^ = const. = c. 

From the initial conditions (starting from rest) it follows 
that c = 0. From the equation just obtained we can find 






or, integrating, 



^ = ^0 — \ 



iP-a2 ^ _, tJ^ 



-tan 



-1 



where ^o is the initial value of ip. When the bug reaches B 
we have t>^ = a, which in the above equation gives the final 
value of ip. 

J Example 2. A particle of mass m (fig. 35) is placed inside 
of a thin hollow ring of mass M, which can roU, toithoiU sliding, 
and in a vertical plane, upon a horizontal line; no friction; no 
initial velocities. Find the motion of this system. 

Let be the initial point of contact, so that originally A 
was at 0; in the time t the contact is at P, so that the distance 
PO = the arc AP = rd, where 6 is the angle corresponding 
to the arc; also let <p be the angle, characterizing the present 
position of the particle m (which originally was higher, say at 
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mi). The only independent coordinates we need are the angles 
6 and <p, since they completely identify the system (as soon 




as they are found in terms of t). The coordinates of m (re- 
ferred to the old axes, x, y) are, 

X = r{d + sin <p)\ y = r{l — cos <p); 

the kinetic energy of the particle, in its own motion inside 
the ring, is 

"2" = "2 (^ + 2/ ) = 2" (^ + ^ + 2«V cos ^) 

(note that here again absolute coordinates have been used, 
although the rules of relative motion would have resulted in 
precisely the same expression; compare with the end of 
Chap. II). 

The kinetic energy of the ring proper will be based upon its 
instantaneous rotation about the point of contact, P (but not 
about C, which itself is in motion) ; the instantaneous value of 
the angular velocity about P is dd/dt or d\ The moment of 
inertia about P will be 2Mr^ (compare with Bowser, Anal. 
Mech., example 9, p. 448, where a = 6) ; hence the total kinetic 
energy 

2T= mr^{e'^ + ^'' + 29V' cos ip) H-^Mr^fl''. 
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The only acting force being gravity, we have a potential 
function, [7 = — mgy = — mgr{l — cos <p) (— , because it 
acts against positive y). 

We now have everything necessary for forming Lagrange's 
equations. Taking partial derivatives of T with respect to 
and <p; also to their derivatives fl' and <p\ we have 

^ = mr2(r + <p' cos <p) + 2MT^e'; fj = 0; 

-r—f = mr{<p' + S' cos <p)\ T~ = — w^^V sin <p 

0(p 0(p 



also 



dU ^ dU 

d9=^' a^=-^FSin^; 



hence Lagrange's equations are 



/ 



j^ [mid' + cp' cos <p) + 2Me'] = 0, — (1) 



mr 



(dip\ de'\ , 

['jf+ ^^^^li )^ mg sm<p = 0. (2) 



Integrating (1), we have m(9' + <p' cos <p) + 2M0' = (the 
constant of integration is = because of the initial conditions: 
starting from rest means ^o' = 0; 9o'.= 0). Integrating 
again, 

m{d + sin ^ + 2M0) = m sin a 

(where a is the initial value of (p when the particle was at rest 
at mi; whence 

TYl 

— cnjT t (sin a — sin ^). (3) 

2M + m ^ ^^ ^ ^ 

On the other hand the equation (2) gives 

dip' dO' g . 

■^+cos^^=--sm^; 
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substituting dO'jdt (derived from (3)), multiplying by dipjdt 
and integrating, we have 



/2 



m 



2<7 



o TIT I " ^ cos^ ^ = — (cos ip — COS a) 



(check back by differentiating this with respect to t) ; whence 



d<p _ hg{2M + ^) / cos <p — cos a 

dT "" ^ \ r \23f + m sin^ ^ 



(4) 



This angular velocity, which was evidently equal to when <p 
was = q: (in the beginning of motion) will again be equal to 0, 
when <p = — a; that is the particle will reach the same height 
on the other side of the contact point, at m2, and will oscillate 
between these two points. Meanwhile [see (3)] the angle 
will have reached a certain maximum (when ^ = — a) 
= (2m sina)/(2Jf + m), and will be back to when (p is 
at mi again; of course y oscillates between r(l — cosa) 
and 0; while x will vary from 

m-2M . 



r sm a 



to 



m+2M 



r sm a. 



Example 3. Two rods, A and B (fig. 36), connected at A 




Fig. 36. 



hy a flexible joint, are placed upon a frictionless horizontal plane. 
The length of each rod is = 21, and the mass of each = M. 
Find the general equations of motion. 
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In selecting the independent coordinates for this case we 
are confronted with the necessity of having not less than 
four such parameters: two coordinates, a and b, of the center 
of gravity G of the system, which is half-way between the 
middle points of each rod, C and Cii one angular parameter (p, 
the angle of the line G-A with the axis of x; and another 
angular parameter a, = half the angle between the rods. 
We can readily satisfy ourselves as to the possibility of recon- 
structing the system from these four parameters: two first 
coordinates locate the center, G; then the line G-A is drawn 
and the length I cos a is laid off upon it; this gives the location 
of the joint A; finally the rods are laid off at the angle a on 
both sides of G-A. 

The kinetic energy will consist of two parts (Koenig*s 
theorem): (1) kinetic energy due to the whole mass 2M, 
considered as concentrated at the center of gravity, G, which 
equals 

and (2) kinetic energy due to the motion about the center of 
gravity; this kinetic energy can likewise be split into two 
separate items, for each rod: the part due to the motion of C, 
the center of gravity of the rod, with the whole mass M con- 
centrated in it = ilf^/2; or, in polar coordinates, r, d (where 
fl is = XiGC = ^ + (7r/2)), this is 



-m)'-^<f)'] 



or, substituting r = GC = I sin a, we have 

-y (Pa'^cos^ a + Pip'Sin^ a). 

The other part is due to the instantaneous rotation of AB 
about its own center of gravity, C, the mass of the rod being 
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M (and the radius of gyration = k) and the angular velocity 
being = {d\l//dt), where yp is the changing angle X2CA, 
which = ^ — a. This part of the kinetic energy will then be 

The corresponding expression for the other rod, ABi will be 
found by substituting — a for a; so that the total kinetic 
energy will be 

2T = MW^ + V^ + (P cos2 a + ie)a'^ + (P sin^ a + ¥W\ 

The mechanism of finding this kinetic energy is extremely 
instructive and should be thoroughly mastered by the reader). 
There are no external forces, therefore U' = and all Q 
will be = 0. 
Writing down Lagrange's equations for a we have 

da' 

-^ = 0, or a' = const.; 

in the same manner V = const., so that the motion of the 
center of gravity is rectilinear and uniform. 
The equation for <p will be 

dt dip' dip ' 
whence (since T does not contain ip) 

d dT ^ dT 

^d^=^' ^" d?=^^^'^- ^" (1) 

(P sin^ a + l^W = C. 

Lagrange's equation for a would be too complicated and 
we shall replace it by the integral of kinetic energy T = [7 + A, 
or, in our case, T = const.; that is 

(P cos2 a + B)a'^ + (P sin^ a + ¥)ip'^ = A\ (2) 



^ 
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where -4 is a constant, absorbing also a'^ and 6'^, both of 
which are constants. From (1) we see that ^' has a constant 
sign, so that the line G-A is rotating in the same direction 
about G, the angular velocity having the limits C/k? and 
C/iP + li?). Substituting fl' from (1) in (2) we have 

a''(F cos2 a + A:2)(p ^^^2 a+Jc') = A'P sin^ a + A^Jc" - (?, 

of which since the first member is positive, the second member 
must be positive as well. Three assumptions can be made: 
(a) If C^ — A^k*^ is negative, the value of a may be any, and 
the angle between the rods will (according as a' is + or — ), 
increase or decrease until they come in contact (that is a = 0, 
or a = tt) ; (6) if C^ — A^lc^ is positive, we can always select 
a constant /3, such that C^ — A^lc^ = A^P sin^ /3; indeed 
C^ — A^](^ is always less than A^P, since at the start a = ao, 
and regardless of the value of ao we always have A^P sin ao 
> C^ — A^lc^. Therefore, in this case, the only condition for 
a will be sin^ a > sin^ /3; so that a can take any value between 
j8 and tt — /3, and the rods will oscillate in relation to G-A; 
(c) if C^ — A^h^ = 0, a can take any value, but the less the 
value sin a, the slower the angular velocity a'; the rods will 
therefore tend to close either way (a = 0, or a = tt) without 
ever reaching that state. 

This discussion is due to M. Appell and is especially in- 
structive. 

Example 4. A system (fig. 37) consists of two rods A-B 
and C-D, each of mass p and of length 2a; also of two rods, 
A-D and B-C, each of mass q and of length 2b; the flexible 
joints at the corners are frictionless. Find the motion of this 
system. 

Here we need but two independent coordinates, the angles 
6 and <p, of the rods A-B and A-D with the vertical. The 
kinetic energy of A-B will be due simply to its rotation about 
0; the moment of inertia being = (mP/12) (Bowser, Anal. 
Mech., p. 431); we have 
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To find that of D-C we must add to the above (which is the 
same as the kinetic energy of D-C about its own center), the 
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Fig. 37. 

kinetic energy due to the motion of its center about 0; and 
that is = (p/2)^Z)2.^/2 = 2pbY^; so that 

K. E.^^ = kpa'd^' + 2pbV'; 

In order to find the kinetic energy of A-D and B-C we must 
add the double kinetic energy KfgtV^)^ due to their instan- 
taneous rotation <p' about their centers of gravity E and F, 
to the double kinetic energy due to the mass q of each, con- 
centrated at E and F, in the motion of these points, that is 

h[2q{bV' + a'e^')]. 

(Remark: This is quite evident; according to Koenig's the- 
orem we first calculate the kinetic energy of F in its motion 
about A, which is qh^<p'^; and then add the kinetic energy 
due to q, concentrated in A, in its motion about 0.) Finally 
the total kinetic energy is 

2T = 2p(26V'' + WO'^) + 2qiaW^ + |6V"). 

We have a potential function U, equal to the mass times 
the distance of the center of gravity of the whole system 
under 0; the mass being 2p + 2q we have, 

U = 2(p+ q)gb cos <p. 

Forming, now, the two Lagrange's equations in 6 and (p we 
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have 

dT _ 2p+6g ar _ 4(3p+2g) , 

de'~ 3 ^^' d<p'~ 3 *^ ' 

^ = 0; 3^ "" ~ ^^ + ?)?& sin ^. 



Hence the equations 
d f2p+Qq 



a^' ) = 0, 



di\ 

^ (|(3p + 2qWip') = - 2(p + g)sf6 sin v^. 

From the first equation we can conclude that V is constant^ 
so that the rods A-B and C-I) turn with a constant angular 
velocity. From the second equation we have 

<P^_ 3 P+g g. . 
(ft2 - 2 3p + 2g 6^^^^* 

comparing this with Bowser, Anal. Mech., p. 458, we conclude 
that the motion will be pendular, equivalent to a simple 
pendulum of length 

6p + 4g 



3p+3g 



6. 



Remark: We already had equations of this sort on two 
occasions; it is well to commit to memory the fact that 

^^ g . 
■j^ = sm ^ 

means pendular motion, equivalent to a simple pendulum of 
length a. 

A few examples will now be given on the motion of rigid 
bodies; their position is generally given by means of Euler's 
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angles^ locating the principal. axes; the reader should^ there- 
fore, revise the end of the first chapter. 

Example 5. A body of revolution is mounted on a straight 
wire R-S (fig. 38), coinciding vrith its aods, one end of which 




Fig. 38. 

can slide along a vertical 0-Z, while the lower end, S, can slide 
in any way on the horizontal plane x-O-y; no friction and no 
external forces except gravity. Find the motion. 

Three angular parameters are necessary to define the posi- 
tion of the body: two will locate the wire and the third param- 
eter will define the body in its possible rotation around (and 
including) the wire. 

The body is not shown on the sketch; its center of gravity 
is at G; through 6? draw axes x, y, z, parallel to X, Y, Z, and 
let N-G be the intersection of a plane, perpendicular to the 
wire R-S with the plane x, G, y. Call xp the angle xGN; 
6 the angle RGz\ and <p the angle between GN and any line 
GP in the plane perpendicular to R-8 and fixed in the body. 
Then, the angle <p will completely locate the body about R-S] 
also 0R8 = 6, thus fixing the position of the wire about 0-z. 
On the other hand, the angle NG8 which is a right angle, 
projects as a right angle upon the plane x-O-y (since one of 
its sides, NG, is parallel to that plane); hence XOS = xp 
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— (7r/2), which fixes the plane ORS in relation to XOZ. 
These three angles ip, and (p, are Eider's angles and the 
position of the system is clearly identified by them. If given 
the angles ip, and <p we could proceed to locate the system 
as follows: From the angle xp we can find xp — (7r/2), and 
thereby locate the plane ROS; the angle will definitely 
locate the wire in that plane; finally, the position of the body, 
in whatever rotation it might have about R-S, will be identi- 
fied by any such line as PG (belonging to the body), through 
the angle <p, from the intersection NG. The motion which 
we are investigating is perfectly general; the end R may move 
down along Z; the end S may, at the same time, describe 
some arc S-Su and, meanwhile, the body itself may turn to- 
gether with its axis about R-S, Whether all these motions 
will or will not take place, in this most general case, will be 
seen from the equations of motion; at any rate the motion 
is perfectly defined by these Euler's angles; they will be our 
independent coordinates. 

The next step is to find the kinetic energy of such a system: 
this will (Koenig's theorem) consist of two parts: (1) The 
kinetic energy due to the motion of the center of gravity, G, 
with the whole mass, M, concentrated in it; also, (2) the kinetic 
energy due to the motion of the body about the center of 
gravity, considered as fixed. Let us calculate these items 
separately: 

1. Let the distance RG = r, and RS = I; the coordinates 
X, y, z, of G, are as is very clearly shown on the sketch 

a: = r sin fl sin ^; y = r sin cos ^; z = (l — r) cos (p. 

Differentiating with respect to t and substituting into 

2T= if (a:'' + / + z''), 
we have 

2T = if[r2(sin2 fl.^'^ + cos^ 0-0'"^) + (l - r^ sin^ 0-0% 
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2. The kinetic energy of the body about its center of gravity 
can be due only to its (possible) rotation about 6?, considered 
as a fixed point. We had a general expression of kinetic 
energy for such a case (see under Body with one point fixed) 

where A, B and C are the moments of inertia about the 
principal axes of the body and coi, C02, C03, are the components, 
of the instantaneous velocity of rotation, on some initial axes, 
in relation to which the principal axes are given by Euler's 
angles. We do not care to have the instantaneous velocity 
itself, but we had the following expressions for coi, C02, C03 in 
terms of Euler's angles 

coi = sin 6 sin ^'^ + cos (p'B', 
C02 = sin cos ^•^' — sin <p'6'y 
(jhz = cos 0'\l/' + <p' 

(see under Euler's angles, at the end of chapter 1). In bodies 
of revolution it is customary to denote the two equal moments 
of inertia by A and B\ in our case A = B are the moments of 
inertia about any two principal axes at right angles to each 
other, in the plane, perpendicular to iJ — iS; while the moment 
of inertia about R — S will be = C. Therefore in our case 

2T= Aia)i' + W2^) + Cws^; 

which, after substitution of wi, u)2, coz, just given, becomes 

2T= A{sin^ S-^'" + O + C(cos fl-^' + (pj 
so that the total kinetic energy of the moving body is 

2T= iA + Mr") sin2 fl.^'' +[A + M{r^ cos^ 

+ (I- r)2 sin2 eW^ + C(cos fli • ^' + <p'y. 

The only external force is gravity, deriving from the poten- 
tial function U = Mgz, that is = Mg(l — r) cos 
10 
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In forming Lagrange's equations we can see at once that 
both T and U contain no ^r or ^ ; thCTefore we have 



that is 



also 



7 = const.; or c' -:- ^' cos 9 = const- = a; (1) 



tt;= const.; or 



f^ + Mr") sin^ 0-^ + C cos diip' + ^' cos ») = const.; 

(2) 
or ^' sin- d + ^w cos ^ = 6 



[ where m=^j^j,^] 



Instead of the third Lagrange's equation (in 0) we can 
substitute the integral of kinetic energy, T = C + A, that is, 

i/^ sin* + (1 + n sin* e)^'* = p cos » + A (3) 

(where 

_ MI Q, - 2r) _ 2Jf^(Z - r) 

^■" ^ + Jfr2 ^"""^ ^" ^ + Jifr2 ' 

eliminating ^ between (2) and (3) we shall have i = JF{B)dB, 
where F(S) is a function of 6 



T.//>v / W- w sin* S V . . /» 

\ (p cos d + A) sm* 6 — (6 — m cos d)* 

therefore, from (2) we have 

b — m cos 6 



*-i 



sin^'fl 



F(fl)<W 
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and finally 

SO that the problem is reduced to integration. 

This example has been given here as an illustration of the 
mechanism of the introduction of Euler's angles; for this 
reason we shall omit the interesting discussion given in M. 
Painlev6's lectures, from which the example has been taken. 

Example 6. A solid of revolution, having one point of its 
axis fixed at 0, and whose initial rotation about its axis is given, 
is subject only to the external force of gravity. Investigate the 
motion of such a top or gyroscope. 

The position of a rigid body of this sort can be given at 
any time through Euler's angles yj/, 6, <p, which we shall take 
for independent coordinates, the axis of spin being z. Let 
(fig. 39) G be the center of gravity and let OG = I; the poten- 




tial function will be 17 = Mgl cos 6 and the expression of 
kinetic energy, 2T = Ao)i^ + Bo)2^ + Cca^, will be given by 
the formulae (see example 5) 

coi = sin sin ^•^' + cos <P'0'\ 
C02 = sin 6 cos (p-^p' — sin ip-B'] 
0)3 = cos S • ^' + ^' 
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from which, considering that here, too, we have equal mo- 
ments of inertia, ^ = J?, we have 

2r = A{e'^ + ^'' sin2 6) + C{<p' + ^' cos dy. 

The three Lagrange's equations, in <p, yj/ and are 

^it ^^' + ^' ^^^ ^^ ^ ^' 
^ (4^' sin2 e + C{ip' + ^' cos 6) cos «) = 0, 

Jar 

A-^- Ail^'^ sin fl cos e + C(^' + ^' cos fl)^' sin S = 

— Mgl sin ff; 
from the first two we have 

^' + ^' cos fl = n; (1) 

A\l/' sin2 e + Cn cos e = ifc; (2) 

n and A; being constants; the third equation gives 

JQf 

A-^- AiP'^ sin e cos e + Cn^' sin » = - M ^Z sin ». (3) 

Eliminating ^' from the last two equations we can obtain 
a rather complicated expression of 6 in terms of t. 

But we can investigate these results in an elementary 
manner, and derive therefrom a few fundamental principles 
of gyroscopic motion. 

We shall begin by replacing (3) by its equivalent, the inte- 
gral of kinetic energy T = 17 + A; in our case 

2T= Aid'"' + ^'' sin d) + C{<p' + ^' cos Bf 

and 

U = Mgl cos 6 (see also (1)), 

and therefore the integral of kinetic energy will be 

A(e'^ + ^'' sin2 0) + Cn" - 2mgl cos = 2 A, (4) 



"^ 
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from which ^' can be eliminated through (2) ; so that 

A sm^ 6 ^ 

Now in the equation (2), let a be the initial value of 0, that 
is let S = a when^< = and when ^' = 0; from this sub- 
stitution we have 

k = Cn cos a and Aip' sin^ 6 = Cn(cos a — cos 6). (6) 

Also k substituted in (5) results in 

./>,2 . ^ w^(cos a — cos S)^ . XV o ^ 7 /» r^T "^"^ '* 

Ad'^ + C^ — — -.--Y-n + Cn^ - 2m5fZ cos = 2h; 

A sm^ S ^ 

so that, in the beginning of motion, when 6 = a and fl' = 0, 
we have 

Cn^ - 2m5fZ cos 6 = 2A. (7) 

This substituted in (5) gives 

(cos 6 — cos a)[2mglA{l — cos^ 6) + C^n^(cos fl — cos a)] = 0. 

This is a cubic equation in cos 0, one of the roots being, of 
course, cos 6 = cos a; and the other two roots can be ob- 
ained from the quadratic 

(cos^ fl - 1) - 2rflgiA ^^^^ ^ "" ^^^ ^^ ^ ^' 
here we can put 

2mglA " ^' 
so that 

(cos^ fl — 1) — 2g(cos 6 — cos a) = 0; 
hence 

cos B = e± Vi + 6^— 2g cos a. 

The greater root will always be > 1 (cos a being always < 1 ; 
if cos a = 1, then only can cos fl be = 1), and therefore should 
be rejected, since of course cos 6 should be < 1. 
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Therefore the only solution is 

cos di= e — Vl + 6^ — 2e cos a, (8) 

This shows that the value of 6 (and therefore the axis of the 
top) will oscillate between 6 = a and 6 = 0i, just found. 

If the rotation of spin is high, we can develop the radical 
and neglect powers of higher than second, so that the simpli- 
fied expression will be 

cos = cos a sin^ a; 

e 

that is to say, the limits of 6 are practically a and 



(sin^ a \ 
cosa--y-y 



cos M COS a 1 . (9) 

Precession. From (7) and (4) we can easily see that 

Ad'^ + A\f/'^ sin2 e + 2mglicos a - cos 6) = 0, (10) 

where, either cos a — cos S = 0, or #= 0. 

1. If cos a — cos 6 is not = 0, then xp' has a maximum 
whenever S' = in that case 

,2 _ 2mgl cos a — cos 
^ "" "^ ^^ ' 

or, substituting 

A sin^ Cn 

cos a — cos ^' 
from (6) we have 

, _ 2mgl 



^'' = 



Cn ' 



this is the value of the angular velocity of the line Xi about 
0-Z (see Euler's angles on the sketch), or, in other words, 
the angular velocity with which the line 0-z describes a cone 
about 0-Z. This is not a constant but an instantaneous value, 
at the moment when the angle is in its furthest position 



System. 139 

from its initial value a and is reversing {0' = 0) the sense of 
its oscillation. 

2. When the axis reaches its initial value, we have cos a 
= cos By and therefore we can conclude (from 10) that at 
that instant hot h 0' and ^' = 0; in other words at that moment 
the axis of the top is temporarily stopped in its double motion 
{away from and about the axis 0-Z). Therefore the angular 
(azimuthal) velocity of the axis of the top about the vertical 
changes from to 2mgl/Cn, while, at the same time, the 
inclination of the axis to the vertical changes from the initial 
value a to its other value, di, which we had derived above (8). 
Therefore the motion of the top on its axis will project upon 
the horizontal plane in a manner shown in fig. 40, where the 




Fig. 40. 

inner circle corresponds to the initial value 6 = a, while the 
outer circle represents the extreme value of 6 (= 6i). 

The azimuthal motion about the vertical is called precession; 
we shall presently show that its average value is mglfCn; 
this is rather slow, since n is always considerable. The 
difference between the angles di and a is almost imperceptible 
(in view of (9) the angle varies between cos a and cos a 
— {2mglA sin^ a/C^n^), where n is generally very great) and 
the oscillation of the axis of the top between these values (so- 
called nutation) is of very short period and can but rarely be 
observed. 

Steady motion. If we assume that the variation of the 
angle d (nutation) is so slight as to be negligible, in other 
words that = const., and, consequently, that S' = 0, we 
have, from (3), dropping dd'/dt and neglecting ^'^, (since xf/' 
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is slow), 

which is half the maximum value, derived above. 

This equation can be arranged in a slightly different form, 
in which it is extensively used in gyroscopic work: multi- 
plying by sin 6 we have mgl sin 6 = Cn\l/^ sin 0; the left mem- 
ber is the external moment, M, whose arm is I; the right 
member is ^' projected upon a direction at right angles to 
the axis of spin of the gyroscope, and to the axis of the ex- 
ternal moment (fig. 41). Using notations customary in 




Fig. 41. 

gyroscopic practice we have M = IQcj, where J is = C, that 
is, the moment of inertia of the gyroscope; Q is the velocity of 
spin, and co, that of the precession about an axis at right 
angles to the axes of both the spin and of the external moment. 
This formula enables us to calculate the moment necessary 
to create a certain precession; also, the moment that will 
result from giving the so-called forced precession to a spinning 
body. It can also be put in the following form 

WR'Nn 
" 2937 ' 

where M is the resulting moment in ft.-lbs.; W the weight of 
the spinning body of lbs.; R is the radius of gyration in ft. 
and N and n are the velocities, of spin and of precession, both 
in revolutions per minute. 

Remark. The reader will now understand the significance 
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of the constant n (see equation (1)). It is not the velocity of 
spin; the latter is = <p\ and forms only a part of n; the other 
part, ^' cos 6, is the component (upon the axis of z) of the 
precession, ^'; on the other hand xp' is never very great and 
therefore the value of n does not differ materially from the 
velocity of spin. 

Example 7. To derive Euler's equations from Lagrange's 
equations. 

The expression of kinetic energy for a body with one point 
fixed is 

2T= Aoii" + Bw2^ + Cuiz\ (1) 

where coi, co2, C03 are the components of the instantaneous 
rotation, co, upon the principal axes (strictly speaking upon 
fixed directions, temporarily coinciding with the latter). 
These components can be given in terms of Euler's angles 
by the transformation formulae which we had before 

0)1 = ^' sin S sin ^ + 0' cos <p, 

0)2 = ^' sin cos (p — 0' sin ^, (2) 

0)3 = ^' cos + <p\ 

where Euler's angles, ^, 0, ip, identify the principal axes of 
the system, taken about the fixed point. Taking the angles 
^, Bj (p as independent coordinates and putting virtual work 
hW = ^8^ + 98^ + ^8^ (where ^, and $ are generalized 
forces, in our case moments) we can write down one of La- 
grange's equations, in (p, as follows 



d^fdT 






dt \d<p 

Now, in order to derive d T/dcp', we might substitute (2) into 

(1) and differentiate the result with respect to ^'; but this is 

not necessary, since, from (1) and the last equation (2) we 

have 

dT dT dcoz 



d<pi dcjz d<pi 



= Co)3 
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and also 

dT _dT du)i dT du)2 _^ j 3coi do)2 

d<p dwi d(p do)2 d<p d(p dip 

Differentiating coi and 0)2 (transformation formulae), we have 

doix , , du)2 

^ — = — 6' s\n<p + ^' cos <p sin = C02, x — = — coi, 

Therefore 

dT 



we have 



X— = 0)10)2 M ~ B) 



That f>, the generalized force corresponding to the gen- 
eralized displacement d<(>, is merely the external moment, N, 
about the axis z, can be shown as follows: Oiu: rule requires 
forming of the equation 

dz\ 
dq) 



«=K4:+4:+^^-^^ 



for every generalized force Q we wish to determine. In this 
particular case Q is $ which we want to find; and, considering 
that we are dealing with the rotation <p, that is about the 
axis z, assuming other displacements constant, we have (for 
rotation about z — ^polar coordinates x = r cos ^; y = r sin ^) 

dx . dy ^2 

-=-rsm^=-2/; ^=rcos^ = x; ^=0; 

and therefore 

which is evidently the moment N about the axis z; therefore 
^ = N, and we have the third Euler's equation 

C^+ (B-^)o)io)2 = N. 
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This equation contains no angles ^, d and (p, and therefore is 
perfectly general in regard to all axes, subject, of course, to 
suitable change in the notations. 

Note, Holonomous and not-holonomous systems. In de- 
riving Lagrange's equations mention has been made of the 
perfect freedom with which generalized coordinates can be 
selected. The only implied condition was that the ordinary 
position-coordinates, x, y, z, etc., can be expressed through 
the generalized coordinates by means of some such equations 
as (2) p. 89, which may or may not contain time. A very 
important remark should be made, however, regarding the 
choice of these independent coordinates. The equations (2), 
just mentioned, should contain no time derivatives; otherwise 
Lagrange's equations cannot be applied. In practical work 
the reader will but seldom meet with difficulties of this nature, 
and in more advanced treatises he will find how corrections 
should be made so that modified Lagrange's equations can 
be applied to such special problems. Systems which can be 
characterized by equations like (2) involving no time deriva- 
tives, are known as holonomous; those, which in the equations 
like (2) would contain time derivatives are called not-holo- 
nomous systems. For instance, in evolving Euler's equations 
by Lagrange's method (example 7, p. 141) what was our 
object in introducing the transformation formulae (2) into 
the expression of kinetic energy (1)? Why could not La- 
grange's equations be obtained directly by selecting coi, C02, C03 
as independent coordinates, and by differentiating (1) with 
respect to them? The reason lies in the fact that the co- 
ordinates of any point of a body can be expressed through 
coi, C02, C03 only by means of relations involving derivatives 
dz/dt, etc. (see p. 48). Therefore we had recourse to the 
transformation formulae (2) taking Euler's angles as inde- 
pendent coordinates, well knowing that by means of Euler's 
angles the system can be completely characterized by expres- 
sions involving no derivatives (see p. 80). 



CHAPTER IV. 
Lagrange's Equations for Relative Motion. 

We have seen why the relative motion of a particle does not 
need any new methods. In dealing with a system, three 
separate methods can be used for forming Lagrange's equa- 
tions of motion. 

First method. Let it be required to find the motion of a 
system in relation to moving axes x, y, z, the motion of which, 
in reference to some fixed axes, X, Y, Z, is known. We can 
select independent coordinates, qi, q2, • • • , etc., in relation 
to the moving system x, y, z, and then deal with the problem 
as if the motion were absolute, since any parameters, say 
Qi, q2, ' ", characterizing the system in reference to moving 
axes, X, y, z, will also define it in reference to the fixed set, 
X, F, Z, the motion of x, y, z, themselves being known. In 
calculating the energy, however, the following should be 
remembered: The expression of kinetic energy is 

2T= SmCa:'' + / + A 

where x, y, z, are absolute coordinates, which should be found 
in terms of qi, q2, • • • , and then substituted, which will give 
2r in terms of q's and g"s. A much simpler way would be 
to form the expression of energy by means of the formula of 
absolute velocity in the most general case of motion (see (2) 
under Relative motion) 

dx 

Vax = ^ + '^rnarx + C022 — COsT/, 



d^ 
dt 






T7 ^^ . I 

"^ di ^'""" ^^^ ~ ^^^' 
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where Vax, etc., are the projections of the absolute velocity; 
dx/dt, etc., are the projections of the relative velocity on the 
relative axes, x, y, z (parallel to X, Y, Z); Vmarx, etc., are the 
projections upon X, Y, Z, of the (mar) motion of the moving 
axes of reference; while coi, C02, C03 are the instantaneous rota- 
tions about the axes z, y, z. 
Hence the final expression of kinetic energy 

2r = Sm ( ^ + 'Omarx + CO2Z — COs^ j 

(dy \2 

^ + Vmary + 0)zX — COiZ I 

+ ( ^ + Vmarz + COiT/ — (a2XJ . 

This expression though apparently long is easily applied in 
practical work. Its advantage is that it involves only rela- 
tive coordinates, which are easily found in terms of q, g', etc., 
while all the velocities dx/dt, • • •, etc., as well as Vmarx, • • •> 
are known. In calculating the generalized forces, Q, the same 
methods will be used as those given before; and if there is a 
force function U, then of course Q = dU/dq for each q. 

Example. The plane x-O-z can turn, about a vertical axis 
0-z, with a constant angular velocity, co. Find the motion of a 




X 



Fig. 42. 



rod A-B, constrained (fig. 42), to move in such a plane, and 

given any initial motion in relation to the axes x-z; no friction. 

The motion of the rod in relation to the moving axes x-O^z 

will be defined by three independent coordinates: a and b. 



146' Lagrange's Equations. 

the coordinates of its center of gravity, G; and <p, the angle 
of the rod with, say, the axis x. The absolute velocity of 
any particle, m, is here made up of the relative velocity, which 
is in the plane x-O-z; and of the {mar) velocity, which is = o)X, 
and is perpendicular to the moving plane at that point, m. 
Hence the absolute velocity Va^ = Vr^ + {vmarY so that the 
kinetic energy 2T = Xmv^ = llmvr^ + llmvmar^. These two 
terms will be found separately. The relative motion of the rod 
consists of translation and of rotation; according to Koenig's 
theorem we have Smur^ = M{a^^ + 6'^ + k^<p'^), where M is 
the whole mass and k the radius of gyration of A-B about G; 
while Xmvmar^ = co^'^moi?; but l^mx^ is the moment of inertia 
about the axis z, which as we know, equals the moment of 
inertia about zi plus mass times the^square of the distance be- 
tween the axes, l^mvmar^ = Mo)\a^ + k^ cos <p). Hence finally 

2r = M{a'^ + V^ + P^'' + co2jfc2 cos2 <p + coV). 

There exists a potential function of external forces, TJ = Mgb 
so that Lagrange's equations for a, b and <p will be 

|(a')-co^a=0; |(6') = (/; 

IT (P^O + Pco^ sin <p cos ^ = 0, 

whence a, b, <p result in terms of time, as a direct result of 
integrating these equations. 

From the first equation we have, a = Ae'^* + jBe"""' from 
the second, b = g{f/2) + Ct + D; and the third equation can 
be reduced to 

<p<p g , 

-TT H sm 2<p = 0. 

or m 

This means pendular motion, equivalent to a simple pendulum 
of the length 2m. 

It will thus be seen that the first method contemplates 
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finding the (absolute) kinetic energy and then expressing it 
in relative coordinates. 

Second method. According to this method we add, to the 
acting forces, the two additional forces (see under Relative 
motion, Chapter 1), the centrifugal force {Fmar), due to the 
motion of relative axes; and Coriolis's force, (Fc); having done 
so we can form the expression of kinetic energy in relation to 
the moving axes as if they were fixed; and, in forming the 
generalized forces, Q, etc., according to the general rules, we 
have to consider the two additional forces. An example will 
make this clear. 

Example. Let us apply this method to the problem, for 
which the solution has already been given in two different 
manners (figs. 31 and 33). 

A particle is sliding in a circular tube (fig. 43) which itself is 
revolving with a constant angular velocity about a vertical axis 0. 




Fig. 43. 

Let 0-x be the reference axis, from which the rotation 
starts, so that the angle (at is sufficient to define the circle by 
its diameter, 0-D) the circle is the moving system, and only 
one parameter, tp, is necessary to fully locate the point M in 
relation to this moving system. According to the principles 
of the second method, all we have to do is to add to the 
external forces two other forces; (1) that due to the {mar) 
motion of the revolving system itself; and (2) the CorioUs's 
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force, of which the expression has been given above. There 
are no external forces; the centripetal acceleration r?lr is 
directed radially; and so is Coriolis's acceleration, so that the 
only force to be considered is the force exerted on M by the 
rotation of the {mar), relative axes, that is of the circle itself, 
with the angular velocity co. This is clearly the centrifugal 
force, M-C, and it equals mco^6, where 6 = 2r cos (^/2). 
The kinetic energy is 

In calculating the only generalized force, f», we reason as 
follows: What force, f», multiplied by the corresponding dis- 
placement, 8^, will give the same work as that, done by the 
force M-C = moy^b, during any possible change of b, that is 86? 
Of course, b = 2r cos {<p/2) so that 

simp <p 
00 = — r ty cos 2 o<p; 

therefore M — C X Sb = — mr^co^ sin ^ • 6^ or 

^S(p = — mr^o)^ sin <pd<p; 

. so that 

* = — mr^cJ^ sin <p. 

The only Lagrange's equation we have to form is (for <p): 

dt dip' d<p 



hence 

^(P 

dt 



|- = — 0)2 sin <p, 



the same as we had before. 

The reader is invited to re-write this solution for a sphere 
of mass M and radius of gyration k (sliding on the revolving 
ring, instead of a particle). This will mean certain changes in 
the expression of T] apply Koenig's theorem. 




tk. 
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iydrd method. This method is due to M. Gilbert, and is 
espeiViAUy interesting. We know that, in order to determine 
the niotion of a system in relation to a set of reference axes, 
^9 Vy *> which itself is in motion in relation to fixed axes, 
X, y, :Z, we must apply at every point of the system two 
additioni%l forces: that is the {mar) motion, and Coriolis's 
force. iLQagine, however, that the motion of the reference 
axes is refi^rred not to fixed axes, X, F, Z, but to an auxiliary 
set of axes, Zo, Fo, Zq (fig. 44), which are always parallel to 




Fig. 44. 

X, F, Z. In other words, the motion of Xo, Fo, Zo is trans^ 
latory in regard to Z, F, Z; while the set x, y, z, can have any 
rotation about the point 0; so that these two motions, trans- 
lation and rotation, enable the system x, y, z, to assume any 
position whatever in regard to X, F, Z; then, all we need to 
do is to apply, at every particle of the moving ssystem, corre- 
sponding (mar) forces, after which the axes Xo, Fo, Zq, can 
be considered as fixed. 

Owing to the translatory nature of the motion of Xo, Fo, Zo, 
the (mar) forces will be the same for every particle of the 
system; if the (mar) acceleration is J, the corresponding 
additional force will be (— mJ). Let the projection of such 
forces upon the reference axes, x, y, z, be — mJx, — mJy, 
— TnJg. To every virtual displacement, 6a:, dy, dz, of any 
particle, there will be a corresponding element of work done 
11 



f 
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by these forces, so that for the whole system this will be / 

- i:m(J^dx + Jydy + J^dz), ' (1) 

where Jx, Jy, Jz, are known functions of t. «• 

If we form an expression 

K = — 2^m(xJx + yJy + zJg), 
which is also 

= — [Jx^mx + Jyl^my + J«S?»a] 

and therefore 

= - M{aJx + hJy + cjr,), /' (2) 

where a, 6, c, are the coordinates of the center .<)f gravity, G, 
in relation to the axes x, y, z (Bowser, Anal. Mech., p. 110), 
we can say that (1) is = dK; it also follows that 

K= - il/.J.OG.cos(J,0<?) 

(from the well-known formula of Anal. Greometry: aee Bowser, 
Anal. Geom., p. 268, formula 6). 

Having thus introduced the only correction necessary (due 
to (mar) forces) we can now consider tbe motion as absolute. 
Of course, in Lagrange's equations, o£ the type 

d ST^dT dK 

dt dq' dq ^~ ^ dg' 

there will appear the last term dK/dq, for every q, as corre- 
spondiag »>hare of generalized forces, for any generalized dis- 
placement, dq, due to the additional function K, which we 
have just introduced. If the actually applied (external) 
forces have a potential function, we simply have the second 
member = (d/dq)(U+ K). 

Regarding the kinetic energy of this motion the following 
remark can be made: We have seen (equation (2) under 
Relative motion), how the absolute velocity of a particle can 
be given through its motion in relation to the moving axes, 
as well as through the rotation of these axes themselves (we 
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can omit the terms (vmarx), {vmary), (vmarz) because our motion 

of the reference axes is about a fixed point, 0); therefore we 

have 

dx 
K ax = ^ + a?22 — cosy, 

dy 

^^al, = ^ + 0)3X — COiZ, (3) 

' ^ di ^^^ ~ ^^^' 

Substituting these expressions into the usual expression of 
kinetic energy 

we have 

2r = Sm[(a:' + (^2Z — casyY + (y' + o)zX — wiz)^ 

+ (2' + o)iy — C02a:)2] ; 
expanding and simpUfying we have 

2T= 2r+2G+2K, 
where 

T = ^SmC^r'' + / + 2''), 

G = |Sm[(co22 — a?32/)^ + {o)zX — coiz)^ + (coiy — co2a:)2], 

K = Xm[x\o)2Z — a?32/) + 2/'(a?3ic — coiz) + z'(wiy — (a2x)]. 

These expressions can be interpreted as follows: 

T is simply the kinetic energy of the system in its motion 
in relation to moving axes, x, y, z. 

G is the kinetic energy of the system in its rotation about 
the instantaneous axis co (see Euler's formulae and fig. 11) 
and is = ilfA:^(w^/2), where Mk^ is the moment of inertia 
about the instantaneous axis. Finally 

V can be transformed as follows: 

V = o)iEm{yz^ — zy') + 032^m{zx' — xz') + cosZm{xy' — yx'). 
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But we have seen that 

S7w(y2' — zy*) = ^«i, 
^m{zx' — xz') = Bo)t, 
7^m{xy' — yx') = Cws, 

are the projections of the impulse axis (see fig. 12) P, or the 
total moment of mx)m£ntum; therefore, in view of the formula 
of anal, geometry to which reference has just been made, we 
see that 

In other words, V is the product of the impulse axis, the 
instantaneous angular velocity and the cosine of the angle 
between them. 

These separate items T, C and V are not difficult to find 
in most practical problems, as will appear from the following 
examples: 

Example 1. A particle of unit mass is constrained to move 
in a vertical plane B (fig. 45) which has a given rotation, co, 
about a vertical axis Z; the only external force is gravity. 

Z alsoZ 




Fig. 45. 



Let X, Y, Z, be the fixed axes; we shall take, as moving 
axes, X, y, z, where 2 can be directed upward, Xy perpendicular 
to the plane, jB, at any point, A ; and y parallel to the plane. 
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The auxiliary axes, Xq, Yq, Zq, would be drawn through 0', 
parallel to X and 7, but they are not shown. Let 0' be any- 
where on the axis Z and let 0-A = I. In applying Gilbert's 
method we can see that the origin 0' not only has no trans- 
latory motion in regard to fixed axes, but no motion what- 
ever, relative to the latter. Therefore the case is of pure 
rotation about a fixed point 0'; we have no function X to 
calculate. The potential function is, as usual, U = —' gz. 
Now the kinetic energy consists of three terms: 

2T = y'^ + z'^ (this is the kinetic energy as if the axes 
X, y, z were absolute). 

2G = (P + y^)o) (this is the same as would he T = ^(mr^)aj, 
or 7(co^/2) ; which is the kinetic energy of rotation of a particle 
about the axis Z). 

V = o)ly' (this will be evident; ca is the instantaneous velo- 
city, which, in this case happens to be permanent; ly' is the 
same as mvl, that is the moment of momentum of the particle 
(in case of a system we would have had the impulse axis or 
the resultant moment of momentum); the angle (co, P) is 
clearly = 0, since both axes coincide). 

The total kinetic energy will therefore be 

2r = / + z'' + (P + 2/2)0,2 + 2o)ly\ 
Finally the two Lagrange's equations, for z and y, will be 

2" = - g; 

and 

|(2/'+M- 2/0)^ = 0; 

that is 2/" — o)^y = 0. 

From the first equation we have 

z= - Y+Ci< + lC2; 

from the second, y = Cae"' + C4e~"'; and these are the final 
equations of motion. 
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Example 2. Foucault'a gyroscope. A spinning gyroscope 
. 46) 13 so arranged in its gimbals thai its axis of spin can 




move only in a korizonial pUme. The balance is suck that the 
center of gravity of the wheel is exactly at the center of suspension, 
that is at the intersection of Z with the axis of spin; find the equo' 
tions of motion under the action of the rotation of the earth, to. 

Let X, y, z, be the principal axes of the gyroscope; z is the 
axis of spin, and the moment of inertia about it is = C, as 
usual; ^ = 5 are moments of inertia about x and y. Through 
the center, 0, let us draw a set of fixed axes, X, Y, Z, so that 
X and Y are in the horizontal plane, in which the axis z is 
constrained to remain; furthermore, let the axis X be chosen 
to coincide with the projection, on that plane, of the axis w 
of rotation of the earth ; this is clearly indicated on the sketches. 
The angle a will evidently be the latitude. The Euler's angles 
will be: 

4', giving the intersection X\ — oi the plane of the disk 
with the horizontal plane X — Y. 



/ 
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0, the angle between ZO and zO; this is constant and = 90° 
throughout the motion. 

<p, the angle swept by any line xO, in the plane of the disk, 
in a given time; so that dcp/dt is the velocity of spin. 

Let us calculate the three parts of the kinetic energy, that 
is T, C and V: 

T (that is considering the relative motion only). From 
the motion of a body with one point fixed we know that 

2T= ^a?i2 + ficoz" + Ccoa^; 
or, in our case (A = B), 

2T= A(a)i^ + 0)2^) + Co)i\ 

but we have the following transformation formulae (see fig. 26) 

dxl/ . ^ . dO 

coi = -^ sin ^\TL (p -\- -rr cos (p, 

d\l/ . dd , 

C02 = "tT" sm cos ^ — ^ sm <p, 

Considering that d = 90° we have 

0)1 = ^' sin (p; C02 = ^' cos (p; C03 = ^'; 

whereby 27= ^^'' + C<p'\ 

C is the kinetic energy /(co^/2), of the rotation of the 
system x, y, z (including the gyroscope) about the instantane- 
ous axis of rotation, which is parallel to the rotation of the 
earth (mention has been made before, why it is possible to 
transfer rotations from one point to another). But how can 
the moment of inertia / be calculated in reference to this 
axis CO ? The ellipsoid of inertia with regard to principal axes is 

Axi^ + Byi^ + Czi^ = 1 



156 Lagrange's Equations. 

(where Xi, j/i, Zi, are the current coordinates of the ellipsoid; 
see Bowser, Anal. Mech., pp. 443 and 444); or, in our case 
A(xi^ + y^) + Cz^ = 1. We also know that the moment of 
inertia about any line passing through the center of ellipsoid 
(for instance the line co) will be 

I = A cos^ a-\- B cos^ j8 + C cos^ 7; 

where a, /S, 7 are the angles of the line with the principal axes. 
In our case I = ^(cos^ a + cos^ /S) + (7 cos^7; now 7 is the 
angle between co and z, and can be immediately shown to be 
cos 7 = cos a sin^ (by projecting co upon z). Then, again, 
from the well-known formula of Anal. Geometry [see Bowser, 
Anal. Geom., p. 261, form (2)] we have 

cos^ a + cos^ /3 + cos^ 7=1; 

so that 

cos^ a + cos^ j8 = 1 — cos^ 7. 

Substituting these expressions in our formula for / we have 

7 = ^(1 — cos^ a sin^ ^) + C cos^a sin^ ^; 

hence 

2G = o)^[A + (C- A) cos2 a sin^ <p]. 

V; this is the product of the instantaneous velocity of rota- 
tion (in our case, co, the permanent rotation of the earth) by 
the impulse axis, and by the cosine of the angle between them : 

K= coX PX cos (co, P). 

The projections of the impulse axis on x, y, and z are (since 
A = B) Ao)u Ao)2 and Ccoa; the latter is (see transformation 
formulae) = C<p^; the other two components can be con- 
veniently projected upon OXi (in the plane x — y), which 
gives -4(coi cos ^ — C02 sin (p), which, in view of transformation 
formulae, = 0; the same two components can also be pro- 
jected upon 0-Z, which gives -^(coi sin ^ + C02 cos (p), which 
in view of the same transformation formulae, = Ax//'. It 
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appears therefore that P is the resultant of A\l/' (along Z) and 
of C(p', (along z) ; hence its projection upon co gives 

V = o){A\l/' sin a + C(p' cos a sin ^) ; 

therefore the total kinetic energy will be 

2T= Ayl/'"- + (7^'' + a?2(^ + (C - ^) cos^ a sin^ ^) 

+ 2(j){A\l/' sin a + C'^' cos a sin ^). 

There are no external forces, since U and -Sl are both = 0, 
because of the suspension at the center of gravity. 
Hence the two Lagrange's equations: 






d^fdT 
dt 



(1) 



The first of these gives 

•T. (C(p^ + o)C cos a sin ^) = 

or 

^' + CO COS a sin ^ = k, (2) 

Instead of expanding the second equation (1), the following 
artificial method will be used: The first equation (1) will be 
multiplied by ^', the second by \[/'; adding the products we 
have 

d / ,dT , ,,ar\ d f ,dr , ^,dT\ ^ 
dA^e;^^+^W-dt[^^d^^ 

The second term is simply = dT/dt; hence, integrating, 
Calculating separately dTjdip' and dTjd^'y multiplying them 
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by ^' and ^' and adding, we have 

A\l/'^ + Ao) sin a-^' + C<p'^ + Co) sin a sin ^•^' = T + A- 

Multiplying by 2 and substituting 2r as originally calculated, 
we have after reduction 

^^'' + C<p'^ - w(A + (C- A) cos2 a sin^ ^) = A. 

Now CO can be calculated from the fact that one revolution of 
the earth takes 24 hours or 86164 seconds; hence, 

2t 



0) = 



86164 



= .000073 



and the square, w^, of this will be insignificant; therefore, 
dropping the third term, we have 

A^^ + Cip'^ = h; (3) 

substituting ^' from (2) and again dropping the term in- 
volving w^, we have 

Axl/'^ — 2Cko) cos a sin ^ = / 

(/ being a new constant, absorbing Ic^). Here, assuming the 
angle zOX = X, that is putting 



we have 



^ = |+x 



.,2 2Cku) ^ f 

A -J- COS a cos A = "7 > 



which means pendular motion with a period 



= 2^ J 



Cko) cos a 
(the reader, no doubt, knows that an equation 



^j = .(cos^-a) 
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means pendular motion equivalent to a simple pendulum of 
length I = 2g/s). 

Therefore the axis of the gyroscope will oscillate about the 
neutral position 0-X; when the oscillations die out, due to 
frictional and other resistances, the axis of spin will place 
itself in the meridian plane, and mil indicate true north; this 
is the principle of the gyro-compass. 

On the other hand we were not limited in location of the 
plane XOY (of the motion of the axis of spin), which could 
have been directed in any manner. The calculations would 
have been the same, and the neutral position of the axis would 
always tend to coincide with the direction of the axis of the 
earth. Thus, if the axis of the gyroscope were free to move in 
a meridian plane, it would finally set itself parallel to the axis 
of the earth; its angle with the vertical would then indicate 
the latitude of the place. 



CHAPTER V. 

Small Oscillations. 

I. General Remarks. — The first mention of vibratory mo- 
tion was made in elementary kinematics (Bowser, Anal. 
Mech., p. 276) in connection with the problem of a particle, 
revolving in a circle, uniformly; it was shown that the linear 
acceleration of the particle, along a diameter, varies directly 
as the distance of the projection of the particle from the 
center. Indeed, if w is the constant angular velocity and a 
the radius of the circle, we have x =^ a cos 6, where x is the 
distance of the projection of the particle from the center 
and 6 is the angle of the diameter with the radius drawn 
from the center to the particle. Hence 

dx . z,^^ 7 . /» 

TT = ■" a sm ^-TT = — aw sm ^: 
at at 

also 

JPx de , 

t;^ = ^ ao) cos 6 Tr= ^ <trx, 
dr at 

that is the acceleration is proportional to the distance ar, 
which is exactly what we tried to show. Note the negative 
sign, which shows that the acceleration is always directed 
toward the center. Later, at the very beginning of Dynamics 
(Bowser, Anal. Mech., p. 295) we had a somewhat similar 
problem: find the motion of a particle under action of an 
attractive force, varying as the distance of the particle from 
the center of force. Assuming the mass to be = 1, and that 
the coefficient of proportionality of the force to the distance 
is IX (that is, at the distance x the attraction is = /ix), we have 
(force = mass times acceleration) 

^ = - MX, (1) 

160 
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the negative sign meaning that the force opposes the motion, 
in other words showing the tendency of the force to diminish x; 
Multiplying by 2dx and integrating we have 

- (fy=M(a^-a^). (2) 

If x = a is the extreme travel of the particle, the distance 
fr om which the particle started when t was = 0, then — (dxl 
Va?" — a:^) = l//irf<, the negative sign being selected because 
X decreases as t increases. Integrating from t = Q to t = i 
we have cos~^ (a;/a) = V^'t or 

t = --=-cos-' -, X ■•■ -^ - -'^(3) 

/ ■ 
so that when x ^ a, ^ = 0; also when .r = 0, t = 7r/2i//x, 
which is the time of one travel from extreme position to the 
center. Having passed the center, where the acceleration is 
= (from (1)) and where the velocity is maximum and 
= aV/i (from (2)), the particle will travel the same distance 
on the other side, until the velocity is again = and the 
acceleration is again maximum and equal — /xa. So that the 
particle will oscillate between the two positions + a and — a; 
the time required for the particle to travel to the other extreme 
and back is evidently^four times the value just found, t= 7r/2 V^/i. 
It is called complete period of oscillation, T, so that 

r-2,i, (4, 

or, since from (1) fi is equal to the acceleration divided by 
the extreme deflection, we have 



^ lextreme deflection 

A acceleration 

The form (5) should not mislead the reader: it is exactly 
equivalent to (4) and therefore shows that the period is 
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independent of the extreme deflection taken by itself, but 
depends only upon the value of the constant /x. A super- 
ficial inspection of (5) might suggest a certain doubt as to 
this being true, until the matter has been thought over. But 
the form (5) is useful in applications of the following nature: 
1. A weight of 10 lbs. rests on a smooth horizontal plane, 
between two pins, to which it is connected by springs (Fig. 
47). The springs are such that a displacement of 1 inch either 
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Fig. 47. 

way calls out an elastic resistance of 2 lbs. ; that of 2" causes 
the resistance of 4 lbs., etc. Find the complete period of 
oscillation of the weight, if pulled to one side and let go. 
We know that the period will be independent of the displace- 
ment (alone), but to check this up we shall calculate the 
period for, say, a deflection of 1" and that of 2". If the 
displacement is 1", the force called into play is 2 lbs. Now 
the mass of 10 lbs. is 10/32.2 and therefore the acceleration 
corresponding to the deflection of 1" will be = force/mass 
= (2 X 32.2)/10; the deflection is = 1" or 1/12 ft.; so that 
the complete period will be 



r = 27r ^1 



1 X 10 

or .73 sec. 



12 X 2 X 32.2 

Now if the deflection is 2", we have an identical result: here 
the deflection is 2/12 or 1/6 ft. The acceleration is = force/ 
mass that is (4 X 32.2) /lO; so that the period 



\6X 



^ ^^ - .73 sec. 



6 X 4 X 32.2 



2. A 5-lb. weight hangs on a light spring; each additional 
pound stretches the spring one-half an inch. Find the com- 
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plete period of oscillation of the 5-lb. weight. Here we can 
use either the form (5) or the form (4). The latter can be 
applied as follows: to find /x we only have to remember that 
at the distance x the attraction is jluc; for instance, if under 
5-lb. load the stretch is x ft., then under 6 lbs. it will be 
a:+ (1/24) ft.; or 

5 = [iXy 



6 = M(a: + ^), 



whence /i = 24. Hence the period T = 27r/l/24 = 1.28 sec. 
In practical applications it is often required to find n, the 
number of complete oscillations per second, or N, the number 
of complete oscillations per minute. It will be clearly seen 
that 

1 1 , XT 60 

n = -fp, also that N = -jp , 

n is often called frequency of oscillation. 

In the above problems and in the equation (1) we have 
shown that, if the force is exactly proportional to the dis- 
placement, the period is independent of the latter and the 
main problem of the theory of oscillations, which is to find 
the period, is then comparatively simple. All we have to 
remember (and this should be firmly committed to memory) 
is that 

dJ'x 



and 



^--'" 



(t)-=M.'-^) 



mean exactly the same thing, namely simple harmonic motion 
X = a cos Viu< (from (3)) with the period T = 2t! Vju. 

The reader must learn how instantly to write the period of 
any such oscillation as (a^ + 6^)^" + cg(p = 0, which is 
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r = 27r V(a^ + V^)lcg\ and in general it would be well for 
him to revise the usual methods of solution of the following 
important differential equation 

a -Tig" + feic = (same as (1)) ; 

he would then see that the general solution of this equation is 
a; = ^ sin V(fc/a)< + B cos ^(b/a)t, where A and B are 
constants. He would also see that this form can be easily 
transformed into a; = if sin ( ^[(bjd)t + N), where M is called 
amplitude of motion and N is a characteristic known as phase; 
so that the solution in either form means pendular or sinitsoidal 
or simple harmonic motion. 

How tremendously was increased the diflBculty of the prob- 
lem when the force is not exactly proportional to the dis- 
placement was shown later, in the problem of the simple 
pendulum (Bowser, Anal. Mech., p. 351). We found that 
the period of such a motion cannot be found at all by integra- 
tion. The expression of dt can be reduced to Elliptic Integrals, 
which form the subject of a special branch of mathematics, 
highly interesting and most useful in applications. 

However, dealing with the problem of the pendulum, we 
limited ourselves to the special case where the oscillations 
considered were very small (4° or 5° each way), which made 
it possible to introduce certain simplifying assumptions, to 
reject certain quantities, as insignificant, under our limitation, 
etc. In view of these assumptions we finally derived the 
well-known formula T = 2t ^(l/g), where T is the period 
of a complete (double) oscillation of the pendulum, and I is 
its length (in feet), g being the acceleration of gravity (= 32.2 
ft. per sec.^). For greater amplitudes of swinging this formula 
is, of course, absolutely worthless: the reader will readily 
imagine that the period can be anything up to (and including) 
infinity (for the case where the pendulum is allowed to swing 
180° each way). The accompanying curve (Fig. 48) shows 
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the rapid increase of the coeflScient before the square root, 
corresponding to increasing amplitudes. It should be re- 
marked here that in problems on oscillatory motion it is often 
customary to express the solution not as period of oscillation, 




V 60" 120° _ 180* 

Fig. 48. 

Ty but as length of equivalent simple pendulum, L The period 
always comes out to be some such expression as T = 2t Vjf ; 
then all we have to do is to equate this to T = 27r ^(hg), and 
to solve for /, which, then, is the required length of equivalent 
simple pendulum. 

Now the problem of the simple pendulum was the first ex- 
ample of the exceedingly large and important class of problems 
known as small oscillations. Here we purposely limit ourselves 
to the small amplitudes, so that all arcs, velocities and dis- 
placements will be necessarily very small; their squares or 
products will be negligible; the sines of small angles will be 
replaced by the arcs proper, etc. 

We shall illustrate the principle of small oscillations on the 
following examples: 

Example L — A homogeneous hemisphere performs small 
oscillations on a rough (no sliding) horizontal plane; find the 
motion (Fig. 49). Let G be the center of gravity of the 
hemisphere, CG = c, and the radius = r; the angle NCA 
will be denoted d. Since N is the instantaneous center of 
rotation we can take moments (due to rectifying action of 
12 



> 
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gravity) about it. Let h be the radius of gyration of the 
hemisphere about its center of gravity G. Then the radius 
of gyration about N will be calculated by taking into account 




the square of the distance GN (see Radiics of Gyration, Bowser, 
Anal. Mech., p. 435). Remembering now that moment of 
inertia X angular acceleration = acting momentf we can write 

m(Jfc2 + (?iV2)tf" = - mgc sin 0. 

Now this equation is absolutely exact; but since we limited 
ourselves to small oscillations we can simplify it by assuming 
thai GN is practically = r ^ c; also that sin = 0, the arc 
being very small. Hence 

[A:2 + (r - cyW + gc0 = 0, 
so that the period is 



= 27r>J' 



Ic'+jr- cf 



Inspection will show that A:^ + c^ = square of the radius of 
gyration about C, which can be shown to be (2/5)r^; also c 
can be found = (3/8) r, so that finally T = 27r ^l{26/15)(r/g); 
this result will not hold good for any but small deflections 
from neutral position. 

Example 2. — The theory of the Dynamic Balancing Ma- 
chine, originated by the author, involves the following problem 
(Fig. 50): a heavy loaded beam, whose one end is hinged at 
A and the other is supported by the spring S, can perform 
small oscillations in a vertical plane. Find the period. The 
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mass of the beam is M, its moment of inertia about A is Mk^, 
k being the radius of gyration, and the distance of its center 
of gravity f rom \4 is b. Let the initial deflection of the spring 
under its load be = 8. The system is slightly disturbed 




Fig. 50. 

through a small amount x (and a small angular deflection <p) 
and let go. The force W acting on the spring when deflected 
through S will result in a moment Wl = Mgb about A; when 
deflected through another small amount x, this will increase 
by the amount 



which is then the acting moment, opposing the deflection, and 
can be written = — mg(bx/d) ; now, since the acting moment 
= moment of inertia X angular acceleration, we have 

- Mgj = MPip'\ ^ ' 

Here the simplifying hypothesis is that, since the oscillations 
are small, we can assume x = <pl, whence ^" = x''/l; substi- 
tuting into the result just found we have 



X 



Mgj=M]^f, 



-Vv^ 



V X 



K '1 ! 



so that a:" + {gblim)x=_0, whence T = 27r <{migbl) ; writ- 
ing this as T = 2Tr ^S/g) • '^Qc^jbl) we can see that if we 
simply assumed the spring to vibrate under its own load, and 
due to its initial deflection 8, the result would be only. 
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r= 2Tr'4bJg (see similar problem above). But since the 
beam oscillates about A, this brings into play an important 
coefficient '^P/bl, depending upon three other factors, char- 
acterizing the system. Since the beam in question is very 
heavy and the deflection 5 is always very small, the oscillation 
cannot very well be any other but exceedingly small, so that, 
contrary to example 1, we do not have to emphasize the im- 
posed limitation, the conditions of the problem take care of 
this automatically. 

Later on these two easy problems will be solved by applying 
a different method. 

In general, therefore, the problem of small oscijlations 
involves a certain amount of experience in introducing simpli- 
fying assumptions. 

We shall now show that by applying Lagrange's equations 
we can make this work much more uniform, almost mechanical. 

Only oscillations about the position of equilibrium will be 
considered here; while in more advanced book is likewise 
included the problem of oscillation about steady motion. 

n. Two important principles of equilibrium will now be 
established. 

(a) When we have equilibrium, the potential function U is 




Fig. 51. 



constant; or, in other words, partial derivatives of U with 
reference to any small displacement from the position of 
equilibrium are = 0. This will be readily understood: from 
the integral of kinetic energy we have T = U + h, where T 
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is the kinetic energy, U the potential function, and A is a 
constant (arbitrary). Therefore, in equilibrium, where the 
kinetic energy is = 0, the function U will naturally be con- 
stant. Thus for a pendulum (Fig. 51) in general U = mgs 
cos <p, which for the position of equilibrium is simply U = mgs, 
A practical way of looking at it is that when the pendulum 
is at rest, then no infinitesimal lateral displacement will 
consume or produce any work at the expense of the force of 
gravity: all such infinitesimal displacements can be only 
horizontal, say Oa or Ob, and therefore the function U will 
remain constant. 

(b) The second important principle, established by Lagrange 
is that for a position of stable equilibrium the potential func- 
tion U must be not only constant, but at the same time 
a maximum. There exists a rigid proof of this proposition 
but we shall omit it here and the working of this principle 
will be explained by reference to the following typical problem: 

A particle (Fig. 52) whose weight is M is attracted toward 
the centers A and B by means of springs AM and BM whose 




Fig. 52. 

constants are X and fi (that is the force of attraction is X X -41/ 
and fjL X BM respectively) ; the particle is constrained to 
move about a fixed center on radius c, being midway 
between A and B. Find the position of equilibrium, also 
establish whether such equilibrium is stable or unstable? 

Assuming <p to be the angle between AO and the present 
position OM, let us form the expression of virtual work of 
such a system for any small displacement of M in its rota- 
tion about 0, say moving downward. We have already found 
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the expression of elementary work dW for this system, starting 
from the following (page 26) 

dW= - \AMdAM - ixBMdBM + Md{c sin ^); 

here the last term of the second part is the virtual work due 
to gravity; the sign — before the first two terms means, same 
as before, that the attraction opposes the increase of the 
distance from the attractive center; if there is any ambiguity 
as to whether this is increase or decrease, the geometric con- 
ditions will at once put this straight. Finally the following 
was found 

dW = c[a(ji — X) sin ^ + if cos (p\d(p = dU. 

To find the maximum (or minimum) value of U we equate 
to this expression of dU which gives 

M 

tan <p = 



a(X — jji) ' 
To find whether this is maximum or minimum we differentiate 



once more 



d^W = c[a(p — X) cos ^ — ilf sin (p]d(p^ 
= c[a(jjL — X) cot <p — M] sin <pd(p^; 

considering the value of tan <p just found, the expression in 
square brackets [ J will be always negative. 



a^ 



so that the value of sin (p will decide the condition of stability: 
if the value of <p, found from 

M' 

tan (p = —^ r , 

a(X — /x) 

is < TT, the second derivative will be negative and there will 
be stable equilibrium. Otherwise, when ^ is > tt, the equi- 
librium is unstable. 



Oscillations. 171 

We shall now somewhat enlarge our conception of potential 
function U. This has been defined as function of coordinates, 
such, that its derivatives with respect to coordinates give the 
value of corresponding forces, as for instance dU/dx = X or 
dU/d<p = $, etc. But the function U may have contained a 
certain constant, C, so that a more general way would have 
been to write U + C, where C may be any constant, and of 
course drops out in differentiation. For instance for a pendu- 
lum we simply had U = vigs cos <p; this really should have 
been written U = mgs cos (p -{- C (where C is a perfectly 
arbitrary constant, positive, negative or zero); this was not 
done simply because we were always concerned only with 
derivatives of L^ in which C does not enter; in other words 
we heretofore implied the assumption C = 0. But hereafter 
we shall assign a definite value to C: we shall select it equal to 
minus the viaximum value of U, so that Umax + C = Umax 
— Umax = 0; in other words we shall select C so that the 
maximum value of U is =0; and therefore all other values 
of U, viewed in this light, are negative. Consider for instance, 
an ordinary physical pendulum (Fig. 51). Here C is chosen 
= — mgs, so that the maximum value of U is zero; the in- 
crease of (p means decrease of the value of (U + C), that is 
mgs cos (p — mgs, until finally the value ^ = tt is reached, in 
which case Umin = ~ 2mgs. This also means equilibrium, 
but it is unstable. For stable equilibrium the value of U must 
be maximum; and we have chosen this maximum value = 0; 
of course, so far as stability itself is concerned, any value of C 
would answer, but the object of selecting C so as to have U 
negative, or at the most = 0, will appear presently. 

The physical pendulum is a system with one degree of 
freedom; as an example of a system with two degrees of 
freedom we shall consider a double pendulum (Fig. 53), for 
which the expression of U is mg{l cos ^ + 6 cos <p) + Mgs X 
cos (p -\- C (check this by considering the elevations of the 
center of gravity of each rod). Here, again, the constant C 
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will be chosen so that in the lowest position U = 0. This 
will be the stable equilibrium and it will be seen that any other 
position, corresponding to any change of ^ or ^ will mean 
negative U, so that here again the condition of stable equi- 
librium means maximum U; when (p and ^ are both = tt, we 
have unstable equilibrium, so that the slightest displacement 




Fig. 53. 

will upset such a system, but will not result in its coming 
back to the former position, as would be the case for stable 
equilibrium, where U is maximum (even if this maximum is 
but = 0). Limiting ourselves to the simplest problems, 
involving permanent constraints (so that the kinetic energy 
is merely a quadratic function of the generalized velocities q' 
etc. ; see (4') and (4") Chapter II) we shall study the effect of 
very small displacements, such as q from the position of equi- 
librium. Only systems with one and with two degrees of 
freedom will be here considered. 

m. Systems with one degree of freedom (for instance a 
pendulum). We shall consider the effect of giving such a 
system a slight displacement, g, from its position of equi- 
librium. Since, according to (4'), Chapter II, the expression 
of kinetic energy can be thrown into the following shape 



2T=^{q)Xq'\ 



(4') Chapter II 



where $(g) is a function of q only, and constants (but does 
not involve q^), we can readily see that, (1), The kinetic energy 
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itself is very slight, since all velocities in S(mD^/2) are very 
small; also, (2), the displacement q is small, in accordance with 
the very nature of our problem; therefore, (3), the generalized 
velocity </' must likewise be small. 

Now suppose that the function $(g) is developed by Mac- 
laurin's formula (Bowser, Calculus, p. 83), 

2/ = /(a:) = /(O) + xf' (0) + •.., 
so that 

*(g) = *(0) + g*'(0)+ •.-. 

Now we shall neglect the term g$'(0) and will keep only the 
first term 0(0), so that instead of the above expression 

2r=$(g)Xg'' 
we shall have simply 

2r= $(0)X/ 

(dropping higher terms). Here $(0) is function of q no longer 
and this is precisely what we tried to accomplish. It is not 
always necessary to actually develop $(g) since the same 
result can be obtained by inspection; but the main object 
must be constantly kept in mind: the expression of T must 
be thrown into a form that does not involve q but only q\ 
The above expression 2T = $(0)g'^ is generally written 
simply as 

2T=aq'\ (1) 

This is our first equation. The second can be obtained by 
developing the expression of U (which as we know is a func- 
tion of q only) by Maclaurin's theorem 

U = F{q) = F(0) + qF^ + ^F'(O) + • • -. 

In this series F{0) = because we agreed to adjust the con- 
stant C so that U is always = in equilibrium (that is when 
g is = 0); likewise the term qF\0) is = 0, because in the posi- 
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tion of equilibrium (that is when g = 0) the value of CT is 
maximum^ so that the first derivative vanishes and the 
second derivative is negative. Finally, dropping higher terms, 
we have U in the following form 

U=-'4, (2) 

where 6 is properly determined from (g^/2)F"(0) = — (b^/2). 
In other words we have outlined the possibility of expressing 
the potential function as constant times the sqiuire of the inde- 
pendent coordinate q. If this can be written down by inspec- 
tion, the reader should by all means do so. He must satisfy 
himself that, in actual practice, if the deflections are small, 
and reckoned from the position of equilibrium as zero, then 
U will always come out to be a homogeneous quadratic func- 
tion of q. Take for instance a weight on an elastic string: 
let the mass be m and the constant of the string = X (that is 
a stretch of z means resistance of Xa; lbs.). The potential 
function is found by computing the expression of elementary 
work 

dW = — \xdx + mgdxy 
so that 

[7 = - -y + mgx + C, 

which is not a quadratic function; but in position of equi- 
librium we have evidently \x = mg. To avoid confusion 
let this value of x be designated 8, so that X8 = m^ or 5 = mgfK. 
In other words what we had designated as z will now be made 
up of two parts: x = d + z, where z is the varying coordinate 
and 8 simply characterizes the position of rest: x = S when 
z is = 0. 

Now substituting the value of x just found into the above 
expression of U (and remembering that X = mg/d) we have 

U= -^(8^ + 2^) + C. 
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We have agreed once for all to select the constant in such 
a way that [7 is = when 2=0 (position of rest). So that 
finally 

which is a quadratic function of the independent coordinate. 
This happens every time we limit ourselves to small dis- 
placements and reckon them from the position of rest. To 
return to the equations (1) and (2): applying Lagrange's 
method, instead of 

dt dq' dq ~' dq' 
we have simply 

ag" = - &g, (3) 

which is the required differential equation of the small oscil- 
lation in question. The reader will see at once that it is the 
same as the well-known differential equation 

of which the integral is y =' A cos (rx + B) ; or, in our case, 

q = a cos {rt+ j8), (4) 

where r^ = b/a and a and jS are constants to be determined 
from initial conditions. As a rule we are not interested in 
them as much as we are in the constant r. The equation (4) 
means periodic motion with the period T = 27r/r. Indeed, 
to the value of t we can add any integer, n, multiplied by 
27r/r and the result will be the same: thus 

g = a cos (r^ + /3) = a cos r ( 1- ^ )+ P 



= a cos 



= a cos 



['(-'•r+O+^J 

Mn-y^+n + ^J = etc. 
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in other words q passes through the same value every 27r/r 
seconds. Now, finding such periods of oscillation is precisely 
our problem. We shall always mean double or complete periods. 

So that, in order to apply Lagrange's method to the problem 
of small oscillations we have to (1) express T as constant times 
g'^; and (2) express U as constant times q^; after which the 
usual application of the method leads to an equation from 
which m can be seen at a glance, thus solving the problem. 

By way of illustration we shall check the well-known formula 
of the simple pendulum, T = 2t '>ll/g (Bowser, Anal. Mech., 
p. 353). We shall suppose that the length of the pendulum 
is I and its mass = m. Let the small deflection from the 
position of rest be <p, so that the current angular velocity 
will be <p\ This small deflection <p will be our generalized 
coordinate (same as q) ; and the kinetic energy will be 

2T= mv^= m{l<py = mP<p'\ 

On the other hand the potential function U will be U = mgz 
+ C = — mgl{l — cos <p), or, developing cos <p by Maclaurin's 
formula (Bowser, Calculus, p. 86) and dropping higher terms 
we have cos ^ = 1 — i<p^/2); so that 11 = — mgl{<p^l2). There- 
fore our only Lagrange's equation will be 



dt 



{mP<p'^) = qz(- ^9l^) or ^". = - ^ <p, 



whence <p = a cos ( ^{g/l)t + j8), — pendular motion of which 
m is = V^ so that the period will he T = 27r/m = 2t '^ijg, 
which is the required formula. ^ 

For ready reference we shall give a few expansions, well 
known from calculus. 

In general 

y = m = /(O) + f(0)x + /"(O) -+■■•, 
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so that 

.3 

6 



sin a: = a: — ^ + • • • (for small values of a: sin a; = x), 



cos a: = 1 — X- + • • • (often simply taken = 1), 

x^ 
tan a: = a: + «- + • • • (for small values of x tan x = ar), 

x^ 
sec a: = 1 + ^ + • • • (often simply taken = 1). 

Example L — Rocking hemisphere, Fig. 49. The height of 
the center of gravity G above the plane is z = r — c cos 0; so 
that the potential function will be t7 = — mgz + C = mg X 
(r— c cos d) + C; here C is selected so that the maximum value 
of U (when 0=0) is = 0; in other words the constant C is 
made = mg(r — c) and finally U = m^c(cos 0—1) or, in 
view of the expanded formula of cos we have 

t/ = — mgc 2" . 

This corresponds to the equation (2) in the general theory: 
[7 = — bq^/2. So far as the kinetic energy T is concerned we 
have T = IUoy2) ot2T= 10'^ = m(P + (r - 0^)0'^ (see how 
the moment of inertia was found before). This corresponds 
to the equation (1) of the general theory: 2T = oxf^. Hence, 
applying Lagrange's equation, or simply by (3) we have 
a^" = — 6g or m(P + (^ — cf)0'' = — mgc0, the same equa- 
tion as we had before. 

In this and the following easy examples the answer could 
have been written almost without any calculations, which 
are given here only for the sake of making matters absolutely 
clear for the beginner. 

Example 2, — A weight mg is suspended on a vertical coil 
spring of which the characteristic is X (that is an elongation x 
calls into play a resistance of \x lbs.). Find the period of 
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the small oscillation of the weight when disturbed from its 
position of rest and let go. 

Supposing that the small deflection from position of rest is z 
we have the following expression of kinetic energy 

2T=mz'\ 
The potential function, as was shown above (see page 175), is 

(where 8 is the initial stretch due to load itself). 
So that instead of 

±dT__dT_dU 
dt dz^ dz dz 

we have mz^' = — mgz/S which means that the period will be 
T = 2TrVilg, which we had before. 

Example 3. — Oscillating beam. Fig. 50. Here we shall take 
as independent coordinate the angle <p, characterizing the 
deflection from horizontal position. Therefore the generalized 
velocity will be ^'; and if the moment of inertia is = Mlc^ 
then the expression of kinetic energy will be 27" = MW(p' . 
Now, the expression of elementary work will be 

dW ^ -\XdX^Mg^jdX 
or 

where X is the total travel of the beam end from the no-load 
position of the spring. 

When the beam is at rest the load on the spring is = M{gh\V) 
and is resisted by X8, where 8 is the initial deflection of the 
spring; so that in the first place X = Mgb/lS; and then X=S+ x 
where x is the characteristic small deflection of the end 



> 
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of the beam from its neutral position, 8 being merely a constant 
quantity. Substituting in the above expression of U we have 

^.M^^ + ^^ + Mf+c 

or, with suitable adjustment of the constant C 

^ 2/8 *^' 

now X = l(p, so that finally U = — (MgbP/2lS)if^, whereby 
Lagrange's equation 

dt d<p' dip d(p ' 
or in our case 

MW = - M^^<p, 
which means periodic motion with the period 



= 2t>J' 



m 

gbl 



or 2TrV {bjg) (k^/bl), which is the result we had before. 

Example 4* — ^ particle of mass m is attached to two points 
A and B by two elastic strings whose constants are X and /x 
and whose natural lengths are L and I. The distance between 
A and B is such that when the particle is in equilibrium the 
respective initial stretch of strings is, respectively, A and 5. 
Find the period of this small oscillation (Fig. 54). 

The kinetic energy will be evidently given by 

2r= mx'\ 

In order to find the potential function U let us find an ex- 
pression of elementary work. If the distance x be taken as 
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independent coordinate, the pull exerted by the right string is 
/i(Z + 5 — a:) ; that exerted by the left string is X(Z + 5 + A 



li.X.A 



-00- 



l.f^.S 






Fig. 64. 
— or). So that the elementary work will be 

dTF = /i(Z + 8 - z)dx - X(Z + 5 + A - x)dx, 
whence 

and the Lagrange's equation will be 

mx" = — x(K -j- ii) + A. 

The solution of this differential equation will consist of a 
trigonometric term plus a certain constant. The reader is 
invited to investigate the nature of this constant. However 
this is quite immaterial for finding the period, which will be, 
according to rules given above. 






The reader will check this result by forming new equations, 
where x has been measured from the position of equilibrium, 
and not from one of the fixed points. 

IV. Systems with two degrees of freedom of which a typi- 
cal illustration is the double pendulum. Fig. 53. Here each 
pendulum has its own period of oscillation, both periods being, 
in general, different from those corresponding to each pendulum 
swinging separately. The general treatment of the problem 
is similar to that of the previous case of one-degree freedom: 
all displacements are small and measured from the position 
of rest as zero. There will be two independent coordinates. 
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say Qi and q2, corresponding to the two-fold freedom of the 
system. The expression of kinetic energy for small dis- 
placements (and therefore small generalized velocities q\ 
and q^) will be given by a homogeneous quadratic function of 
generalized velocities, such as 

2r=agi'' + 26gi(/.+ cg2''. (1) 

Likewise the potential function TJ will be given by a homo- 
geneous quadratic function of the generalized coordinates q\ 
and q2, such as 

t^ = - i(«gi' + 2j3gig2 + 792'). (2) 

This expression cannot very well contain any constants 
because 17 = for gi = and ^2 = 0; neither can it contain 
any terms of first degree in q\ or q^ because for gi = and 
ga = the potential function t7 is a maximum (compare with 
explanation for one-degree freedom). 
Higher terms are neglected in expressions for both T and U, 
Having T and TJ we can write down Lagrange's equations 
for both coordinates q\ and qi 



d dT 
dt dqi' ' 


dT dU 

dqi dqi ' 


d dT 
dtdqi' 


dT dU 

dq2 dq2 ' 



(3) 



(4) 



this will mean in view of (1) and (2) 

aqi' + bqz'' = — {aqi + fiq2), 
bqi" + cq^' = - (/3gi + 7^2). 

The usual method of solution of these linear differential 
equations of second order with constant coeflScients consists 
of assuming 

q\ = Pi cos (r^ + A:); 72 = P2 cos (ji + A:), (5) 

where pi, p2, ^ and h are constants. What we are interested 
13 
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in are the values of r, of which as we shall presently see there 
will be two, fi and r2, and from which we can determine the 
required periods Ti = 2x/ri and T = 27r/r2. 

In order to find the constants ri and r2 let us substitute the 
tentative solutions (5) into Lagrange's equations (4). We 
get 

apir^ + 6p2^ = (xpi + Pp2, (6) 

bpir^ + cp2f^ = jSpi + yp2. 

Eliminating the ratio pilp2 we have an equation of second 
degree in r^ 

(a - ar^Xy - cr^) - (fi - hi^f = 0. 

Among the four roots of this equation, it can be shown that 
two values, r\ and r2 are positive. They are the required 
results. 
Example 1. — Find periods of small oscillations of the double 




Fig. 65. 

pendulum. Fig. 55, consisting of two masses, m and mi, sus- 
pended by inextensible strings a and b. 
Here the x, z coordinates of m and mi are 

X = a sin d; Xi= asind + b sin <p; 
z = a cos d; Zi = a cos 0+6 cos <p; 
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so that 



(i)'+(iy=---''^ 



also 



( — )'+( ^ )' = a^6'' + Wv'' + 2abdy cos {<p - 6) 

and therefore the kinetic energy is given by 

2T= maW^ + mxWff'' + 6V'' + 2a6eV' cos (^ - »)]; 

while the potential function will be found from the foUowing 
expression of elementary work 

dW = — mg sin ddJd — mig sin Odd — mig sin <pdip. 

(This can be written down at once by taking tangential com- 
ponents of the weight of each mass, m and mi, and then by 
allowing the system a small displacement dd, maintaining <p 
constant; also by allowing a small displacement d(pi keeping 
constant; then adding the results.) 
Hence the two Lagrange's equations, 

dt 66' ee " dO' 

dt dip' dip dip ' 
become 

a{m + mi)S" + hmiip" cos {ip — 0) — bnii sin {ip — 0)ip'^ 

+ (m + mi)g sin = 
and 

a cos {ip - 0)0'' + hip" + aO'^ sin (^ - S) + ^f sin ^ = 0. 

So far no simplifying assumption has been introduced and 
these equations are strictly correct, although much too diflScult 
to integrate. 
Limiting our problem to small oscillations we shall replace 
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the sines by corresponding arcs; also all cosines of small angles 
we shall simply put = 1 ; likewise all terms involving squares 
or products of velocities will be neglected. Hence 

a(m + mi)0" + hm\ip" + (m + mi)g6 = 0, 
aO" + hip'' + g<p=Q. 

The first equation can be simplified by means of the second; 
so that finally 

amd" — miQip + (m + m{)gO = 0, 

aO" + h<p" + ^e = 0. 

Assuming the usual tentative solution 

6 = pi cos (jt + k) and ^ = p2 cos {rt + k) 

we finally have 

amr^ — (m + mi)g + p^m^g = 0, 
ar^ + {hr^ — g)^^ = 0. 

From these equations we can find the two required positive 
values of r : ri and r2. Considering the special case where 
a = 6 and m = *mi we have 

ay^ — 2g + p2g = 0, 
ar^ + (ar^ — 5r)p2 = 0, 
so that r^ = (g/a){2 db V2); in other words 

ri2 = 3.414^, r22= .586-, 
a a 

whence the periods 

2w la 

ri = — = .541X27rV-, 

n ^g 

T2 = — = 1.307 X2x\/-. 
^2 ^ ^g 

It will thus be seen that the upper pendulum will swing much 
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quicker and the lower much slower than what would be their 
natural rate To = 2t ^a/g. 

Again, if we have a pendulum of the length 2a, the rate 
will be still slower, T3 = 1.414 X 2x Vo/^. 

The investigation of the pendulum of Fig. 53 would lead 
to rather complicated results; the theory of such a pendulum 
is closely connected with the celebrated problem of Bell and 
Clapper, so ably given by Dr. Slocum in his book The Theory 
and Practice of Mechanics. 

Example 2. — An investigation into the subject of Dynamics 
of the Automobile has brought the author to the following 
elementary problem: the Fig. 56 represents an imaginary 
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carriage, consisting of one tire T supporting a platform, 
representing the so-called " unsprung '' weight mg. The 
" sprung " weight Mg rests on the spring S mounted on 
the platform. If the platform is steady, the natural period 
of oscillation will be 



= 2.yjf, 



where A is the initial deflection of the spring S under its load 
Mg; this we had before. As a matter of fact, however, the 
tire itself is elastic and the platform is susceptible to its own 
up-and down oscillations, although the elastic constant of the 
tire is much greater (the tire being stiflFer) than that of the 
spring, say Xi = 15X, where Xi is the constant of the tire 
and X that of the spring. It will thus be seen that the problem 



y 
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is reduced to a well-known problem of double spring, as shown 
on the right, Fig. 55. 

We shall take as independent coordinates q\ and q<i deflec- 
tions of each mass m and M from its neutral position. Then 
the kinetic energy is evidently given by 27" = mq\ + Mq^^. 

In calculating the potential function U from the expression 
of elementary work we shall first consider the variation of q\ 
alone, quite disregarding q2\ then the variation of q2, although 
allowing for the fact that the variation of qi may have dis- 
placed the zero of q^. For the upper spring s we have 

f^i = y- + mgxi + Ci, 

where, as we had before, x\ is the total stretch, including the 
initial stretch 8 and the deflection q\ itself: X\= 5 + gi. 
Simplifying and selecting proper constant we finally have 

TT ^3 2 

For the lower spring S we have 

It should be observed here that X2 again consists of two terms: 
X2 = A+ q2, where A is the initial stretch of the spring S 
and q2 is the small displacement. But the zero-position itself 
of M is dependent upon that of the support m of the spring S, 
which itself is characterized by its coordinate qi. So that 
instead of X2 we shall substitute 0:2 = A + ^2 — ^i, which 
after reductions and considering proper choice of C2 gives 

C^2 = - § iqi' + q,' - 2gig2) 

and the combined potential function is 

_. mgq,^ Mg . , , ^ . . 

U= - -gg- - 2X (91 + 9' ~ 29i92)» 
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so that the two Lagrange's equations in gi and g2 will be as 
follows: 

mqi" = — Xigi + X(g2 — gi), 

Mq^" = - X(gi - gi). 

Dividing the first equation by m and the second by M we can 
write them down as 

gi" = - agi + 6(^2 - gi), 

(a) 
g2" = — c(g2 — gi). 

These are simultaneous differential equations of second order, 
linear and homogeneous. Their solutions will involve four 
constants. 

It is not unreasonable to suppose that such a motion in 
general will be rather complicated. But let us see if we can 
find such a mode of combined oscillation in which each weight 
will perform simple harmonic vibrations (up and down) of 
the same period and phase, although without restrictions as 
to amplitudes. To test this out let us suppose that the 
solutions can be found in this form 

gi = Vi cos {ri + k)y 

03) 
g2 = P2 cos {rt + k), 



/^ 



where q\ and g2 are the distances from the middle position of 
each weight; pi and 2>2 are the amplitudes; the fact that r is 
the same meaning that we seek that mode of motion where 
both periods will be alike (T = 27r/r); and k being the same 
in both equations meaning that the phase is the same, in 
other words that both weights will be in their respective 
extreme positions, or, say, in their middle positions, at the 
same time. If we differentiate the equations (jS) twice with 
respect to time and substitute in (a), the cosines cancel out 
and we have 

— Pir^ = — api + &(P2 — Pi), 

(7) 
_ pgr^ = — c(p2 — Pi). 
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Dividing by p2 we can eliminate the ratio pi/p2, so that finally 

r^-7^ia+b + c) + ba= 0. 

This is an equation of second degree in r^, which can readily 
be shown to have two roots, ri^ and r2^, both real and positive. 
This means that we have obtained two distinct solutions of 
the form (j8), each solution having its own constants p and k: 

qia = Pi cos (nt + ki)y 
qib = Pi cos (r2< + ki), 
q2a = Pi cos ijit + ki)y 
q2b = P2' cos {r2t + ^2), 

these are the particular solutions of the differential equa- 
tions (a). The general solution will consist of these particular 
integrals, superimposed : 

qi = Pi cos (rit + ki) + px cos {r^t + ^2), 
q2 = P2 cos {rit + ki) + p^' cos (r<ii + ^2), 

where the four constants pi, pi', k\ and ki are arbitrary 
(determined from initial -conditions); while the constants 
Pi and Pi are not independent but connected with p\ and px y 
respectively, by either equation (7), say the second, whereby 

Pi g— r^ pi' c— r^ 

— = and — , = . (f ) 

Pi c Pi c 

It can be shown by theory of equations that this ratio will be 
positive for one and negative for the other of the roots ri 
and Ti. So that gi and qi may have the same sign or may be 
opposite in sign, while being of the same period. The initial 
conditions (that is the manner in which the system is started 
off) can be selected such that one of the constants, say pi is 
= 0; then the other corresponding constant pi is likewise = 
(by f ) ; and we have, instead of the general, more complicated 
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motion, simply 



qi = Pi cos (fit + ki), 
Q2 = P2 cos (fit -\- ki). 



It is easy to construct a model to demonstrate this sort of 
vibratory motion; both weights can be made to move in the 
same sense, or in opposite directions, although with the same 
period; only in the latter case the period will be much quicker 
than in the former. This difference in the periods can easily 
be seen from equations (f)« Indeed, placing pi/p2 = A we 
have (since this may be either + or — ), 

ri^ = c + cA, 
r2^ = c — cA, 

so that the periods will be Ti = 27r/ri and T2 = 27r/r2 respec- 
tively. 

Example 3. — This was likewise suggested by the author's 
investigation into dynamics of the automobile. The car-body 
can have three main degrees of freedom and, accordingly, 
three modes of oscillation: vertical or plunging, angular about 
a longitudinal axis, or rolling; and angular about a transverse 
axis, or pitching. Assuming, however, that we have a two- 
dimensional car, possessing only length and height, let us 
find the periods of its two oscillations. The sketch relative to 
this much smplified problem is given in Fig. 57; a beam of 
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Fig. 57. 



mass m is supported by two identical springs ^/ and B, the 
constant of each spring being X (in other worjois X5 = half 
weight of beam). Such a beam possesses tWo degrees of 
freedom: it can vibrate vertically (independent coordinate z), 
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and can also oscillate about the transverse axis C (independent 
coordinate <p). 

Therefore the expression of kinetic energy will be given by 
27 = mz'^ + ml^<p'^, when k is the radius of gyration about C. 

The potential function U will consist of two parts, one, 
Uu depending upon the gravity; the other, U^y upon the work 
of springs, due to a virtual angular displacement of the 
system. 

Forming expressions of elementary work we have 



hence 



or 



Also 



hence 



or 



so that 



dWi = — 2X(5 + z)dz + mgdz, 
Ui= - MS + zy + mgz + Ci 







dWi = - X(5 + x)dx + X(S - x)dx; 



U2= -'Kx'= -?oV; 





U= -fiz'+aV) 



The Lagrange's equations will be 

and the periods will be T^ = 27r ^[^J2g and T^ = 27r ^{SI2g) • kfa- 

As a rule k is always less than a, therefore the free period of 

plunging is slower than that of pitching. But it is possible 

to artificially increase this value of k as well as to decrease a 
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to the limit which may be imposed by practical considerations, 
and thereby to make the period T^ much slower than is the 
usual case. The object of this is to increase comfort in riding, 
which depends upon the period of pitching as well as upon 
the linear acceleration at the extreme points of the car, A 
and B, which may be even more remote than the distance of 
springs proper from the centerplane C 

If the oscillation is given by ^ = a cos {rt + e) where a is 
the amplitude and r can be calculated from T = 27r/r and is 






i k 

Now, the maximum angular acceleration can be found by de- 
riving (p twice and by making cos (r< + €) = 1; so that 

(Pmaz' = ar^ = a-r -j^ 

and the maximum linear acceleration on the distance a from 
the center-plane will be proportional to a^, hence the further 
importance of reducing a. 

Example 4 (Loney) . — ^A uniform rod (Fig. 58) of length 2a, 




Fig. 58. 

can turn freely about one end, which is fixed on a vertical 
axis. Initially the rod is inclined at an acute angle a to the 
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vertical and is set rotating about the vertical axis with angular 
velocity co. Show that during the motion the rod is always 
inclined to the vertical at an angle which is > a or < a 
according as co^ is greater or less than 3^/4a cos a, and that in 
each case its motion is included between the inclination a and 



cos 



1 ^ [— n + Vl — 2n cos a + n\ 



where n = (aco^ sin^ cc)/Sg. If the rod be disturbed, slightly, 
when revolving steadily at a constant inclination a, show 
that the period of the small oscillation will be 



= 2^4 



4a cos a 



3g{l + 3 cos2 a) ' 

Let, at any time, t, the rod be inclined at the angle 6 to the 
vertical; also let the plane through the rod and the vertical 
be characterized by the angle <p from the initial position. 
Consider an element of the rod, d{, at P, the distance from 
being {. The mass of the whole rod being m, the unit-mass 
will be m/2a; and the mass of the element d^ will be (d{/2a)m. 
The circumferential velocity of P will be evidently = {co sin 
(in its motion about the vertical). The velocity of P in its 
(possible) motion in the plane AOV will be = {&', so that the 
kinetic energy of the element d^ will be given by 

2T = ^m(eO'^^ + ^co^sin^ff) 

and that of the whole rod will be given by 

2r = g (»'' + co2 sin2 6) r ed^ = ^ (»'' + co^ sin^ B). 

Now the potential function U will be easily found from the 
elementary work of the center of gravity of the rod in dis- 
placing it up or down 

JJ = mga cos tf + C 



